U 
Linn? oy TAWA) 


The Deo Jy 


PHILOSOPHICAL 
MAGAZINE 


FIRST PUBLISHED IN 1798 


L. 42 SEVENTH SERIES No. 333 October, 1951 


| A Journal of 
Theoretical Experimental 


and Applied Physics 


EDITOR 
PROFESSOR N. F. MOTT, M.4., D.Sc., F.R.S. 


EDITORIAL BOARD 
SIR LAWRENCE BRAGG, 0.B.E., M.C., M.A., D.Se., F.R.S. 
ALLAN FERGUSON, M.A., D.Se. 
SIR GEORGE THOMSON, M.A., D.Sc., F.R.S. 
PROFESSOR A. M. TYNDALL, C.B.E., D.Se., F.R.S. 


PRICE 12s. Od. 
Annual’Subscription £6 Os. Od. payable in advance 


Neen e ce aaa 
) AND PUBLISHED BY TAYLOR & FRANCIS LTD., RED LION COURT, FLEET ST., LONDON, E.C.4 


A New Publication added to 


THE PHILOSOPHICAL MAGAZINE 


(A Journal of Theoretical, Experimental and Applied Physics) 


QUARTERLY SUPPLEMENT 


ADVANCES 
IN 
PHYSICS 


VOLUME 1 JANUARY 1952 NUMBER 1 


CONTENTS 


Theory of Crystal Growth. By Dr. F. C. FRANK (Bristol University). 


On the Generation of Vacancies by Moving Dislocations. 
By Professor F. Seirz (University of Illinois, U.S.A.). 
Electrical Conduction in Thin Films. 
By Dr. E. H. SONDHEIMER (Cambridge). 


PRICE per annum £2 15s. Od. post free 


Editor: PROFESSOR N. F. MOTT, M.A., D.Sc., F.R.S. 


Editorial Board : 
SIR LAWRENCE BRAGG, O.B.E., M.C., M.A., 


D.Sc., F.R:S. 
ALLAN FERGUSON, M.A., D.Sc. 


SIR GEORGE THOMSON, M.A., D.Sc., F.R.S. 
PROFESSOR A. M. TYNDALL, C.B.E. D.Sc., F.R.S. 


[ 1065 ] 
CIX. Theory of Barium Titanate—Part IT. 


By A. F. Devonsuirge, 
H. H. Wills Physical Laboratory, University of Bristol*. 


[Received June 7, 1951.] 


SUMMARY. 


The phenomenological theory of barium titanate given in an earlier 
paper has been extended, and expressions obtained for the piezoelectric 
constants, elastic coefficients for constant field and dielectric constants 
for constant strain in terms of other physical constants of the material. 
Curves are drawn showing these quantities as functions of temperature, 
and some comparison made with experiment. The relations of the 
constants of the ceramic to those of the single crystal are then briefly 
discussed. 


§ 1. INTRODUCTION. 


In an earlier paper we discussed the properties of single domains of 
BaTiO, in a phenomenological way, that is we used some of the known 
physical quantities to predict the values of others. We were able to 
account for the existence of three transition temperatures, and to 
correlate the observed dielectric constants, spontaneous polarization and 
spontaneous strains. In this paper we extend the calculations to other 
physical constants : first we calculate the piezoelectric constants ; then 
we use these to calculate the clamped dielectric constants, which are 
very different from the free ones ; finally we calculate the elastic constants 
for constant field, which are very different from those for constant 
polarization. These calculations are on similar lines to those made by 
Mason (1948) for the ceramic. Finally, we discuss tentatively the 
relations between the constants of the single domain and the ceramic. 


§ 2. THEORY. 

The method we used in Part I. (Devonshire 1949) was to expand the 
free energy as a function of strain and polarization, taking it to be zero 
when the crystal was in equilibrium and had cubic symmetry; the 
tetrahedral, orthorhombic and rhombohedral states were treated as strained 
forms of the cubic state. The expression for the free energy given in I., 
9.21 was 
A=$eh, (eet yp +22) Cho(Y Zot Het Cay) + de4s(te+Y, +2) 

og tb) Ps ope, be tbe sex (bat bat ee) 
+$Ef,(PS+P)+P2)+ dot0(PGP2+ Pz + PePy) + o6"(Pe+ Py + Py). (2-1) 
* Communicated by the Author. 
SER. 7, VOL. 42, NO. 333.— OCT. I95I 4D 
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The c?’s are the elastic constants of the cubic crystal and x” its reciprocal 
susceptibility. If we regard the strains as of order P? equation (2.1) 
includes all terms of order P* or lower and some of order P® (the latter are 
omitted in Part I.). Equations for the stresses and the electric field 
can be derived from (2.1) by the usual thermodynamic formule : 


oA 
Xx,=— aE —ch a, — choy y—Cho%,—GirPs—G12P3—G1aP 2, 
x 
oA ; 
Y= aa ou. = —ChYo—Gaak yP 
2 (2.2) 
oA 


B= Se =Polx! + EyP2+ ElaP2 + aP24L'P) 
+2P (911% 2 +9 19Yyt9 12%) +P yGsstytP Gases 


If we solve these for x, etc., then we have 


C= —si,X,—8}o¥ y—8}2L,+ Q11P?+Q,,P?+ Qi2P%, 


eo et (ets eo) Wey, @ TS. +6" Vee ee) wt ee ey 6) heme 


ee cinco tit ATA Ss (2.3) 
B,=P2(x" + £41P2+ fjoP24+ €)P2+ 6'P3) 
+2P,(Qi1X p+ Qi2V y+Q12Z,) +P, Qua XytP.QuZe, 


e/) 0: 20) (eo gey Fe se) \e, 16,0e ey Comal el ie) 6.5.06 nes 


(8t;—8i2)(¢ii— Cia) = 1, (st, + 2sj2)(ct1-+ 2cts)=1, sicyu=1, 
Qi1— Qi2= — (81 —$12)(911—- 919) 
Qi t+ 2Q12= — (811+ 28}2)(G11+ 2919), Qa= — 859 a4» (2.4) 

S=Srr + 2911911 + 4912Qi2, 

S1a= S10 t+ 2(911Q12+-912Q11 +912Q12) + Gas Qas- 

These equations are now in a similar form to those given by Mason (1948) 

for the ceramic; but in that case the substance when unstrained and 

unpolarized is isotropic instead of cubic, so that 
Siu= 2(811—S9), ] 
Quu=2(Qi1—Qis). | 

Mason also uses as a variable D/47 (where D is the displacement) instead 


of P; this, however, makes little difference, since P is always large 
compared with E. 


(2.5) 
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Now as barium titanate is cooled it passes successively through cubic, 
tetragonal, orthorhombic and rhombohedral phases; and in the last 
three it has permanent polarization along an edge, face diagonal and body 
diagonal, respectively, of the original cube. We shall, therefore, consider 
in turn the physical constants for each of these phases. 

In the tetragonal phase let the permanent polarization P, be along 
the z-axis. Therefore in equations (2.3) let us replace P, by P,+-P, (so 
that P, is now the difference of the polarization from its saturation value) 
and neglect all terms in P?, P? and P?. The equations then become 


Uo =Q12P}—8},X,—Si2Y y—SigZ,+ 2Q1oP.P,, | 
Yy =Q12Pi—sppX,—St1Y y—8i2Z,+2Q,5P,P,, 
Bz =Q1iP3—sipX,—Sp2¥ y—S11Z,+ 2Q 1,P,P,, ‘ 
Yz =—84Y,+QuP,P,, 


> 


ee =—81,Z,,+Q P,P, te (2.6) 
Ly =—Sisgdry, 


E,= XuPs Sir QuP.Z,; 
E,= XuPy 2 QuP.Y., 


i, =XgsP,+2Q12P(X2+ Vy) +2Qi1P.Z., j 
where 
ae eae | (2.7) 
Xs =X" +3€1,P7+5C'P¥. 


Now the linear relations between stress, strain, polarization and field 
for any crystal can be written in the Voigt tensor notation : 
X= — chit +dm Pm (a) | 
Eo), Xen (0) 
Xp=— Cpt item rim, (¢) 
Pa emetat tints OL 8 eg) 
Ey =— SH Xi tOm Pins (6) 
E,,= Om 1X pele ee) 
Ly — SR Fn Emmy (9) 
P= —Amn Xn t+NimErrs (2) J 
where the summation convention has been used, that is each product is 
summed over all values of any suffix occurring twice in it. In the 
Voigt notation 


tg Ca Y ys V3— Fe Chea: (2.9) 
LY ry, Up 2g Xz) Le=XLy=Ya 
with a similar notation for the X’s, so that the suffixes ) and 7 run from 


1 to 6, and k and m from 1 to 3. The physical meanings of the various 
constants are implicit in the equations; thus the d’s are clearly the 


4D2 
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piezoelectric coefficients as normally defined. Between the constants 
there are numerous relations, which can easily be derived from the 
equations: if, for example, we eliminate P between (e) and (h) and 
compare the resulting equation with (g) we find that 


os nase (2.10) 
“ d pn=Om WMkm* 
Other relations that we shall need are : 
B= enChis 
Me Xtm to k 2m h 
= Ogi a CH Ba ar (2:10) 
Xkm kivmhvhi? ° e 
Cni8ig—=O np: 
XKIMIm= 9 tems 


where 56 is the tensor defined so that it is zero if the suffixes are different 
' and 1 if they are the same. 
By comparing (2.8 e) with (2.6), we see that in the tetragonal phase 


by ;=b12=b13=b4=b,,=0, by5=QaP,, 
b3=022=b93=bo5—=bog=0, bya=QuP,, : (2.12) 
b313=b32=2Q10P,, b33=Q11P,, bg =b3,=b34=0. 


We also have 


M=Ne=UxXi1» N3s= 1/ X39 (2.13) 
2=31=N23= 0, nS tea 
and 
C}1=Cp=C33, Clo= C51 Cg, Chy= CLs —= Co, 
(ci, -+ 2ci,)(st, + 2sf)=1, 
(2.14) 


(C11 —Ct2)(811 —S}o)=1, 
chuSiy= Ll, 
while all the other c’s vanish. Hence from (2.10) and (2.11) it follows that 


dy, =, .=d,,=d,,=d,,=0, ds5=QuPenii 
dy) =da.=dy3=dy,=dy4=0, da=Qu Poni (2.15) 
d3)=d32=2Q12P 733, d3g=2Q11P 033, d34=d3,;=d3,.=0, 
and 
S11 =822= 81, +4Q?,P ngs, 833= 811, +4Q}, P2nss, 
812=8jo+4Qi,P?ng3, 813=8)3= 819+ 4Q,1Q1 oP 23s, - (2.16) 


Bes eh apes iets 
$44=855= 844+ Qi,P i711, Seg=S44, 
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and 
X11 = X22= X11 +64, Q7,P2, | (2.17) 
X3s=Xast-4 {cr Qh +4c12Q11Q19-+ (CF, +e%) 12} Ps. 

In the orthorhombic phase let the permanent polarization have: 
components P, along the y and z-axes of the crystal. In equations (2.3) 
let us replace P, and P, by P,+P, and P,+P,, and neglect all terms in 
A ea wae Bal Bes ‘The equations then become 
aoe Cae $— S11 g— Sto ¥ y—StoZ,+ 2Q1 P,P, + 2QiaPP., 

=(Q11+Q12)P3—sppX,— str Y y—SioZ,+2Q1,P, P+2Q;2P.P,, 
Ry | Se —Sj2X,—8i2Y y mpat +2Q12P,P,+2Q,,P,P., 
Ye =QuPi—siiY,+QuP, sree ts st 2 
ig > —8iZ,+QuP,P,, (2.18) 
Ly = —si,X,+Q,,P.P,, | 
B=xXnP e+ QuPZe+QuP,Xy, 
By=x33P yt xX23P-+ 2Q 12PX2t+ 2Q11P,Vyt+2Q12P,Z,1+ QP Y, 
E, =X25Py+xssP2+2Qi2P,X_+2Q12P,¥ yt 2Q11P,Z.+ QuP Ve 
where 
Xiu= xX" +2f oP 5, 
Xas= x" + (3E1 + £49) P51 50'Py, - + . (2.19) 
Xo3= 2612P 5 
By comparing these equations with (2.8) e or f we see that 


by4==51,=b13=b14=0, by5=016=QauP, 
bo3=b3;=bo3=b32—=2Q1 oP, bog=b33=QiiP,, a 40 (2.20) 
Doy=b34= QasP Dys=0og=b35=036=9. 


We also-have 
mi=1/xi1 2=N113=9, Nee=Nss> | _ (2.21) 
(33-+753)(x33-+ X23) =L, (33—"23)(Xs3— X28) = 1, 
and the c?’s are given by (2.14), so that by substituting in (2.10) and 
(2.11) we have 
y= .=yg=dy4=9, dy5=d16=QaaP 11 
d.g3=d31=2Q4 oP. (433-723), 
do3=dgo—=2P,(Qi 233+ Q11793) . (2.22) 
dog=dzg—=2P,(Q1 1733+ QieN2a); 
do4=A34=QuaP (33-193); 
dy5=dag=d3,=d34=0, 
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and 
2 p2./ y Hien yey 
s¥ =8}, + 8Qi2P3(n33-+N23), S59—= 833 = 811 | 


+ {8Q11Q1 9793+ 4( 71+ Qio)n33 }P%, 
si = st = sto +4Q1o(Q11+Q12)PF(733+723), 
853—=Sjo+ {8Q11Q19733 +4 71+ Qis )nos}P | 


514=4Q2Q44P3(733-+7123)5 S34=S51—=2(Q 11 +Qi2)QuaP5 (n33-+ 723) Cae 
815=835=855—=S816= 836 836= 0 
83, = 8}, +2Q,2P?(33t+193), 855-=860= 844+ QP emia 
$= QPP s=8=0, | 
and 
Xu=Xit2¢f,Q2,P%, 
X22= X33= X33 
+[4{QiicF, +4Q11Q1 oC} + 2Q7(ct, +cf2) }+-cF,Q7,]P2, (2.24) 


Xer=Xi2=0, 
X23= X23 : 
+[4{Qhicto t+ 2Q11Q10(¢t +¢i3) + Qio(Ch + 8cts) } + QF,ch JP 

In all these equations for the orthorhombic phase it must bs remembered 
that the resultant saturation polarization is P,,/2. 

In the rhombohedral phase the permanent polarization will have the 
same component P, along each axis of the crystal. Hence in equations 
(2.3) we replace pe P, and P, by P,+P,, P,+P, and P,+P,, and 
neglect terms in P2, P? and P2. The equations then become 


%_ =(Qyi+ 2Q,2)P?—s},X, —sto¥ y—S}oZ,+ 2Q11P.P, 
+2Q,2P,P,+2Q,.P,P,, 


ete ee ee ee oe en AT eer ty te. Ty U4 


(2.25) 
E,=xpyPetxiePytxieP+2Q11:P,.X_+ 2Qi2P,Vy 
at 2Q; aP,Z,+ QuP.Z,+ QuP,X,, 
where 
Xu=X" + (8&1 +2€),)P2+50'P4, 296 
Xia= 26 ;oP 2. pee 


By comparing these equations with (2.8) e or f we see that 
by, =b22=b335=2Qi1P,, 
by 2=b13=b2;=bo3=bs;=b3o=2Q) oP, 
bb eben, Pad 
boy =b34=b) 5=b35=b1 «= bog=Qy,P, 
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We also have 


ies Nes, "23—=Ns1—Nies 
(2.28) 


(M1 +2019)(x11 + 2x12) = iF (X11—X12)(M11— M12) = 1, 


and the c*’s are given by (2.14), so that by substituting in (2.10) and 
(2.11) we have 


diy =dy.=d33= (2011711 +4Q, amie) Pes 
d9=d13=d9)=d.3=d3,=d39= {2Q1 9711 +2(Q11+Q12)ni2}P,, 


dy 4=d5=Agg=2Q gan ioP ss a) 
d45=4g=d y= dy g—=dgy=dg5= Qa, (N41 +-712)P } 
and 


. $11 =83=835= 8], + 4{ (Qi + 2Qjo)n11 + (4Q11Q19 + 2Q7) N19 }P, 
Gee) oP 
Oeee Sta O12 O10 
+4{(2Q5,Q) 0+ 935): +(Q71+2Q11912+3Q?,)n12}P2, 
Sh4=855= S36 + 4{Q1 2Qaa ti + (Qi1 + Qi2) Qa 12 }P?, 


S15= S16 S24= 826= 535= 836 oe 
=2{(Q11 + Q12) Qua + (Q11 +3Q12)Qua712}P 5, 
844= 855 =S86g—=S4r + 2Q (N11 +712)Ps, 
856=S§1= S15 = QD +3 19)P ss 
and 
Xt1=X22=X33= Xi 
+ [4 {cf Q71+ 4c72Q11Q19 + 2(ch, +Cj2) Qiao} + 2cLQUIPS, '(2.31) 


RSS 
soles {c7,Q7,+ 2(cf, +¢7)Q1,Qi2+ (cy + 3cf,) to}+ef,Qia]P2. 


Tn all these equations for the rhombohedral phase it must be remembered 
that the resultant saturation polarization is P,/3. 


§ 3. NumERIcAL RESULTS. 


In the previous sections we have found formule for the piezoelectric 
coefficients (d), the elastic compliance coefficients for constant field 
(s"), and the reciprocal susceptibilities for constant strain (y”). To 
evaluate these formule we must know the permanent polarization 
P,, the compliances and elastic constants for constant polarization 
(s? and c?), the susceptibilities for constant stress (y’) and the electro- 
strictive coefficients for the cubic form (Q). In Part I. we showed how 
to calculate P, and y’ (and hence 7’) as functions of temperature by 
using a limited number of experimental results, and we shall use these 
calculated values here. The electrostrictive constants Q can be determined 
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by comparing the measured spontaneous strains (Kay and Vousden 
1949) with the measured spontaneous polarization (Merz 1949 a); in 
this way we find that 

Q.4= 29x 10-2, 

Qi — 1 x Og) Pe Piety. ea 8 (3.1) 
and Oi a Oe 


in c.g.s. units. These values can be checked by using Merz measurements 
of the shift of transition temperature with hydrostatic pressure ; for it 
can be shown that if T, is the transition temperature 


a, 2x" 
op OT 
Now, according to Merz, the upper transition temperature is lowered 


by about 1° C. per 170 atmospheres, and measurements on the dielectric 
constant in the cubic phase show that dy"/0dT is about 10-4, so that 


Qi1+2Q,.=0°6 x 10-”, 


which agrees reasonably well with (3.1). Elastic constants for the single 
crystal have been measured by Bond, Mason and McSkimin, who found 
that 


zeny SCTA , a 


cl, =2-07 x 1022, 
<: SeeRU O 


ch, = 1-27 x10”. 


There are no measurements of ci,, so we have assumed it equal to 
ck, ; we then find that 
s?, =0-90 x 10-22, 
s?, = 0-34 x 10-22, pi i eee re 
si, =0-79< 10-22. 


Both the Q’s and the s’’s must be functions of temperature, but the 
temperature dependence has not been measured, so we have treated 
them as constants, which is probably not far wrong. 

It is interesting to compute the quantities g and é”, since they are 
directly related to the forces between the ions. We gave values for these 
in Part I. (Devonshire 1949), but these now need revision as more 
experimental data are available. From (2.4), (3.1) and (3.4) we find that 


Ju= —3-2, 
J12=9, 
Ju= —3-4, 


and using the values for €;, and €;, given in Part I. equations (7.14) and 
(7-15) we have 

Ey =141x 10-2, 

Eo 7°4X 10-22, 
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We can now calculate the components of d, s® and n” as functions of 
temperature, and the results are shown graphically in figs. 1, 2 and 3. 

The piezoelectric constants are seen to be very large and to vary 
rapidly with temperature. Piezoelectric constants are usually of the 


Fig. 1. 


rhombohedral orthorhombic 


(45) 


5x40 


cs 


-10 


-15 ee os 
700° - —160° —|20° —80° —40° 0 40° 80° 1210° 

Piezoelectric coefficients d of BaTiO, calculated as functions of temperature. 
The suffixes in all cases refer to the principal axes of the cubic or 
tetragonal phase; these are not, of course, principal axes in the 
orthorhombic and rhombohedral phases. The directions of polarization 
are (001), (011) and (111). Values in ¢.g.s. units. 


order of 10-7, so these are about a hundred. times larger; these large 
values are also found in other ferroelectrics such as Rochelle salt or 


KH,PO,. 
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The compliance coefficients for constant field differ greatly from those 
for constant polarization ; they may, indeed, be several times larger. 
Their extreme anisotropy is remarkable if we remember that the departure 
of the axial lengths from cubic symmetry is only about 1 per cent. The 
peaks near the transition temperature in some components can easily be 


Fig. 2. 
6xI0- 
rhombohedral orthorhombic tetragonal cubic 
} 
Sy 
4 
3 
ie ea? 


Ir Sé& = — 
Si Ski 
= pn Seg Soe ae 257 
Sk 
)/_—— eS + + + 
E Si 
St Se ee Siz 
E gah nd ME aah GR gore Reem S 
Ps Sia Si2= S\3 Sis SE 
oA 
2 
E 
S33 
—3 
— 160° —120° —80° —40° 0 40° 80° {20° 160° 


Elastic compliances for constant field of BaTiO, calculated as functions of 
temperature, Axes as in fig. 1. Values in c.g.s. units. 


explained : for example, the rise in s!, and sas we approach 0° from above 
means that the substance can be easily sheared in the yz or xz planes, which 
we should expect, since at the transition it shears spontaneously in one 
of these planes. 
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There are not many measurements of the d’s and they are confined 
to the tetragonal phase. The coefficient dj, has been determined by 
Merz and Caspari by measuring the extension produced by a known 
electric field; they found that d,, is about —10-® at room temperatures 
with a very sharp rise to —7 x 10-* at the upper transition temperature. 
These values are less than the calculated ones ; this is probably partly, 
but not entirely, due to the fact that they were measured with a strong 
biasing field in order to keep the crystal as a single domain. The 
piezoelectric coefficients can also be measured by a dynamic method ; 
by measuring the resonance and anitiresonance frequencies of a crystal 
oscillator both the piezoelectric constants and the elastic compliances 
can be determined ; the theory of this method has been given by Mason 


Fig. 3. 


1000 


00} Phombohedral orthorhombic tetragonal 


= 
—— 
-_— 


=— 
= 
<== 

_— 


1-0 
—160° = i20g — 80° —40° 0 40° 80° 120° 160° 


Dielectric susceptibilities for constant stress (7’) and for constant strain. (7") 
of BaTiO, calculated as functions of temperature. The suffixes refer to 
the same axes as in fig. 1, so that 733-73, etc., are the values for the 
principal axes of the orthorhombic or rhombohedral phase. The 
susceptibilities are related to the dielectric constants (1) by the equation 
e=1+4z7n. The vertical scale is logarithmic. 


(1950). Dynamic measurements of d;, have been made by Caspari and 
Merz 1950), who found a sharp increase near the upper. transition 
temperature as in their static measurements, but they did not measure 
absolute values. Dynamic measurements at room temperature have been 
made by Bond, Mason and McSkimin (1951) who found that 

dg, 3-1 x L0r®, 

dlag= 10-4 x 10-8, 
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These values are higher than the theoretical ones ; the crystals were 
not, however, strictly single domains. Measurements of s® were also 
made, and these agree fairly well with the theoretical values. ; 

Fig. 3 shows that the “ clamped ” susceptibilities (7”) are smaller than 
the free susceptibilities y’ by a factor that varies from 2 to 10; this is 
a very large effect. The ‘“‘ clamped ” susceptibilities can be determined 
by making measurements at frequencies greater than the resonance 
frequencies. Measurements made by Bond, Mason and McSkimin (1951) 
show that at room temperatures 733 is halved when measurements are 
made at frequencies higher than the resonance frequency. Similar 
measurements have been made on ceramics, but we shall consider these 
in the next section. 


§ 4. PHysicaL CONSTANTS OF CERAMICS. 


As most of the experimental work on BaTiO,, and all the earlier work, 
has been done on the ceramic, it is desirable to consider its properties. 
It must in principle be possible to derive the constants of the ceramics from 
the constants of the single domain—one would expect the ceramic 
constants to be some kind of directional average of the single domain 
constants. The average may, however, be made in many ways; for 
example we can average either the susceptibilities (7) or the reciprocal 
susceptibilities (vy); again we can average either the compliances (s) or 
the elastic constants (c), and since s® is the same as s? for the unpolarized 
ceramic we can average over s", s?, cP or cP. For a nearly isotropic 
substance it matters little how we take the average, but for a substance 
so very anisotropic as BaTiO, it matters much. Now the method of 
averaging depends on the boundary conditions ; if the boundary condition 
were that stress is continuous then it would be correct to average over the 
s’s and not the c’s. Boundary conditions are not, however, so simple ; 
if there is no slip then three components of stress and three of strain are 
continuous ; if the boundary is perpendicular to the z-axis then the 
continuous stress-components are T,,, T,,, and T,,, and the continuous 
strain-components are 8,.,, S,,, and 8,,. Again the electrical conditions 
at a boundary are that the normal components of polarization* and the 
tangential components of field are continuous. Thus we cannot take a 
simple average for any of the physical constants, and the correct way of 
taking the average is not apparently known. 

The directional average of a symmetrical tensor of the second order is 
given by 


Ay=$(Ay, +Ago+ Ags). <te! Cote Leer as tae 


For a fourth order tensor having the symmetry of the elastic constants 
there are two independent averaged quantities, namely : 
SE Eee 


# More precisely of displacement, but the distinction is of little importance 
for BaTiOg. 
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Byii1=3(Bi111 + Byo92-+Bssa5) + 3 (Bagge + Bgg4+ By 129) 


+7,(Bogo3+ Bsis1-+Byar9), Tan (4.2) 
and 


Byj22=75(Bir11 + Bose2+ Bsgg3)-+4(By 122+ Boogst+ Basi) 
— 2 (Bisi2+Bose5+ Bai31); 


Byore=2(By111—By, 22): 
When translating these into Voigt notation we must remember that 


while 


Ssa= 481019, 
- but that 
C4g—=Cy219- 

As an example of the effect of different methods of averaging we shall 
calculate s,, and s,, in four different ways. We pointed out in the last 
section that we could average over s® or s?, or alternatively average over 
c® or cP and then calculate s by using the usual relations between the s’s 
and thec’s. We give the results below, where we have used the theoretical 
values for s® and s? at 20°C. The experimental values are those of 
Mason (1948), obtained by resonance methods. 


Average Average Average Average Experi- 

over s? over s= over c? over c# mental 
81, 10}* 0-56 1-50 0-41 0-83 0-88 
8,5 X 10% —0-17 —0-61 —0-13 —0-29 —0-236 


The problem of averaging the electro-restrictive constants is even more 
difficult, since the strain corresponding to a given polarization is not 
necessarily always the same. The polarization of the ceramic may be 
produced by altering the direction of some of the domains or by changing 
the magnitude of their polarization, and the corresponding strain will be 
different in the two cases. Mason (1950 b) found that for the ceramic 


O= 3°6 x 10-22, 
Qpe 1-35 K 10-22. 

The values obtained by averaging the constants for the cubic crystal are 
Qoa== 0-6 10577, 

Fig. 4 shows the susceptibilities for both the free (n’) and the clamped 
(7”) material averaged in two different ways—directly and over the 
reciprocal susceptibilities. The differences are sometimes large ; and for 
7’ the discontinuities at the lower transition temperatures are reversed. 
Measurements of the free susceptibilities for the ceramic have been given 
by many authors. In particular an extensive series of measurements 
have been made by Von Hippel, Breckenridge, Chesley and Tisza (1946) 
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and measurements for temperatures down to a few degrees absolute by 
Blunt and Love (1949). The measured values follow a curve which, for 
the most part, lies in between the two mean curves for 1. 

The clamped constant 7” would presumably be obtained if measure- 
ments were made at such high frequencies that the strain could not 
follow the polarization. Powles (1948) found that at a frequency of 
9:45 x 10° cycles the dielectric constant had dropped considerably (by a 
factor of about 5). A further drop occurred at a frequency of 
24-5 10% cycles. These frequencies are extremely high for piezoelectric 
resonance frequencies (e. g. those found by Mason for his single crystals and 
polarized ceramics were about 10° cycles). However, in an unpolarized 


Fig. 4. 


1000 Ws 


tetragonal 


160° 120°. ~~ —g0® 40° 0 40° 80° 120° 160° 


Dielectric susceptibilities of BaTiO, for constant stress (y’) and for constant 
strain (y”) averaged over all directions. The full lines represent direct 
averages, and the dotted lines averages made over the reciprocal 
susceptibility. ; 


ceramic, which consists of polarized domains pointing in various directions, 
there will be no vibration of the material as a whole in an alternating field, 
and the effective resonance frequency will be approximately that of a 
single domain ; since resonance frequencies are inversely proportional to 
linear dimensions it is possible to have high values if the domains are small 
enough ; Powles results suggest that the domain size is of the order of 
10° mm. ; this agrees with measurements of the “‘ Barkhausen ” effect 
by Newton, Ahearn and McKay (1949). 
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A different explanation of Powles’ results has been put forward by 
Mason and Matthias (1948). In their model of barium titanate the 
titanium ion alternates between a number of positions of equilibrium ; 
there is naturally a relaxation frequency associated with the movement of 
the titanium ion, and if we make the right assumption about the size of the 
potential hill separating the equilibrium positions it will correspond to the 
frequency at which the dielectric constant falls. This explanation is 
supported by the fact that Powles found a drop in dielectric constant even 
above the Curie point ; here the crystal is not piezoelectric, so that there 
should be no difference between the free and the clamped dielectric 
constants. The fall in dielectric constant is, however, much less above 
the Curie point, and in a ceramic there will probably be internal strains, so 
that there may exist tetragonal regions even above the Curie point, as has . 
been suggested by Caspari and Merz (1950). Moreover the difference 
between the “free ”’ and ‘“‘ clamped ”’ constants certainly exists and is 
about the right amount to account for Powles’ results ; its existence does 
not depend on any particular model, and it must be observed if we go to 
sufficiently high frequencies. This does not mean of course that there 
may not be other relaxation effects associated with a particular model such 
as Mason’s. 
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ABSTRACT. 


In §1 it is shown that Gorter’s assumption that the forces acting on 
the normal and superfluid are equal and opposite is true only if there is 
no entropy of mixing: §2 shows that his derivation of the fountain 
effect pressure assumes also that entropy is proportional to concentration. 
In §3 we point out an additional assumption implicitly made by Gorter, 
Tisza, and by Nakajima and co-workers—that effective mass is pro- 
portional to concentration. §4 is devoted to a proof that whatever the 
energy spectrum the entire entropy and internal energy and all the 
excited atoms (or excitations) follow the motion of the tube in which they 
are contained; this is not in agreement with Temperley’s suggestion 
that other states besides the very lowest one contribute to the superfluid 
density. In §5 we calculate the extra terms arising from the creation 
of new phonons during the motion of the normal fluid. The status of 
the relation p, «<S in the two-fluid theory is discussed in §6. The next 
section is concerned with the theoretical condition determining the 
lambda-point, and §8 with the expected differences in behaviour of 
Het and He’. 


§1. THe Forces AcTING ON THE NORMAL AND SUPERFLUID. 


GORTER (1949) has recently attempted to generalize the equations of the 
two-fluid model by discarding the usual assumption that there is no 
entropy of mixing. In his derivation he assumes that when there is only 
a temperature gradient in Helium II, the normal and superfluid com- 
ponents diffuse into each other, and that the effect may be described by 
a diffusion force equal in magnitude but opposite in sign for the two 
fluids: this means that K,(T)=—K,,(T) in his notation. From this 
assumption the values of K,, and K, are found by an argument analogous 
to the standard derivation of Kelvin’s equations for thermo-electricity. 

Let us now examine the possible bases of the assumption that K,=—K,,. 
First, we note that it does not follow from Newton’s third law—that 
action and reaction are equal and opposite—since each force vanishes 
under equilibrium conditions, 7. e. when there is no temperature gradient. 
Secondly, we observe that it implies that the rate of change of ‘momenta 
of the two fluids are equal and opposite, so that 0(p,v,)/0t= — O(p,,v,,)/Ot, 
whence p,v,+p,¥,=0 and p,x,+ P,tX,=0. Thus Gorter’s assumption 


ee een 


* Communicated by the Author. 
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implies that the centre of gravity of the two fluids remains fixed if a 
temperature gradient is imposed. Although irrelevant to our argument, 
we may note that there seems to be no reason why this should be generally 
true ; even in ordinary gaseous systems thermal diffusion may lead to 
its disease. 

In a previous paper (Dingle 1949 b) we showed how the forces K, and 
K,, may be calculated separately. Generalizing the model discussed in 
that paper, we now determine the conditions under which K,=—K,,. 
In the notation of that paper, the force acting on the superfluid per unit 
volume of Helium II is, for any two-fluid model (loc. cit. § 4 (i)), 


a sw a 88 a 8W a 8 a 8W 
oe (5 V. Tme,7 (5 Tm)e, .: (5 an 7 (5 ve, v(G =) 


(1.1) 
where 5W is the change in potential energy. Since the total mass of the 
fluid is conserved, 


Sp/p+V .x=0, ree Sass hae (12) 
where x is the coordinate representing the centre of gravity of the two 
fluids. As.we have seen, this is stationary on Gorter’s assumptions 
provided there is no pressure gradient, so that 6p=0, and the second 
term in (1.1) vanishes. Repeating the argument for the force on the 
normal fluid (loc. cit. §4 (ii)), we see that K,=—K,, only if 


05S 068 
(soa)e--(sveje 0 


Let us suppose that when the normal fluid moves a distance x,, a certain 
quantity of entropy p,S,, moves with it, and similarly for the superfluid. 
No assumption has been made as yet that there is no entropy of mixing, 
or that these entropies borne along by the fluids bear any relation to the 
entropy of the fluid considered as a whole, pS say. If no entropy is 
destroyed by irreversible processes, 


9 = 205) + ,8,V. Bee Orrin Sas. (14) 


Any entropy of mixing is here expressed by the fact that pS is not 
necessarily the simple sum of p,S, and p,S,, but may be a more 
complicated function. Now (1.4) may be written 


0= s{z + pyV.tn+psV - i} 
0s ; ; 
- {055 — (S—S,,)p, V - %—(S—S,)p,V . i} 5. ar ae Se ta) 


The terms in the first bracket together vanish, since the total mass 
p=p,»tp; is conserved*. Hence 


pdS=p,(S—S,) V - tm+p(S—S,)V.%, - - - + (1.6) 
* See the final paragraph of §4 of the present paper. 
SER. 7, VOL. 42, NO. 333.—OCT. I95I 4E 
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and (1.3) is satisfied only if 

pPS=prSntpsde oe ae 
The condition for the validity of Gorter’s assumption is therefore simply 
that entropies are additive. Gorter’s results do not therefore appear to 
constitute any valid generalization of the theory, since they can only be 
correct when there is no entropy of mixing, the assumption made by 
previous workers. 


§ 2. THe Fountain EFFECT. 


Assuming that K,=—K,,, Gorter finds that the force per unit volume 
is given by 


een ie 6) VI, Si) eee 


where X is the relative concentration of normal fluid present. This 
force leads to a fountain effect of magnitude 


os 
Vp=pX aX VIL. ie ee ee ee 
This reduces to the customary expression (London 1939) 
Vp= ps VEO Sa Dey, 7 ee 


only if X(0S/0X)=S, 7. e. if SaX. We shall show in §6 that for this to 
be so the normal fluid must retain its identity when the entropy of the 
fluid as a whole changes, which is true only for an ideal gas. Gorter has 
indeed apparently assumed this in his derivation of the magnitude of 
the forces, for it corresponds in his thermo-electrical analogy to the fact 
that the same charge travels right round the circuit. We shall return 
to this point in §6, merely noting here that no such assumption is 
necessary to obtain (2.3), since the concentration of normal fluid does 
not enter London’s argument. 


§3. CONCENTRATION AND EFFECTIVE Mass. 


In writing down the equations of motion resulting from the forces K,, 
and K,, Gorter makes the further implicit assumption that the effective 
masses of normal and superfluid are pX and p(l1—X) respectively ; 
v.e. that their effective masses are proportional to their respective 
concentrations. In other words, it is assumed that the effective mass 
—defined as the quantity by which the velocity must be multiplied to 
obtain the momentum, and hence capable of measurement (or calculation) 
only by an essentially dynamic method such as sound propagation or the 
Andronikashvili experiment—is equal to the actual material mass of the 
component present, a quantity which could only be measured when the 
system is essentially static. These two masses are identical for a perfect 
gas, but not otherwise. This point is well illustrated by the model treated 
by Temperley (1951 a) in which the energy spectrum is assumed to be 
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of the form* E=Ap”. The concentration of excited atoms is then found 
to be proportional to T?/", and the effective mass to be proportional to 
T°"? For a perfect gas n=2 and the fraction by effective mass is equal 
to the concentration. 

The work of Landau (1941), Ward and Wilks (1951), and the author 
(1948, 1949), is free from the implicit assumption that the effective mass 
is identical with the actual mass. On the other hand, Gorter (1949) has 
implicitly made such an assumption by ascribing the same symbol X to 
each, whilst Tisza (1947) and Nakajima and co-workers (1950 b) have 
done so by first attempting to prove that entropy is proportional to the 
concentration of normal fluidt and then replacing this concentration by 
p,/p (the fractional effective mass) without giving any reason. 


§4. Tue Fiux oF ENTROPY AND ExciraTions. 


Assuming an energy spectrum of the form E=Ap”, Temperley (1951 a) 
compares the value of p,/p calculated as the fraction of the number of 
excited atoms with the value given by Landau’s method. Finding that 
when n+2 (i. e. for systems other than a perfect gas) the methods give 
different results, with the experimental evidence in favour of Landau’s 
treatment, Temperley concludes that this is ‘‘ tantamount to saying that 
if liquid Helium II can be represented as a condensed Bose—Kinstein 
assembly, other states of the material particles besides the very lowest 
one contribute to p,’’. We shall now show that this suggestion is unsound. 
Let F be the Helmholtz free energy, S the entropy, U the internal energy, 
and N,,,. the number of excited atoms (or excitations), all measured per 
unit volume. Now 


4nkT 


F=+ ” In (LEAexp{—E(p)/KT})p? dp . Saye 41} 
he) 


the upper signs being taken for the case of Bose—Hinstein statistics, the 
lower for Fermi-Dirac. §S, U and N,,,, may then be determined at once 
from the relations S=— 0F/éT, U=—T?0(F/T)/oT and kTN,,,= —AdF/aaf. 
If we now suppose the liquid to be contained within an infinitely long tube, 
and move the tube with velocity v along its own length, we may calculate 
the amount per unit volume of free energy, entropy, internal energy, and 
the number of excitations carried along with the tube. Let the flux of 
these quantities be F’v, S’v, Uv, and Nj,.v, respectively. To allow for 
the motion of the tube, E must be replaced by E—pv cos 0, @ being the 


* The usual statistical formulae involving the energy levels of individual 
atoms can then be used only if the interatomic forces are adequately described 
by a self-consistent field. 

+ Which is itself unsound, as we shall see in $6. 

+ We are taking into account only the free energy contributed by the 
excitations, the lowest state p=0 being automatically excluded from (4.1) by 
the presence of the statistical weight factor p’. 


4E2 
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angle between p and v. Since an excitation moves in the direction of 
its momentum p with a (group) velocity given by Hamilton’s equation 
as 0H /dp=dH(p)/dp, the flux of free energy is 


Font | [ |" In {1A exp (pv cos 0—E)/kT)} 
p=0/ 0=0"¥ g=0 


he 4 

Bs 2 si 6.d 42 

[ate 8 DEAT UD Gee, +t ge en eee 
Ak LU sh Soe ee ee ee ae ieee, _nyet) 2E 

=F = d = In (1 Ae P**) — pw d 
=. rea) cos? @ . sin 0 a8 | b| FTO (1 Ae Ui p? dp 
(4.3) 

to the first order in v. On integration by parts we see that 

Ff’ = spats = : 2m | In (1-FA exp {—E/kT}) 3p? dp=F, (4.4) 


the integrated term vanishing provided that E(p)>oco as p>oo. By 
differentiation with respect to T and A, we deduce that S’=S, U’=U, 
and N‘,.=N,x- also. We have therefore proved that for any energy 
spectrum the entire entropy and internal energy and all the excited 
atoms (or excitations) follow the motion of the tube in which they are 
contained ; 7. e. none of them contribute to the superfluid. 

The lack of proportionality between the effective mass and concentration 
is therefore not due, as Temperley suggests, to the fact that other states 
besides the very lowest one contribute to p,. Instead, it appears to 
be a consequence of the fact that if there are interactions each energy 
level has associated with it its own effective mass. For if we take 
E=Ap"=}p?/(p?-"/2A), it is plausible to suppose that the effective mass 
per state m(p) is proportional to p?-” ; averaging over all momenta then 
gives 

- | pnp? dp 
Pn A-T exp (Ap"/kT)+1 


which is to be compared with 


el btietemnemem AL oo Yh 


PEED 3 
Nee Iya exp (Ap"/kT) +1 COMMS A eS (4.6) 
This result for p,, is just that found rigorously by Landau’s method. 
In the notation used above, the effective mass, defined as the quantity 
by which the velocity must be multiplied to find the momentum, is, for 
small velocities, 


=” phd 2 | \ 
po he 8)? °P OE texp (Eel) zly °° 7 MS) 


per unit volume. The effective mass cannot therefore be obtained from 


the free energy simply by a process of differentiation, as can the other 
variables. 
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Landau’s method of determining the value of p,, involves a steady flow 
of normal fluid moving forward with the tube at a velocity v, the super- 
fluid remaining stationary. If the centre of gravity of the two fluids is: 
required to remain at rest, as in second sound, we must imagine that a 
piston is introduced into the tube and moved with a velocity u=—p,,v/p 
relative to the tube. Such a piston will move the liquid as a whole, 
giving it a momentum pw, which must be considered as the sum of the 
extra momentum p,w imparted to the normal fluid, and a momentum 
p;v transferred to the superfluid. Thus we have p=p,+p,. The sum 
of the effective masses of normal and superfluid is therefore equal to the 
actual material mass of the fluid considered as a whole. 


§ 5. THe CREATION OF PHONONS. 

In the particular case of phonon excitations, E=pc, it is possible to 
determine EF’, 8’, U’, Ni... and p,, without approximation. The results are 
Payee Uae No DY 

ee el ee rate, glk Ot | 
Fos =U-Ne=('-¥) oy 
4U Ta 
and m= 5a(1— 7) Ree eric... (0.2) 


3c? 
where v, is the velocity of the normal fluid. The terms in v2 are due to 
the creation of new phonons during the motion. Such terms need not 
be considered in second sound propagation, as they are of higher degree 
in v, than those leading to shock waves: (cf. Temperley 1951 b, Osborne 
1951). On the other hand (5.2) gives a critical velocity above which there 
can be no superfluidity ; this critical velocity, 


4U \ 113) 1/2 
tea={1 — (55) ; 5 0: WG0) Gly Ge ace oe (5.3) 


might be a determining factor close to the lambda-point. 


§6. THE STATUS OF THE RELATION p,,«S IN THE Two-FLuIp THEORY. 


Apparently not accepting Landau’s method of calculating p,, Tisza 
(1947) attempted to prove that p, oS. He imagines a cell fitted with a 
piston impermeable to normal fluid. On moving this piston inwards, a 
certain quantity of superfluid leaves the cell through the pores of the 
piston, but no normal fluid can leave. The total entropy in the cell 
remains unchanged since the process is supposed adiabatic. According to 
Tisza the concentration of entropy will therefore change in the same 
proportion as the concentration of normal fluid, so that p, «S. 

This argument appears to be in error on three counts sai 

(1) It is implicitly assumed in the derivation that since no normal 
fluid leaves the cell, the total amount in the cell does not change. This 
is not generally true, since the density of normal fluid is fixed by the 
temperature of the contents of the cell; this is decided by the entropy 
concentration, which in its turn is determined by the essential condition 
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that the total entropy within the cell remains constant. (If the normal 
density is not the same as that given by Tisza’s relation, a relaxation 
effect will probably take place, its value assuming at once the value given 
by the relation, and then tending towards its equilibrium value determined 
by the entropy concentration. ) 

(2) The process is assumed to be an adiabatic one, and any relation 
derived merely provides us with a method of calculating the change in 
temperature. within the cell. Without further assumptions (such as 
that the normal density is independent of pressure) we cannot use the 
relations obtained to determine equilibrium values for given temperatures 
(cf. Dingle 1948, §12). 

(3) It is implicitly assumed that the effective mass-density p, 1s 
proportional to the concentration of normal fluid. We have shown in 
§ 3 that this is true only for an ideal gas. 

Nakajima and Shimizu (1950 b) also make these assumptions. They 
attempt to justify the supposition that normal fluid mass is conserved by 
remarking that ‘‘ the two-fluid theory is formulated self-consistently with 
respect to this assumption, for we can derive the following equation from 


this theory : 
1 ¢/d8 g.do 
| 4 ; i - S; | (Fe Ag ¥ J 45 } 


where J’ is the rate of transformation from the superfluid to the normal 
fluid per unit volume...”. On referring to the proof of this expression 
(1950 a), one finds that the Tisza relation has been used in its derivation, 
so that the argument given by Nakajima and Shimizu is merely a circular 
one. The particular two-fluid theory with [=0 and p,«<S presented by 
these authors (1950 a) is cerfainly self-consistent ; but so also is a general 
theory with [40 and p, not proportional to S. 

Nakajima and Shimizu (1950 b) point out that it is curious that the 
normal fluid should be conserved in an adiabatic change, but show that 
this is indeed the case for a condensed perfect Bose-Einstein gas. We 
remark only that if a general energy spectrum is introduced to allow for 
interactions, the normal fluid is no longer conserved. 


§ 7. Tur Lamppa-Pornt. 


All authors have tacitly assumed that on the two-fluid model the 
lambda-point is given by 


Pall pSe. 1.16 ply) =e0.5 a en een en 


This is in agreement with the experimental fact that the velocity of second 
sound w,=(TS?p,/Cp,,)!/2 tends to zero at the lambda-point. 

At first sight it appears that (7.1) has no theoretical justification. As 
we have seen, p, is an effective mags which can only be relevant when there 
is relative motion of the normal and superfluids. How then can it be 
involved in any relation determining the lambda-point in the specific 
heat, a phenomenon observed when the liquid is stationary ? The 
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lambda-point must be due to the breakdown of the two-fluid model, and 
should therefore occur at the temperature for which N,.,==N rather than 
that for which p,—p. In the general theory these two conditions lead to 
quite different values for the lambda-temperature. 

The experimental fact that the velocity of second sound tends to zero 
at just the same temperature as that of the thermodynamical lambda- 
point provides strong circumstantial evidence in favour of the author’s 
theory that the high entropy and specific heat of Helium II are due to 
supersonic waves of second sound (Dingle 1949 a, 1949 c). Such a theory 
leads to a critical point when the velocity of second sound tends to zero, 
and provides an immediate explanation for the appearance of the effective 
mass in the theoretical expression for T,. 


§ 8. ExprcTED DIFFERENCES IN BEHAVIOUR OF Het anp He’. 


Under some conditions a system of fermions (such as He*) may be 
described in terms of an equivalent assembly of bosons (Tomonaga 
1950). It seems plausible to suppose that this will be the case if the inter- 
action energy between the particles is much greater than the thermal 
energy, for then the actual statistics obeyed by the particles should be of 
little importance. If this is true, the behaviour of He® at sufficiently 
low temperatures should be similar to that of He*. For instance, second 
sound propagation should be possible, with its velocity w, tending to 
u,//3 at absolute zero, just as for Het. 

At higher temperatures the difference in statistics should become more 
marked as the thermal agitation of the atoms overcomes their mutual 
interactions. He* will continue to obey a two-fluid model because the 
normal and superfluid densities are comparable in a Bose—Kinstein 
condensation, but He® will cease to obey such a model because there is 
no condensation on Fermi-Dirac statistics and the superfluid density 
therefore becomes negligible much more rapidly as the temperature is 
raised. He® may therefore have some sort of lambda-point at an extremely 
low temperature—perhaps at about 0-2° K., the temperature at which 
the phonon model for He* probably breaks down (Atkins and Osborne 
1950)—as compared to that of Het at 2-19° K.—the temperature of the 
Bose-Hinstein condensation phenomenon. 

The velocity of second sound in a mixture of He® and Het will also 
tend to u//3 at very low temperatures, where 1, is the velocity of first 
sound in the mixture. This disagrees with Pomeranchuk’s (1949) 
prediction that the value of wu, in mixtures should be very small near 
absolute zero. To obtain his result, Pomeranchuk has used the relation 
(pn/p)r-o= eum, Whexe ¢ is the fraction of He® atoms, p the effective mass 
of a roton, and m the actual mass of an Het atom. This amounts to just 
the same thing as assuming that all the H® atoms contribute to the 
normal fluid (their effective mass being 4), which could hardly remain 
a valid approximation at the very low temperatures for which any 
description in terms of normal and superfluid atoms becomes quite 
inadequate. 
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Berkeley Cyclotron. Cosmic Ray Protons and Mesons of 5 to 50 MeV. 


By K. Gorrstzin* and J. H. Munvey, 
H. H. Wills Physical Laboratory, University of Bristol}. 


[Received May 30, 1951.] 


SUMMARY. 


The first section describes a determination of the scattering constant 
by measuring the multiple scattering of 336 MeV. protons in Ilford 
G5 nuclear research emulsions which had been exposed in the 184” 
Berkeley cyclotron. 

Then follows an account of an experimental determination of the 
relation between the scattering parameter and residual range for protons 
and mesons of the cosmic radiation. A second value of the scattering | 
constant was found, using the range and energy of »-mesons from the 
decay of 7-mesons stopping in the emulsion. The result was used to 
obtain a range energy relation for protons of energy up to 200 MeV. 


§ A. Protons oF 336 MeV. 


Inrorp G5 emulsions 400 » thick were dropped through the proton beam 
of the 184” Berkeley cyclotron. This method was chosen because it reduced 
the density of tracks in the emulsion and so made easier the identification 
of individual tracks. The beam was highly collimated and its energy 
was given as 340-+1 MeV. No tracks were employed which differed 
noticeably in direction from that of the main flux. The tracks were inclined 
at an angle of 0-5° to the plane of the emulsion and entered at the surface. 


Experimental Procedure. 


Scattering measurements were made with a Cooke, Troughton and 
Simms M4000 type microscope, using x95 objective and x 15 eyepieces. 
The ‘ coordinate method” as described by Fowler (1950) was applied. 
Details of the technique have been discussed by Menon et al. (1951). 
To minimize the effects of emulsion distortion, the measurements were 
confined to regions of the plate at least 1 em. away from any edge. The 
total length of track measured was 20 cm., and no individual track was 
less than 8 mm. in length. 


* On leave of absence from the Max-Planck Institut fiir Physik, Gottingen. 
+ Communicated by Professor C. F. Powell, F.R.S. 
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Results and Discussion. 
Second differences D were obtained for cell-sizes of 200, 400, 600, | 
800 and 1000 microns. Fig. 1 shows the distribution of |D|s from 
330 non-overlapping cells for a cell-size of 600 p. The distribution function, 
according to Moliére’s theory (Gottstein et al. 1951, referred to as I.), 
and a normal curve are drawn for comparison. Both have been 


Fig. 1. 


¢ ID lin SCALE DIVISIONS 


Distribution of second difference, from 330 non-overlapping cells. A normal 
curve and Moliére-distribution curve -—-— are given for comparison 


d:D, (no cut-off), 

Cc: Dz (corresponding cut-off : C), 

a: Ds; (corresponding cut-off : A), 

6: D for overlapping cells applying “‘ single value cut-off ” (B), 
g standard deviation of normal curve. 


normalized to the same number of | D|s ; the standard deviation of the 
normal curve is indicated by arrow “g”, 
determined in three different ways : 


(1) the mean value of all |D|s was taken. 


The mean value was 


Result : D,=1-70. 


This includes a number of rather | 

rhis inc é arge values (see fig. 1) of |D 
ane if included in a sample of limited size, would vield os ies 
a value for the mean, thus giving rise to a systematic error in the 
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estimation of the energy. It has been found convenient in previous 
work (Goldschmidt-Clermont et al. 1948, Fowler 1950) to apply a 
“4D cut-off’. Although in the present case the sample is fairly large, 
we have applied the same “ cut-off” in the following two ways, in order 
to study its influence on the scattering-constant K,. 


(2) The mean was taken of all | D[s remaining after removal of those 
values, |D|>4D,. This was done by successive approximation. 


Result : D,=1-55. 


(3) A value of | D|>4D occurring in a cell of the smallest size employed 
(200 ») indicates the position of a large angle scattering. In the larger 
cells all the values of | D | related to this ‘“‘ single scattering ’’ were removed. 
This is equivalent to “ cutting ”’ a track at a position of a single scattering 
and measuring the scattering of each side independently. In the case 
of overlapping cells the effect is to remove a block of Ds from the column 
of values as recorded, and so we refer to it as a “ block ” cut-off. 


The method yielding D, is called “ single-value cut-off”. 


Result (for overlapping cells) : D,=1-45. 


These values, D,, D,, D; are shown in fig. 1 with the corresponding 
“ cut-offs ”’. 

The “ block cut-off ”’ gives a smaller D than the “ single value cut-off”. 
This is because the former tends to exclude from the statistics the values 
of |D| which, being just below 4D, would therefore not have been 
eliminated by a “single value cut-off”. In the case of a small sample 
(short track) a “ block cut-off’? seems to be preferable to a “single 
value cut-off ’’, since it eliminates more rigorously such values of | D| as 
are well above the average and so might falsify the result if included. 

A mean has also been calculated for the distribution of 990 | D|s from 
overlapping 600, cells. This was obtained using a single value cut-off. 
The mean is indicated in fig. 1 by the arrow “b”’. 

For 600 » cells the “ noise ’’ was found to be negligible compared with 
the “true scattering ”, so noise elimination was not needed (fig. 3, I.). 
(The microscope used here was M40135.) 

In fig. 2, D2 is shown plotted as a function of s*. Within experimental 
error it is a straight line. This indicates that distortion of the emulsion 
is negligible, since the presence of distortion would have given this curve 
a pronounced upwardly concave shape (unpublished observations). 

The energy lost by the protons as they passed through the emulsion 
was taken into account, and the mean energy corresponding to the 
multiple scattering measured was calculated to be 336 MeV. The values 
of the scattering constant (as defined by I., eqn. (12)) obtained from the 
total distribution for 600 » cells are 


(a) with single value cut-off 29:2-+1, 
(6) without cut off........ 30-721. 
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Energies were measured in MeV. and cell-sizes in units of 100 He The 
comparison of these results with Moliére’s theory has been discussed 
ingle 

§ B. Prorons AND MESONS OF FROM 5 To 50 MeV. 
Experimental Procedure. 


Ilford G5 nuclear research emulsions were exposed to the cosmic 
radiation during a high-altitude balloon flight. We selected for 
investigation nineteen tracks with total length of 20 cm. ending in the 
emulsion ; each track was at least 8 mm. long. 


Fig. 2. 


10 
10090 
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DIVISIONS ) 


(SCALE 


Variation of D3 with (cell-size)’. 


Tracks produced by protons were identified by measuring their grain 
density and scattering (Fowler 1950). Tracks of particles giving rise to 
typical co-stars, or suffering u-decay were accepted as those of 7-mesons ; 
v-mesons were identified by the characteristic electron decay at the end 
of their range. 

Of the 19 tracks selected, 11 were protons, 4 were z-mesons and 
4 w-mesons. Scattering measurements were made as described in $A 
with a primary cell-size of 40 4. As before, no tracks lying near the edges 
of the plates were used. 

Results. 


Second differences D were obtained for cell-sizes of 40 and 80 microns. 
Protons 7- and .4-mesons were considered as three distinct groups and 
|D| was calculated for both cell-sizes as a function of range. This was 
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done by dividing each track into sections of 1:4 mm. length, starting from 
the point at which the particle came to rest. Then, adopting a single 
value cut-off, the mean of all | D|s belonging to sections of eee residual 
range was found. 

All emulsions used in this investigation were processed together. 
Measurements of the curvature of steeply dipping tracks indicated that 
the contribution of distortion to the scattering observed was negligible. 
However, with these small cell-sizes, the contribution of ‘ noise ’’ cannot 
be neglected. We have applied a method of noise elimination assuming 
noise to be independent of cell-size, and neglecting the variation of the 
L-term (see I.) with cell-size. If noise be taken as proportional to 
S14 (Fowler, private communication), any error in « due to these two 
approximations is, in this case, less than 2 per cent. 


Thus the equation used for obtaining the “true” value of D for each 


residual range was 
ee? 
T) (D.—D, 
Dawn [Se ME ees eh) 


(comp. eqn. (3)) 
(Menon et al. 1951.) 
From this « was calculated. 

In fig. 3, log « is shown as a function of log R. Here R is the effective 
residual range obtained when the energy loss by the particle along the 
track is taken into account. In our case it differs very slightly from the 
range at the centre of each section, and is given by 


0-15. (Ry—R,)] °*?, 
R=| ane | ery cae (2) 


where R, and R, are the ranges of the extreme points of the section 
considered. This is obtained from the range energy relation for protons 


OR Lots wee em te) | (3) 


where C is a constant, the numerical value of which is not required here, 
and the equation 


a Bae “aR 
= eer ee had oak ee) 
ae R, (Ra—Ry) Ri): 


where «7 is the mean multiple scattering parameter observed for the 
section, and « is taken to be inversely proportional to the kinetic energy. 

The values of « and R obtained from the z- and »-meson tracks have been 
replaced in fig. 3 by their equivalents for a proton of the same velocity, 
assuming the masses of the mesons to be 


m,=210+4, (Smith et al. 1950, 
m,=276+46 Lederman et al. 1951.) 


in units of the electron rest mass. 
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The theoretical formula relating log % with log R is that of a curve 
which deviates from a straight line by less than 2 per cent in the energy 
interval considered, as follows from eqn. (12) (I.) and the range-energy 
relation. This deviation is smaller than the statistical spread of the 
experimental points and so a straight line was fitted as a convenient 


approximation. 


Determination of the Scattering Constant. 


Using the decay scheme 
a7—>p-+ neutrino (Lattes et al. 1947, Marshak 
and Bethe 1947, 
O’Ceallaigh 1950.) 
Fig. 3. 
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anise in MILLIMETRES 
a as a function of range for protons. The measurements were made on protons 
(...), m-mesons and p-mesons. Instead of the values obtained for 
the mesons their equivalents for a proton of equal velocity have been 
plotted (0 00 for 7-mesons, x for x.-mesons). 


and the above values of the meson masses, the kinetic energy of a 4-meson 
arising from the decay of a z-meson at rest can be shown to be 4:1 MeV. 
The corresponding range of a p-meson in the plates used for this 
investigation we found to be 595-45 microns. These values of energy 
and range were transformed to their equivalents for a proton of the same 
velocity, and the « corresponding to the transformed range was obtained 
from fig. 38. Then (with z=1 and s in units of 100 Lt) 


K,=«(100 p)pv=26-1-+0-7 
for cell-size of 80 1 and B=0-26, This value of K, is included in fig. 1 (1.). 


MeV 


ENERGY in 
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Range Energy Relation. 


K,, as determined in this manner, is very insensitive to variations in 
the value of the exponent in eqn. (3), which was used only for the 
application of a very small correction to the range R. From fig. 3 an 
independent range-energy relation for protons may be derived in regions 
which have not yet been covered by range-energy measurements in the 
photographic emulsion. Fig. 4 shows the extrapolated curve of Bradner 
et al. (1950), together with our experimental points. The uppermost 
point (at 336 MeV.) was obtained from the measurements described in 
§ A, making use of fig. 3. 

The agreement between the curve of Bradner et al. and our experimental 
points is within the limits of our statistical errors. 


Fig. 4. 


RANGE in MILLIMETRES 


. Range-energy relation as derived from fig. 3. Notation is the same as in 


fig.3. The divided circle denotes the point obtained from the experiment 
described in §A (protons of 336MeV.). The full line gives the 
extrapolated range-energy relation of Bradner et al. (1950), the upright 
crosses denoting experimental points of these authors. 
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Department of Natural Philosophy, The University of Glasgow*. 


[Received June 11, 1951.] 


ABSTRACT. 


An account is given of B- and y-ray energy measurements. Various 
coincidence techniques appropriate to scintillation detectors are examined 
experimentally, and an investigation is made of the application of a single 
erystal, with source buried inside, to the elucidation and_ closer 
examination of decay schemes. Work is described on standard sources, 
on Hg”, and in some detail on La’. New information is given for 
La™® concerning the 0-093 y-ray, the high energy y-ray and the decay 
scheme. 


§ 1. INTRODUCTION. 


THE detection of individual particles, or quanta, by means of the 
combination of phosphor and photo-sensitive electron multiplier has been 
the subject of considerable research and development since it was first 
used by Curran and Baker (1944) for «-particles. As in the case of other 
instruments of fundamental importance much of this early research has 
been of a rather tentative and exploratory nature. The instrument 
lends itself to studies of this type since it contains at least three major 
elements which may be varied over wide ranges of performance. Thus 
we have :— 

(1) The fluorescent crystal ; obviously several criteria of performance 
can be exhaustively studied such as the conversion efficiency 
of the kinetic energy of the particle, or particles, to light of suitable 
frequency, the transparency of the crystal to such light, together 
with many other factors such as excitation and decay times, ease 
of growth, ete. 

(2) The photo-sensitive layer ; clearly the efficiency of conversion of 
the quanta excited in the chosen fluorescent substances to useful 
photo-electrons is the most important consideration, but this has 
to be regarded closely in relation to the emission of thermal electrons 
from the cathode, since such emission determines, to a considerable 
extent the ultimate sensitivity of the device, particularly for 
particles of low energy. 

(3) Ihe multiplier ; the chief objective here is to secure uniformly 
high and stable gain per stage. 


* Communicated by the Authors. 
SER. 7, VOL. 42, NO. 333-— OCT. 1951 4¥ 
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The extensive researches made in recent years with the scintillation 
detector show that the best commercial photo-multipliers satisfy largely 
the requirements (2) and (3) and gradually the data acquired on (1) is 
clarifying the situation. Thus, at present, many workers with solid 
crystals have chosen to regard naphthalene, anthracene and stilbene 
as standards among the organic materials, while sodium iodide, thallium 
activated, is preferred among the inorganic substances. In the present 
work we report results obtained with crystals of anthracene and activated 
sodium iodide, particularly in the field of y-ray spectroscopy. It is in 
this application that the detector perhaps acquires its greatest significance, 
owing to the fact that its high absorption gives it a very high efficiency. 
We discuss our results in the light of the possibility of using large volumes 
of certain liquids containing phosphors as the fluorescent substances. 
Tt is shown that, in certain circumstances, new methods of studying the 
decay schemes of radioactive substances are practicable. 


§ 2. Gamma Ray SPECTROSCOPY. 

y-rays passing through substances are absorbed or scattered mainly 
in the following processes :— 

(i) Photoelectric Effect—The quantum ejects an electron from one of 
the K, L, M, . .. shells. There is a high probability of conversion of 
the energy of the X-ray quantum excited in this process to kinetic energy 
of electrons, and a single peak may be observed corresponding to the full 
energy of the y-quantum. To secure maximum sensitivity it is obviously 
useful to employ crystals containing elements of high Z value such as 
iodine in Nal (Tl). Ifthen the photo-electric process greatly predominates, 
as in the region below hv~-:25 MeV., say, we can use a single crystal 
(~1 cm.*) as a spectrometer. For greater hy values the second process 
becomes evident :— 

(ii) The Compton Process—The scattered quantum and the recoil 
electron share the energy of the incident quantum, and a large range of 
energy values of the recoil electron results. However, the electrons 
projected in a definite direction relative to the direction of incidence are 
homogeneous. This principle has been made the basis of a method of 
measuring the quantum energy by Hofstadter and McIntyre (1950). 
We shall consider some of the features of this method in detail below. 

(ili) The Production of Pairs.—The process of creation of pairs has a 
threshold at an energy hy=2 m ¢?, and it is not appreciably evident in 
light materials up to hyv~2 MeV. It offers an attractive method of 
measuring quantum energies at high hv values, since the pair of electrons 
have constant total energy hv—2 m 2. Using the fact that the 
annihilation of a positron leads to the emission of two oppositely directed 
quanta of energy hyv=m,c?, Johansson (1950) selected pulses due to 
electron pairs from a background of pulses, and so secured good resolution. 
However, the straightforward exploitation of the process of pair 
formation requires crystals containing a high Z value element, large 
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enough to confine a high percentage of electrons entirely within the 
volume and to trap the bremsstrahlung emitted in flight. This is a 
formidable requirement. At present it would appear useful to explore 
the possibility of replacing the two arrays of Geiger Miiller counters in 
the method of magnetic spectrography due to Walker and McDaniel 
(1947) with two large crystals, or liquid counters, from which the separate 
light signals due to the positrons and negatrons could be piped to a pair 
of multipliers, adjusted to give equal output pulses for equal input 
signals. This is indicated in fig. 1. The two output voltages would be 
fed to a fast coincidence circuit and then simply added in amplitude. 
In this way the rather elaborate circuitry required to correlate all the 
pair distributions would be avoided. The magnetic field would be set 
at a suitably high fixed value. It need not be uniform or regulated since 
it is a mere particle separator. A pulse analyser would directly give the 
spectrum of the added pulses and hence of the y-radiation. The crystal 
or liquid counters need to be thick enough to stop the incident electrons 
and wide enough to satisfy the solid angle requirements of the pair 
spectrometer. 


Fig. 1. 
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‘S 
Pair spectrograph with scintillation detectors. 


Anthracene has been shown by Hopkins (1950) to give a light output 
to the multiplier proportional to the energy of incident electron in the 
region 100 KeV. to 1 MeV., but deviating slightly from linearity at lower 
energies. Activated sodium iodide has been investigated by West, 
Meyerhof and Hofstadter (1951) and has shown a linear response in the 
range 2 KeV.—500 KeV. We can corroborate the linear response of 
activated sodium iodide in the range 20 KeV.—24 MeV., from experiments 
listed below, and other experiments. 


§ 3, ExpERIMENTAL Resutts aND Discussion oF WorRK 
WITH SINGLE CRYSTALS. 
I. y-ray and B-ray Energy Measurem-nts (In*™ and La), 
Spectroscopy with single crystals offers many outstanding advantages, 
provided that the resolution can be made sufficiently good. This is 
particularly true if collimation of the y-ray beam is not required. In - 


4F2 
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this event the crystal can be used as a spectrometer of very large solid 
angle, and y—y and B—y experiments can be performed, with two 
scintillation detectors, in which particular y-rays can be related to specific 
B-particles and y-rays. This is, however, not readily achieved with 
complex sources and higher energy y-rays. 

The results which are shown in figs. 2 and 3 refer to work with crystals 
of anthracene (2 mm. x 4 mm.”) and Nal (TI) (1 cm.*) with the radiations 
of In!“ and La”. The crystals were coated with paraffin and mounted 
on an E.M.I. tube 5311. The crystal of Nal (Tl) was covered with an 
aluminium reflector. Thus, in fig. 2 we have the peak, due to photo- 
absorption, superimposed on a Compton background, as observed in 


Fig. 2.5% 
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Photo-peak due to ~22 KeV. X-rays as observed in thin anthracene exposed 
to radiations of In14, 


anthracene with an uncollimated In™* source. The X-rays emitted 
after the internal conversion of the y-ray of energy 190 KeV. show as 
a photo-electric peak in the range 22 to 25 KeV. as expected. 

In fig. 3 we show the results obtained with a single crystal of Nal (T1) 
and a source of La. More will be said of the y-rays from this source 
below, but the figure shows in some detail the results for the rays of 
energy 0-335 MeV. The Compton edge and photo peak are both evident. 
Actually, a Compton edge due to a hard y-ray (0:49 MeV.) is noticeable, 
though obscured by the 0-335 MeV. photo-peak. For this complex 
source of radiation the resolution achieved is fairly satisfactory for some 
purposes. The last peak of fig. 3 is due to hard y-rays saturating the’ 
amplifier. 
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It has been shown by Beach, Peacock and Wilkinson (1949) that 
La should emit y-rays of energy 0-093, 0-335, 0-49, 0-82, 1-62, 2-5 MeV., 
but they did not directly observe, by means of a photo peak, the y-ray 
of energy 2:5 MeV. and the evidence for the 93 KeV. radiation is not well 
founded. The decay scheme proposed by them is sketched in fig. 4. 


Fig. 3. 
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Fig. 4. 


Decay scheme of La’®, 
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Here we discuss the radiation of low energy. A number of separate 
investigations were made, and the extension of the curve of fig. 3 shows 
a rather smooth variation over the range below and above 90 KeV. 

From the nearly 100 per cent efficiency of the detector for such radiation 
we estimate that y-radiation of energy ~90 KeV. cannot be emitted in 
more than 15 per cent of the disintegrations. This matter is discussed 
further in the double crystal method below. Since, according to the 
proposed decay scheme, the majority of the disintegrations result in 
transitions of this energy, a thin crystal (1 mm.) of anthracene was used 
to detect the conversion electrons from a very thin source of La™® (several 
ugm./em.2). The results of the analysis are shown in fig. 5. A peak 
of very low intensity was observed of energy 55+5 KeV. ; presumably 
corresponding to internal conversion in the K shell of a y-ray of energy 


Fig. 5. 
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Conversion elect ons, energy ~55 KeV., due to weak y-ray, energy 90-95 KeV., 
as observed in thin anthracene exposed to electrons from La!?, 


about 90-95 KeV. The intensity is found to correspond to about 
| per 100 disintegrations and seems to establish the necessity of modifying 
the decay scheme. Using a 1 cm. crystal of anthracene the B-spectrum 
was found to extend up to the region of 2-26 MeV. A weak counting 
rate in excess of background was found up to ~2-4 MeV., but this could | 
be due to accidental pile up of lower energy pulses which are relatively 
very numerous. ‘The K conversion peak was detected for the 0-335 MeV. 
y-line, but that from the 0-49 MeV. y-ray was not strong enough to 
appear, in this particular run. . 

Il. Coincidence Measurements (using La), 


i is frequently necessary in the investigation of decay schemes to 
determine which yTays are in cascade with one another. By placing 
a source (with the B-particles screened off) between two scintillation 
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counters, the simultaneous arrival of y-rays in the counters can be 
recorded in a coincidence unit. The pulses arriving at either counter 
can, after amplification, be displayed on a cathode ray tube and the 
brightness of the display can be made negligibly weak unless a coincidence 
occurs. In such circumstances determination of the pulse height 
distribution from either counter on the cathode ray tube gives information 
concerning the y-rays which are in coincidence with one another. For 
instance, if only two y-rays appear in the distribution, the two y-rays 
are evidently in cascade. If several y-rays appear, it is desirable to 
determine specifically how they are associated. The output pulse 
distribution of one of the counters can be passed through a single-channel 
“ kicksorter ” so that only pulses associated with one particular y-ray 
are passed on to the coincidence unit from that counter. The y-rays in 
coincidence with this one y-ray can then have their pulse heights displayed. 
In this method care has to be taken to keep random coincidences low. 
This latter technique can be employed advantageously for B-y processes 
also, if the decay scheme is not too complex. 

B-y and y-y coincidences with La!° have been reported by Osborne and 
Peacock (1946), Mitchell, Langer and Brown (1947), Mandeville and 
Scherb (1948), but without any direct evidence being obtained for the 
energies of the y-rays in coincidence. 

y-y coincidences have been investigated for La, using anthracene 
as the material of the crystal in the counter connected to the display 
unit. In the first case no attempt was made to select specific y-rays in 
the other counter. The pulse height distribution curve in the anthracene 
crystal at coincidence is shown at (a) in fig. 6. The response ‘of the 
crystal (where all y-ray pulses are recorded whether in coincidence or 
not) is shown for comparison at (b). The 1-62, 0-82, 0-49 MeV. lines appear 
strongly in the coincidence curve. 

In the next case the second counter was connected to a discriminator, 
and only pulses greater than 0-8 MeV. were allowed to pass on to the 
coincidence unit. The curve so obtained is shown at (a) in fig. 7, the 
Compton type curve giving the normal single crystal response to y-rays 
from the source (whether in coincidence or not) is shown for comparison. 
- No counts were observed above 0-8 MeV. now, while below this energy 
the coincidence curve follows the same trend as the Compton type curve. 
The coincidence rate was too low for any detail in the curve to be seen. 

The 0-49 and 0-82 MeV. y-rays are, therefore, in coincidence with the 
higher energy 1-62 MeV. line. It was similarly shown that the 0-335 MeV. 
y-line was in coincidence with the 1-62 MeV. line. 


Ill. Integrating Method (using Hg?®*). 

As already indicated, one of the great advantages of single crystal 
spectrometry is to be found in the analysis of decay schemes, by means 
of coincidence experiments. The use of crystals makes it possible to 
introduce new methods of observing the inter-relationships of the particles 
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and quanta. Thus, suppose we first consider the type of de-excitation 
process indicated in fig. 8, where two y-quanta occur in cascade. With 
a source at some distance from the crystal we can observe the relative 
intensities of the y-radiations by taking a histogram. Thus, if the 
transition of energy hv, be of negligible intensity no corresponding peak 
is observed. However, if the source be placed very near the crystal the 
number of coincidences of the cascade quanta will increase, and a peak 
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coincidence control. 


corresponding to the sum of the energies hy,+-hv., or hy, will be observed. 
This will establish that hv, and hy, are in coincidence and if the source 
be Placed within the crystal (between two flat-faced blocks, say) the 
relative intensities of the peaks due to hv,, hv, and hv sie can be 
used to determine the number of decay processes that are 4 eoincidentes 
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Fig. 8. 


De-excitation process from two excited states. 
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if the y-rays are detected with 100 per cent efficiency a single 
Merits to fe ai hv,+hvgis obtained. If, on the other hand, the number 
of quanta of energy hv, exceeds those of energy hy, (for example in 
complex B-decay) then two peaks, one corresponding to hv quanta (not 
in coincidence with hy,) and one to jv; and hv, hv will be obtained. 
In this way complex schemes may be readily investigated. . Variation 
of the resolving time of the coincidence circuit can be introduced to see 
whether any of the states have measurable life-times. 

Above we considered y-rays only, foil shielding the source when B-rays 
were emitted. Fig. 9 serves as an example of a decay giving rise to 
B-y coincidences. Since it is typical of many disintegration schemes, we 
will indicate the possibilities of the method by applying it to Hg?®. 
Saxon (1948) and Slatis and Siegbahn (1949) have shown that Hg? decays 
with the emission of a f-ray of energy 208 KeV. and a y-ray of energy 
278 KeV. This y-ray is internally converted in the K and L shells 
(~20 per cent and 5 per cent respectively), giving rise to photo-electrons 
of energy ~200 KeV. and ~270 KeV. respectively. 


Fig. 9. 
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De-excitation following B-emission. 


A thin source (~1 pgm./em.?) of HgS was evaporated on to an aluminium 
foil of thickness 1-8 mg./em.?, and this was placed between two blocks of 
Nal (TI) each 1 em.x 1 em.x} cm. A single multiplier, type 5311, was 
arranged symmetrically with respect to this assembly. The efficiency 
of each of the two halves of the crystal was more than 99 per cent for all 
X-radiations following internal conversion of the y-radiation. Thus, 
conversion in the K, L, M, . . . shell, giving rise to photo-electrons of 
different energies, resulted in the same amount of electron energy release 
in the crystal. 

Hence the integrating effect of the crystal is such that we expect 
mainly two types of pulses :— 

(1) B-rays—those accompanied by y-rays which are not counted and 

which escape from the counter. 

(2) B-rays—those accompanied by photo-electrons. These photo- 

electrons can arise in several ways, so that this group can contain :— 
(i.) B-particle+L conversion electron+-photo-electrons released in 
the crystal by capture of X-radiation. 
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(ii.) B-particle+-K conversion electron+ photo-electrons released by 
the capture of X-rays. 
(iii.) B-particle+photo-electrons arising from photo-absorption of 
y-rays in the crystal. 
These two main groups should thus result in two spectra, one due to 
particles of the normal £-spectrum ending at 208 KeV. and another due 
to particles of the 8-spectrum carried forward by photo-electrons, so that 
the group starts at 278 KeV. These two groups are clearly separated. 
The second displays clearly the form from zero energy of the particle 
right up to the end point, lying on the energy scale between 278 KeV. 
and 486 KeV., where its form can be examined readily over the whole 
range. In the main, the group (2) particles of fig. 10 show the advantage 


Fig. 10. 
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Spectrum of Hg? as observed for source clamped between 
two blocks of Nal (Tl). 


of the crystal in stopping all electrons or hard X-rays within the volume 
of the detectors. Furthermore, the method can be applied very widely 
when high efficiency of detection of high energy quanta can be achieved. 
Even with the fairly small crystal used in the above work with Hg, the 
y-ray quanta are detected with an efficiency of ~50 per cent. With still 
larger crystals almost 100 per cent absorption of the radiation could be 
expected. We see from fig. 10 that the Compton process 1s appreciable 
and it results, unfortunately, in a filling-in of the gap between 
group (1) and group (2) particles. It would, therefore, be very 
advantageous if the Compton process could be made to contribute 
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usefully ; this can be done only if scattered quanta are not allowed to 
escape, so that a pulse of definite amplitude corresponding to the full 
energy of the incident y-radiation arises from the multiple Compton 
scattering process. 

Clearly, it is not at all easy to satisfy the last requirement except with 
very large crystals of good optical quality. Liquid phosphors would 
appear to offer greater promise. With a sufficiently large flask of liquid 
the escape of scattered quanta could be made negligible for incident 
y-rays of several MeV. energy. However, the efficient collection of light 
from a volume of some litres introduces considerable difficulties, and a 
possible method of overcoming this obstacle is being explored. In this 
method a large bath of the phosphor is placed in close proximity to an 
efficient large high voltage photo-cell. The anode of the photo-cell is 
covered with an efficient phosphor such as ZnS(Ag), so that the energy 
of the photo-electrons accelerated to this anode is converted again to 
light quanta, which are then piped efficiently to the cathode of an orthodox 
commercial photomultiplier. The device is actually a type of head 
amplifier. With an acceleration voltage of 10 kV. across the unit, which 
is in effect a simplified image tube, a stage gain yield of quanta of about 
4 can be expected. It has been shown by Krisek and Vand (1946) 
that at 20 kV. a unit similar in principle to the above has achieved a gain 
of 20 in light flux. A system like this seems to offer many advantages, 
but has one disadvantage, e. g. that efficient phosphors such as ZnS(Ag) are 
relatively slow, and the major advantage of rapid response is lost. 


§ 4. DouBLE CrystaL MrerHop or MEAsuRING ENERGY 
(ustne La}4°, Co®, ThC”). 


This method was tested with several sources of gamma radiation, 
among them La!°, Co®, ThC” and some new information was obtained 
concerning La1*°, 

This experimental method, based on the Compton effect, is made 
possible by the fact that backward scattered quanta do not differ much 
in energy. The energies of electrons corresponding to quanta scattered 
at angles of 90°, 120°, 150° and 180° for various values of incident quantum 
energy are tabulated :— 


0-5 MeV 1 MeV 2 MeV 
90° 25 67 1:59 
120° 30 “75 Megas 
150° 32 ‘78 1-76 
180° 33 80 1-78 


It is seen that there is a good degree of homogeneity, and within the range 
120 to 150°, which applies to most of the work described here, the 
theoretical resolving power is ~3 per cent. This is the main component 
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of the total width of the lines as determined by this method. Obviously 
the ideal geometry would be with the source placed between two crystals, 
utilizing an angle of scatter ~180°, where the uncertainty in energy is 
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Fig. 11. 


Two-crystal Compton method of y-ray analysis. 
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Gamma ray spectrum for La™° with 1 cm.’ of anthracene as scattering crystal. 
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least and the intensity greatest. This arrangement is generally not 
practicable for complex sources because of the y-y coincidences due to 
the emission of time-related quanta. The arrangement shown in fig. 11; 
using collimation, works extremely well. For maximum sensitivity 
the crystal 2 should be toroid shaped since this does not interfere with 
resolving power. However, in practice, crystal 2 consisted of a 
3-inch diameter crystal of Nal (Tl), crystal 1 was a l-cm. cube of 
anthracene or of Nal (Tl). The lead collimator was 4 inches in length and 
such that the y-rays incident on crystal 1 lay in a cone of semi-angle 3°. 


Fig. 13. 
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Gamma ray spectrum for La! with 1 em.’ of Nal (Tl) as scattering crystal, 


It is proposed to replace crystal 2 with a specially formed vessel filled 
with terphenyl in xylene. Results for the different sources are shown 
in figs. 12, 13. Figs. 12 and 13 were obtained with a source of La“ using 
crystals of anthracene and sodium iodide (TI) as scattering eryethls and 
using Co® and ThO” as calibration sources. Peaks were obtained 
corresponding to y-rays of energy 1-62, 0-82, 0:49 and 0-335 MeV., two 
of which are shown in these figures. 

The weak high energy y-line was also found, which by. comparison 
with the ThC" line (2-62 MeV.) had an energy of 2:55 MeV.+2 per cent 
We have to note that this radiation has been detected by Wattenberg 
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(1947) with an entirely different arrangement (nuclear photo effect) and 
similarly by Hanson (1949) and Bishop, Wilson and Halban (1950). 
No trace was found by this method of a y-ray of energy 0-093 MeV., 
using an anthracene crystal as scatterer. 

An assessment of the relative intensities of the y-rays emitted by 
La? as measured by the double crystal method can be made. These 
are listed as follows :— 


y-ray energies in MeV. Relative intensity in °% 
0-093 <u Ah 
0-335 3 
0-49 22 
0-82 16 
1-62 56 
2-55 3 


and agree in general trend with the results quoted for the higher energy 
y-rays by Rall and Wilkinson (1947). The method, though very 
effective, suffers as noted above in respect of solid angle considerations. 

As a result of the work with La! it appears that the general features 
of the decay scheme of fig. 4 are correct, but the transition giving rise 
to the soft radiation (hv=0-093 MeV.) is in error and cannot be a member 
of the main cascade. 


§ 5. ENERGY EXPENDED PER QUANTUM OF FLUORESCENT 
RADIATION. 


The use of a large crystal (or large mass of fluorescent liquid) with 
an «, 8 or y source placed within its volume is similar in principle to the 
technique of placing an « (or occasionally a §) source within a. proportional 
counter or ionization chamber. The multiplier performs the role of the 
amplifier, and the important corresponding quantities are the energy 
expenditure per fluorescent quantum, and the energy per ion pair produced 
in the gas. The former can be readily determined with sources of known 
energy or average energy. Thus if we employ a f-active material then 
the total number of f-particles emitted per sec., say No, is counted, and 
the current I, in the anode circuit of the multiplier is measured in 
electrons/second. Then we have the energy K expended per fluorescent 
quantum excited in the crystal given by 


K=N)VaeM/I, Day eR RY ERE ar ore 8) 
where V=energy or average energy of the primary particle in electron 
volts. 


a—collection efficiency for quanta at the photo surface ; 
including all effects of solid angle, transparency of crystal, 
etc. 
«=photo-efficiency of multiplier cathode. 
M=current magnification of multiplier. 
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Provided that K is nearly constant over the range V, say, the relationship 
(1) can be utilized in many investigations. Thus, calibration with 


radiation of known V, say, V’ gives us 
Nowvie— 15. 
where 6=aeM/K=constant. 
For an unknown radiation, we have 
No V’o=l55 
so that : Vie VENGIa/ Nolo MPEP ee ee es eS 


This direct procedure can be of considerable value in preliminary 
investigation of sources, since the total energy release per disintegration 
can be deduced. This information can be of considerable value in setting 
up a decay scheme, although it is limited, of course, by virtue of the 
escape of neutrinos. 

In conclusion, we wish to thank Professor P. I. Dee, F.R.S., for 
encouraging the prosecution of the work, and Mr. J. G. Balfour for his 
work in connection with the kicksorter. 
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ABSTRACT. 


In a cloud-chamber investigation at 2867m. altitude 70 V-tracks due 
to the decay of neutral V-particles have been observed. An analysis has 
been made on the basis of the two-body decay processes: V9 P+-+-- 
and V3—7+-+7- and it is concluded that the present limited data are 
consistent with these two processes. The ratio of the numbers of the 
two types of V-particles is found to be Nyo/Ny.=1-6+0-5 and accurate 
mass estimates have been obtained for 12 examples of the first type 
of decay and for eight examples of the second type. The results are 
consistent with the unique mass values (2203+12)m, and (796+27)m, 
respectively. The possibility that there are neutral secondary particles 
among the decay products cannot be excluded. However, some evidence 
is presented that seems to exclude an explanation of the observations in 
terms of two modes of decay of a heavy particle each giving a nucleon 
and two mesons. Finally, it is concluded that there are probably two 
different types of neutral V-particles. 
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$1. INTRODUCTION. 


V-sHarep tracks, usually associated with energetic nuclear interactions, 
have been found in many cloud-chamber investigations. Seriff et al. 
(1951), Bridge and Annis (1951), Astbury e¢ al. (1951) and the authors 
have found these tracks at mountain altitudes and Rochester and Butler 
(1947), McCusker and Millar (1951), Fretter (1951), Thompson et al, 
(1951) and Leighton et al. (1951) have observed them at sea-level. The 
majority of these tracks are due to the decay of hitherto unknown particles, 
which have been provisionally called V-particles. The present paper 
contains a detailed analysis of 70 V-tracks, due to the decay of neutral 
V-particles, obtained at an altitude of 2867m. on the Pic-du-Midi in the 
French Pyrenees. In a preliminary publication (Armenteros e¢ al. 1951) 
it was established that the neutral V-particles are not homogeneous and 
that at least two separate decay processes are involved. In one process 
the charged decay products are a proton and a negative 7-meson while 
in the other both charged secondaries are probably z- or p-mesons. 
These conclusions have been confirmed by Thompson ef al. and by 
Leighton et al. The latter authors have considered the possibility that 
a nucleon and two mesons may be emitted in each decay process and 
have concluded that their results can be interpreted in this way. The 
significance of this interpretation is discussed in §7 of this paper. 


§2. THe EXPERIMENTAL ARRANGEMENT. 


The cloud chamber, 28 cm. in diameter, has been used in a field of 
about 7500 gauss, produced by an 11-ton electro-magnet. The equipment 
is basically the same as that used in previous investigations of penetrating 
showers at sea-level by Barker and Butler (1951) and by earlier workers. 
Since accurate momentum- and angle-measurements are required, special 
precautions have been taken to avoid distortions in the chamber and in 
order to measure the angles azcurately the events are carefully 
reconstructed in three dimensions. 

(i.) Some important features of the arrangement. 

The apparatus was operated in a permanent stone-built laboratory, the 
air temperature being regulated to within 1° C. The cloud chamber was 
isolated thermally from the core of the electro-magnet and its temperature 
maintained constant at about 15° C., to within }°C., by a water-jacket 
fed from a thermostated supply. This close control of temperature is 
essential in order to reduce convection currents in the chamber to a 
minimum. The time interval between expansion and photography was 
0-07 sec., as short as possible consistent with adequate intensity of the 
photographic images. 


(ii.) The counter selection system. 
The first V-tracks to be discovered were associated with energetic 


nuclear interactions or penetrating showers and therefore counter systems 
which detect such interactions efficiently have since been used. Further- 
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more, owing to the discovery that the majority of the V-tracks are due 
to the decay of particles with short lifetimes, the material in which the 
interactions occur should be very dense and as near to the cloud chamber 
as possible. The six-fold counter system used on the Pic-du-Midi (fig. 1) 
can be operated by only one ionizing particle above the chamber. This 
ensures that energetic interactions, produced by protons in the lead 
absorbers immediately above, or inside, the cloud chamber, will be 
recorded. The average energy of a shower required to operate the counter 
system is greater than 101°eV. The counting rate was (7-2-40-1) per hour 
at 2867 m. altitude and it is estimated that more than half of this rate 
was due to nuclear interactions. With a resetting time of 6 minutes 
after each expansion, the rate of photography was 4-2 per hour. 


Fig. 1. 
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Side view of cloud chamber on Pic-du-Midi. 
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(iii.) Measurements of the curvature of cloud-chamber tracks. 


The cloud-chamber tracks were photographed by a pair of Leitz Hektor 
lenses (focal length 2:8 cm. ; aperture f/6-3) on 35 mm. Ilford 5G91 film, 
at a mean magnification of 0-1. Several methods have been used for 
measuring the curvatures of the track images, the choice of method 
depending upon the magnitude of the curvatures. For curvatures smaller 
than 1m.-! on the film the prism compensator, devised by Blackett 
(1937), was used. For greater curvatures, the track images have been 
measured with a travelling microscope or have been projected on to a 
screen and compared with a set of standard circular arcs. 


AG2 
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The minimum-detectable curvature on a counter-controlled track is a 
function of the length of the track, for given condensation conditions, 
chamber illumination and specific ionization of the particle giving rise 
to the track. It has been found that for good track images of particles 
at minimum ionization, longer than about 1 cm., the limit is set by the 
curvature due to gaseous distortions. The curvature produced by chamber 
distortions has been determined from an examination of .-meson tracks, 
10 cm. long in the chamber, which were photographed in zero magnetic 
field. The probable curvature due to distortion derived from these 
measurements was 0-2 m.~! and corresponds to a maximum-detectable 
momentum of 10! eV./c., for a mean magnification of 0-1 and an average 
magnetic field of 7500 gauss. For track images shorter than 1 cm., the 
random deviations of the drop images from the true trajectory become 
important. Thus, for example, the maximum-detectable momentum for 
track images 0-5 cm. long is reduced to about 7x 10%eV./c. A lead plate 
across the centre of the chamber was used during part of the running time 
and slight distortions sometimes occurred close to the plate. Sections of a 
track that are slightly distorted in this way can be rejected when making 
curvature measurements with the prism compensator. 


(iv.) Estimates of momenta and the measurement of angles. 

In order to derive the momentum of a particle from the curvature of 
its track image, average values are found for the magnetic field over the 
trajectory and for the magnification of the track on the film. These 
values are obtained from calibration measurements, the position of the 
track being determined by reprojecting the stereoscopic photographs 
through a pair of Leitz Hektor lenses, nearly identical optically with those 
used to photograph the chamber. A correction is made for the uniform 
dilation of the track due to the expansion of the cloud-chamber gas. 

The angle between the branches of a V-track is also obtained from 
reprojection measurements, the accuracy of the results depending on the 
magnitude of the angle and the track curvatures. The accuracy obtained 
in the present series of measurements was usually better than 5 per cent 
but, in some cases for which the angle is small and the plane of the fork - 
is inclined appreciably to the plane of the chamber, the error may be as 
high as 10 per cent. 


§3. An ANALYSIS OF EVENTS OCCURRING IN THE Gas 
OF THE CLOUD CHAMBER. 
About 105 V-shaped tracks and six events in which more than two 
tracks diverge from a point in the gas have been found among 13,500 
photographs obtained on the Pic-du-Midi. 


(i.) Nuclear interactions in the gas. 


Of the six events in the gas involving more than two ionizing particles, 
three are typical of nuclear interactions produced by neutrons of moderate 
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energies and one is an energetic penetrating shower produced in the gas 
just below the top wall of the cloud chamber. The remaining two events, 
each consisting of three energetic ionizing particles, seem to be examples 
of specially simple types of interactions. The simplest interpretations 
are that they represent the following interactions : ; 


Ne baad Nem tae ees Pe eg T) 
and 


Tage ere ee Coes Seal. wea eae G2) 


It is extremely unlikely that both of these interactions, due to charged 
particles, involved collisions with free hydrogen nuclei in the gas of the 
cloud chamber, the geometric cross-section of which is less than 5 per cent 
of that of the argon nuclei present. If heavy nuclei took part in the 
collisions the emission of evaporation particles would be expected, 
whereas none are observed. It is possible, however, that both of the events 
were associated with recoil tracks less than 1 mm. in length. The Y-shaped 
event, which can be interpreted by scheme (1), is reproduced in Pl. X X XIII. 

An estimate of the number of interactions in the gas of the chamber 
can be made in a variety of different ways and is of the same order of 
magnitude as the observed number. Hence, it is concluded that not 
more than one or two of the V-tracks discussed in (ii.) below can be due 
to nuclear interactions. A detailed discussion has been given by 
Armenteros e¢ al. (1951). 


(ii.) V-tracks observed in the gas of the cloud chamber. 


The stereoscopic pairs of photographs of the V-tracks have been 
carefully reprojected to ensure that, in each case, the apex of the event 
is within the illuminated region of the chamber. This procedure is 
essential because V-tracks are quite frequently simulated by two tracks 
crossing very close to the edge of the light beam. The angle of each 
V-track, namely the angle between the two ionizing tracks, has been 
measured and the distribution of the angles is shown in fig. 2. The 
majority of the events have angles smaller than 90°, while the remainder 
have angles close to 180°. In general, their orientation with respect to 
accompanying tracks suggests that the V-tracks in the first group are 
due to neutral particles and that those in the second group are due to 
charged particles. 

There are about 20 V-tracks, with angles smaller than 1°, which are 
characteristic of electron pairs produced by photons and most of them 
contain identifiable electron tracks. It is probable that a few electron 
pairs have been missed during the examination of the photographs; the 
smallness of the observed number, however, is sufficient to exclude the 
possibility that any of the V-events were due to to nuclear interactions of 
photons. The ratio of the cross-sections for nuclear interaction and for 
pair production by energetic photons is less than 10-*. 
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The distribution of the observed angles cannot be explained by the 
scattering of charged particles. Furthermore, only one V-track has a 
recoil track at the apex, although an argon recoil track would be 1 mm. 
long and just visible even if its momentum were as low as 5X 107 eV./c. 

The majority of the V-tracks must be due to the spontaneous decay 
of unstable particles, provisionally called V-particles, and the remainder 
of this paper is concerned with an analysis of the 70 V-tracks produced 
by neutral V-particles. 


Fig. 2. 
7 y Hiei 
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The distribution of the angle between the two tracks forming a V-track. 


§4. THE PROPERTIES OF THE DEcAY PRopucTS oF NEUTRAL V-PARTIOLES. 


The properties, such as ionizing- and penetrating-power, of the decay 
products of neutral V-particles must be determined before the nature of 
the decay processes and of the unstable particles can be elucidated. 


(i.) Measurement of heavily-ionizing secondary particles. 


Sixteen secondary particles produced heavily-ionizing tracks in the 
cloud chamber, eight of which were positively charged and eight 
negatively charged. Estimates of the specific ionization and of the 
momentum can be obtained for nine of these tracks; the data enable 
good estimates to be made of the masses of the particles. The ionizations 
were found by a procedure which has been discussed in detail by Butler 
et al. (1950). Five of the slow positive particles have estimated masses 
close to that of the proton; the tracks of two of these particles are 


Properties of Neutral V-Particles 1119 


reproduced in Pls. XX XIV. and XX XV. and the data for all five are given in 
Table I. Probably the remaining three slow positive particles were also 
protons. No identifiable negative protons have been found. No evidence 
has been found for the existence of slow secondary particles of mass 
about 1200m,. It is probable that the heavy particles found by 
Leighton et al. (1951) were protons, despite the rather low mass estimates 
in several cases. 


Taste I. 
Momentum and estimated ionization for five proton secondaries. 


Event no. 33 43 47 56 92 


Momentum 

(10% eV. /e:) 2:-4240-2| 4-8+0-5 |3:-5-+0-35| 4:50-+0-45] 6-0+0-6 
Kstimated ionization 

(Inin) 8-15 3-4 4-6 oe 2-3 


Mass range (m,) 1550-2400) 1700—2400)1450—2300)| 1600-2250 |1450-2500 


Calculated ionization 
for proton (Imin)  - 9-3 ol 5-0 3°3 2:3 


The four negative particles that have been measured have masses close 
to those of the 7- and p-mesons, the data are given in Table II. The 
measurements are not sufficiently accurate to permit discrimination 
between 7- and p-mesons. A study of the dynamics (cf. §5) of the four 
events shows that the positive secondary particles, which accompany the 
four slow these negative particles, were probably protons. The 
remaining four negative particles, on which accurate measurements 
could not be made, were all less massive than a proton. It may be 
noticed that no slow positive mesons have been observed among the 
secondary particles. 


Taste II. , 


Momentum and estimated ionization for four meson secondaries. 


Event no. 


re | a 


Estimated ionization 


(Tin) 2-3 2-4 3-4 2-5-3-5 
Mass range (m,) 270-430 210-380 290-370 230-330 


Calculated ionization 


for 7-meson. (Ipin) 2-0 2-7 2°8 3:0 
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Seriff et al. (1950) have identified one slow secondary particle as a 
meson. Hopper and Biswas (1951) found that one of the tracks forming 
a V-track in a photographic emulsion was due to a proton, and Thompson 
et al. (1950) have identified two protons among the secondaries of a sample 
of neutral V-particles. 

The behaviour of one of the slow negative secondaries observed on the 
Pic-du-Midi is of particular interest. A neutral V-track, shown in 
Pl. XXXVL., was produced just above the thin lead plate and the negative 
secondary particle was just heavily ionizing above the plate and much 
more heavily ionizing below the 7mm. plate. Five centimetres below 
the plate the track of this particle is apparently deflected through 23° 
in the gas, without any marked change in ionization. 


TABLE IIT. 


Momentum and estimated ionization for the decay 
of a negative secondary particle. 

Track No. on Pl. XXXVI. Track (2) | Track (3) 

Momentum (108 eV./c.) 7:8+0:4 | 6-1-L0:3 


Estimated ionization (I;pjn) 3-4 3-4 


Calculated ionization for 7-meson 2-8 4-0 


Calculated ionization for u.-meson 2-0 2-8 


From the measurements on the track, summarized in Table III., it is 
concluded that the change of direction of the track is most probably due 
to the decay-in-flight of a negative 7-meson. The angle of deflection is 
close to the most probable value for the disintegration of a 7-meson 
having the measured momentum. Furthermore, it is estimated that, if 
all the observed low-momentum negative secondaries of neutral V-particles 
were 7-mesons, the probability of having observed one decay is about 
25 per cent. Since the probability of a large-angle single scattering in 
the gas is extremely small, this event is regarded as convincing evidence 


for the existence of negative 7-mesons among the decay products of 
neutral V-particles. 


(ti). The momentum spectra of the secondary particles. 


Momentum measurements have been made on 43 negative secondary 
particles and on 34 positive particles; the frequency distributions are 
shown in fig. 3. A striking feature of these spectra is that the average 
momentum of the negative particles is markedly lower than that of the 
positive particles ; this disparity indicates that a considerable fraction of 
the neutral V-events has positive secondaries which are much heavier than 
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the negative secondaries (cf. §5). The contributions of heavily-ionizing 
particles are indicated on the histograms by shadéd areas. Although no 
heavily-ionizing positive mesons have been found, there are positive 
particles at approximately minimum ionization in the momentum range 
below 7x 108 eV./c. where protons would appear heavily-ionizing. This 
indicates that positive mesons are present, since it is unlikely that there: 
are any electrons among the secondary particles (cf. section (iii.) below) 


Fig, 3: 
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Frequency distribution of measured momenta of secondary particles. 


Some of these mesons which produced tracks at minimum ionization could 
have been of 1000m, type, but for the present this possibility is excluded 
and they are assumed to have been either 7- or 4-mesons. 


(iii.) The penetrating properties of the secondary particles. 
During part of the experiment a lead plate was placed across the centre 


of the chamber in order to examine the penetrating properties of the 
constituents of the V-tracks. In all, the secondary particles traversed 
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20 em. of lead without any evidence of radiative energy loss ; on the other 
hand two probable intetactions occurred. This is in good agreement with 
the results of Seriff et al. (1950) and suggests that electrons do not occur 
among the secondary particles. 

Pl. XX XVII. shows an interaction by a negative secondary particle, of 
measured momentum 3:5 108 eV./c., in the 2cm. plate with the pro- 
duction of a slow proton. The angle of the fork is 3° and the possibility 
that the V-track was due to an electron pair can be excluded. It is 
probable that the particle responsible for the interaction was a negative 
a-meson. The positive secondary had a considerably higher momentum 
and was probably a proton. The second interaction, produced by a 
positive secondary, is reproduced in Pl. XX XVIII. ; this secondary particle 
was deflected through 15° in traversing the 2 cm. lead plate. Unfortunately 
the momentum of the particle cannot be measured accurately, but it is 
probable that it was greater than 5x 108 eV./c. The momentum of the 
negative secondary was (3:3-0-3) x 108 eV./c. ; it was probably a meson. 
Since the momentum of the secondary particles transverse to the direction 
of the decaying particle was greater than 1-6 108 eV./c. and using the 
arguments outlined in §5, section (i.), it is probable that both secondaries 
were mesons (cf. fig. (6)). 

A tentative value of the cross-section for the nuclear interactions of 
the secondaries of V-particles, on the basis of the four interactions found — 
in the present work and in the work of Seriff et al., is about 10 cm. of 
lead, whereas the geometrical cross-section of the lead nucleus corresponds 
to 15cm. of lead. This supports the conclusion that at least a considerable 
fraction of the secondary particles, produced by the decay of neutral 
V-particles, have strong nuclear interactions. The observations, however, 
are not yet sufficiently extensive to exclude the presence of an appreciable 
number of .-mesons among the secondary particles. 


(iv.) Neutral secondary particles. 


No conclusive evidence has yet been reported for the existence of 
neutral particles among the decay products. The most direct way of 
testing for their presence would be to find a discrepancy between the 
direction of the V-particle, as deduced from its point of origin, and the 
direction of the vector sum of the momenta of the two secondary particles. 
Thus, the conservation of momentum, transverse to the plane of 
the V-track, may be examined by a coplanarity test and that in the plane 
of the V-track, by testing whether the transverse momentum components 
of the secondary particles are equal. 

In the present experiment, the majority of the V-particles originated 
in thick lead absorbers above the cloud chamber and it has not been possible, 
in any of these examples, to locate the origin accurately. Nine V-particles 
originated in nuclear interactions in a lead plate placed across the middle 
of the chamber, and, for two of these examples, reasonably accurate 
measurements have been made; one is reproduced in Pl. XXXIX. The 
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angles between the branches of these V-tracks were 46° and 77° re- 
spectively and coplanarity was established to within the experimental 
error, estimated at about 5°. The distances traversed before decay were 
41cm. and 1-9 cm. respectively. The momentum balance in the plane 
of the V-particle, perpendicular to the line of flight, was also examined. 
The differences in the transverse momentum components, for the two 
examples, were (0-1+0-3)x108eV./c. and (0-4-+0-2)x 10% eV./c. re- 
spectively, and thus there is no evidence for the existence of a third neutral 
particle. In each of these cases the secondary particles were both mesons. 
Another method of testing for the presence of a third particle is to look 
for experimental evidence of interactions in the lead plate, due to neutral 
decay products. Photons or neutral 7-mesons, which decay extremely 
rapidly into two photons, can initiate electromagnetic cascades, and 
neutrons can give rise to nuclear interactions. Four photographs, in 
each of which a V-particle decayed immediately above a 0-5 cm. lead plate 
across the centre of the chamber and both the ionizing secondary particles 
penetrated the plate, have been carefully examined. Two of these events 
contained low-momentum negative secondaries, and, in the other two, the 
secondary particles were of comparablemomentum. No evidence has been 
found of clearly associated electromagnetic radiation (cf. Pl. XX XVI.). 


§5. Tor Masszs oF NEUTRAL V-PARTICLES. 


Since protons and both positive- and negative-mesons have been found 
among the decay products of neutral V-particles, it must be concluded 
that there are at least two decay processes involved. Assuming that no 
third neutral particle is produced, the most likely two-body processes 
are : 

Vo > Pt+ta- Se Wa a ec AO) 
and 
Ve see ee ss ge in (4) 


Some evidence has been found for 7-mesons accompanied by protons 
but the nature of the mesons in scheme (4) is uncertain; they may be 
p-mesons. Schemes (3) and (4) may be modified by the presence of 
neutral particles, but no direct evidence of their existence has yet been 
discovered. In the present analysis the consistency of the experimental 
data with the assumed two-body processes (3) and (4) will be examined. 

The analysis of the data is complicated because many of the measured 
V-tracks have both secondary tracks at minimum ionization and therefore 
they cannot be directly assigned to one or other of the two schemes. 
Some events, however, can be directly assigned to either scheme (3) or 
to scheme (4) and, using the measurements on these events, a procedure 
has been devised for separating statistically the remaining data between 
the assumed schemes. This classification procedure is outlined in 
section (i.) below and the data for each decay process are examined in 
detail in sections (ii.) and (iii.). 
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(i.) The classification of V-tracks. 


The decay-in-flight of a V°-particle, mass M, into two secondary 
particles with masses m, and m_ is shown schematically in fig. 4. If the 
momentum of each secondary particle in the rest system of the unstable 
particle is denoted by p* and if the secondary particles are emitted in 
the rest system at an angle 6* to the direction of the V°-particle, then the 
transverse momentum component is invariant and equal to p* sin 6*. 
Further, it can be shown that, for random emission of the secondary 
particles in the rest system, the probability that the transverse component 
lies between py, and p,+dpy is 


i) = poh EE Pics: d; . 
V°-particle oe Ee 
Pepe 


7_sp.» B. 


Laboratory system. Rest system. 


Thus, for a given decay process, the probability that the secondary particles — 
have transverse momentum components close to p* is very high. The 
experimental distribution of », for a number of V-tracks can be used to 
check that the directions of emission in the rest system are at random. 

It is useful to define a quantity « given by : 


gases p+ —p_ cos p_ 
D4 C08 $y, FP 008 f.’ 
where p,cos¢, and p_cos¢_ are the longitudinal momentum com- 


ponents of the positive- and negative-secondary particles (cf. fig. 4). It 
can be shown that 


9 


ey. 2 2 
ral fe Pa Mm 


Leeruclate 
0 Pp ye +2p* cos 6* ae = mt , 
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where M and P are respectively the mass and momentum of the V-particle. 
For simplicity, the last equation may be written in the form 


a=%+f(P, cos 6*) 


and thus the longitudinal momentum components of the secondary 
particles can be related simply to the invariant quantity o, where 
%y—=(m2_—m? )/M?, because the mean value of f(P, cos 6*), for a number of 
decay events randomly oriented in their rest systems, is zero, 7. e. the 
mean value of « equals x». 

Values of « and py are used when deciding by which of the two likely 
decay schemes a particular V-track can be interpreted. These values 
can be calculated for the 29 measured V-tracks; the quantities p,, p_ 
and ¢ are known for each of these events and the value of p, can be 
obtained conveniently from the relation pp=p,p_ sin ¢/P. 

Five of the 29 V-tracks have slow proton secondaries, which have been 
identified with considerable certainty. These five events can only be 
interpreted by scheme (3); the average of four of the values of « is 
+0-60. The fifth value is very small and is neglected in the preliminary 
determination of «; it can be shown that the latter event (Catalogue 
no. 33) is an improbable one although consistent with scheme (3). The 
average value of p* is (1-10+0-10) x 108 eV./c. 

The secondaries of three of the other neutral V-particles were all 
lighter than the proton; they cannot be identified, however, because 
the ionization of their tracks is close to. the minimum value, but the 
evidence, given in §4, suggests that they were probably 7- or u-mesons. 
Thus these three V-tracks can only bé explained by decay scheme (4) and 
are due to the decay of V2-particles. . The average value of « for the three 
V-tracks is +0:-08; this value is not significantly different from zero, 
the value expected if the assumed decay scheme is the correct one, The 
average value of p* is (1-40+0-15) x 108 eV./c. 

It will be noticed that the average values of «, namely a, for the two 
schemes are different and the experimental-yalue of ‘p* derived from 
only three cases, for scheme (4) is probably greater than the value of p* 
for scheme (3). These experimental values of « and p*, for the two 
schemes, can be used to decide on the- interpretation of the remaining 
21 V-tracks. The positive-secondary tracks of these.events are sub- 
stantially at minimum ionization and, with one exception, may have been 
produced by either protons or mesons. One of the positive particles 
cannot have been a proton and was probably a,meson.. The majority of 
the negative-secondary tracks must have been produced by mesons. 

Values of « and‘p, can be calculated for the.21 V-tracks, because they 
are independent of the constants of the decay scheme. A preliminary 
classification of the events can be effected by deciding with which of the 
two values of «these values of « <are. most consistent. The 29 values 
of « and their probable errors are shown in fig. 5, where they are plotted 
against 1/P. A number cluster around the line «=0-65, which is close 
to the value of x found for four of the tracks which probably decayed 
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according to scheme (3). In addition, a number of the values cluster 
around zero, the expected value for scheme (4). 

The preliminary classification may be checked by considering the value 
of py, in addition to the value of «, foreach event. Ifa Za, then 6* ~7/2 
and therefore pp=p*. Now it has already been shown that, for a given 
scheme, the majority of the values of p, are close to p*; this result may 
be used to advantage when assigning data to one or other of the two 
schemes. For example, if «=0-7 and py,=0-2x 10% eV./c., the event 
cannot be interpreted by scheme (3) but can probably be interpreted by 
scheme (4). 


Fig. 5. 
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A detailed examination of the 29 sets of « and py reveals that 16 of the 
V-tracks can be explained by scheme (3) and ten by scheme (4); 
the interpretation of the remaining three events is uncertain. For 
example, for one of the latter events «=0:84 and p,»=0:46 x 108 eV./c 
The value of p» is small for the value of a, if the event decayed according 
to scheme (3) ; alternatively, the value of « is large for scheme (4). For 


the latter scheme p,—0 when «0-7, since P, the momentum of the 
particular V°-particle, is 1-7 x 10° eV./c. 
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The average value of « for the 16 V-tracks due to the decay of V!-particles 
is 0-65 +0-02 and is shown by a dotted line in fig. 5; the corresponding 
value of the mass of the V{-particle is (2250-35)m,. The distribution 


Fig. 6. 
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of the transverse momentum components, py, is shown in fig. 6 (a) and is 
consistent with the theoretical distribution, indicated by the curve, for 
a V-particle of mass 2200m, decaying into a proton and a 7-meson. 
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Considerably more data are required before the theoretical distribution 
can be checked in detail. ye 

The average value of « for the decay of V9-particles is + (0-05 +0-06) 
and the distribution of the values of py is shown in fig. 6 (6), where the 
curve has been calculated for a value of pp=1-5 x 108 eV./c., corresponding 
to about 800m, mass. 

The sample of 29 V-tracks is very small and therefore it is difficult to 
assess the efficiency with which their classification has been made. If P, the 
momentum of the V°-particle, is less than about 6 x 10% eV./c., the proton, 
if present, will be heavily ionizing and therefore the classification of these 
V-tracks can be effected directly. For V-particles with momenta above 
6x 108 eV./c., it can be shown in the case of V?-particles with «%)=0-65 
and p*=1-1x 10% eV./c. that their « values cannot be negative. Thus 
the total number of V$-particles, with momenta greater than 6 x 108 eV./c. 
can be assessed, within statistical fluctuation, by doubling the observed 
number of events with negative values of «. Finally, while there may 
be doubt about the interpretation of individual V-tracks with P greater 
than 6 x 108 eV./c., it can be shown that the overall separation of the data is 
at least 85 per cent efficient, provided that the individual values of « and 
zy are known to an accuracy of about 15 per cent. 

When a larger number of accurately measured V-tracks are available, 
a detailed examination could be made of the two distributions of the 
transverse momentum components obtained for the two schemes, which 
would provide a check on the classification procedure outlined above. 


(ii.) The decay process V9 > P*+-7-, 


It has been shown that 16 of the measured V-tracks were probably 
due to the decay of V{-particles. The secondary tracks of 12 of these 
have been measured accurately and the mass values calculated, using 
the laws of conservation of energy and momentum. An explicit formula 
for M, the mass of the V-particle, in terms of measurable quantities may 
be written : 


M?=m?_+-m?_+2p,p_ [t+ (m2) PP {a+ Gale —cos 4] : 
P Pp 


where m, and m_ are the masses, in energy units, of the proton and 


m-meson respectively. The energy of disintegration, or Q-value, can be 
calculated using the value of M since 


Q=M—(m,+m_). 

The values of M, Q, « and 
Table IV. 

The average mass value of the V°-particles is 2203m, and the corre- 


sponding average Q-value is 46 MeV. The errors quoted in columns (4) 
Bare Terns ‘ : : 
and (5) of Table IV. are probable errors, which were obtained by combining 


py for the 12 V{-particles are given in 
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the known distortion errors in the chamber and the various measurement 
errors. The possibility that sources of systematic error may have affected 
individual mass values cannot be excluded. For example, unexpectedly 
large gas distortion, which escaped detection, may have occurred on one 
or more photographs, since it is not possible to measure the magnitude 
of the gas distortion on an individual photograph. It is considered that 
the errors quoted must be regarded as minimum values. 


TaB_eE IV. 
Mass values of 12 V°-particles. 


(Mass of proton=1836m, ; mass of 7-meson=276m,.) 


Catalogue pr Mass value | Q-value 
no. (108 eV./c.) (m,) (MeV.) 


33* 
37 
42 


2277440 | 82+20 
2998140 | 58420 
2186 +20 
2218 +10 
2228 +10 
2160120 
2183 +20 
2181-+10 
2188 +10 
2157 +30 
2169 +10 
2256 +50 


‘78 
“94. 
‘00 
“15 
“31 
82 
“92 


* Events containing an identified proton secondary. 


The probable error of the mean mass value is 11-6m,, calculated from 
the scatter of the 12 mass values, and 8-0m, calculated from the individual 
probable errors. Thus the average mass is (2203 +12)m, and the average 
Q-value (46-6) MeV. With such a small number of events and with 
individual probable errors of the order of 20m,, it is not possible to 
decide whether the events give a unique mass value, or whether there 
is a significant distribution of values. For example, photograph 88 
apparently gives a low Q-value, namely (23 +15) MeV., but the difference 
(23-16) between this value and the average is not significant, particularly 
as the error quoted is a probable error. Similarly, photograph 33 
apparently gives a high Q-value of 74 MeV., but again the difference 
between this value and the average is hardly significant. Leighton et al. 
(1951) have obtained seven Q-values from events which probably decayed 
according to scheme (3); their average value is 49 MeV. and the individual 
values range from 32 to 65 MeV. These values are consistent with those 
given in Table IV. and, as their errors are also comparable, it is hardly 
possible to reach their conclusion that the Q-values are distributed over 
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a range of values. Thompson et al. (1951) have reported two Q-values 
of (31-L5) MeV. and (28+10) MeV. for events which, if they were two-body 
decays, probably followed scheme (3). 

(iii.) The decay process V3 > a+-+-7-. 

Eight of the ten events, which are due to the decay of V$-particles, 
have been measured reasonably accurately and their masses and Q-values 
are given in Table V. These events are less satisfactory than those 
described in section (ii.) above, since the majority involve measurements 
on high-momentum secondary particles. Most of the values in columns (4) 


TABLE V. 
Mass values of eight V$-particles. 


Vo >at+r. 


Catalogue p SS Q-value 
: (MeV.) 


796-4130 | 122+65 
888+ 50 | 165425 
872+ 50 | 160-425 
841+ 60 | 144430 
820+ 50 | 134425 
700-- 30 74.115 
673-100 61-50 

— 50 —25 
785+ 30 | 116-415 


* Events in which both secondary particles must have been less 
massive than the proton. 


and (5) of Table V. are compatible with the average mass value of 
(796+-27)m, and average Q-value of (122+13) MeV., if the secondary 
particles are 7-mesons. If the secondaries are j-mesons, then the 
corresponding average values are (705+32)m, and (142+16) MeV. 
respectively, the errors being probable errors of the means. 

Pls. XX XTX. and XL. show examples of the decay of a V$-particle 
into two mesons. In the first event, both secondaries were less massive 
than the proton and were probably mesons. In the second case, both 
secondaries could be of protonic mass, but the dynamics of the decay 
show that the event can only be interpreted by the decay of a V3$-particle. 


§6. Momentum Spectra AND RELATIVE FREQUENCIES 
oF NEUTRAL V-PARTICLES. 


An estimate has been made of the relative numbers of the two types 
of neutral V-particles that decayed in the chamber; the result obtained 
is Ny,/Ny,=1:6+0-5. There are insufficient data to decide whether 
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this ratio varies significantly with distance from the lead absorbers ; 
such a variation would, in general, be expected if there are two types of 
V-particle with different lifetimes. 

The differential momentum spectra of the V?- and V°-particles are 
given in Table VI. (a) and (b); the data for the V8-particles has been 
normalized to the observed number of V)-particles. These results can 
only be regarded as preliminary, but it appears that the average momentum 
of the V‘-particles is about the same as that of the V2-particles. A possible 
interpretation of this result is discussed in the next paragraph. 


TABLE VI. 


Differential momentum spectra of V°-particles. 


Momentum range (108 eV./c.) | 0-8 | 8-16 | 16-24 | 24-32 | 32-40 | >40 


(a) Observed spectrum of 


V°-particles 6 3 2 2 1 2 
(6) Observed. spectrum of 
V$-particles oo 6-4 1-6 1-6 — 3°2 


(c) Calculated spectrum of 
V°-particles 


$7. DiIscussION: 


The preliminary data that have been presented in §5 are consistent 
with the existence of two types of neutral V-particle, each with a unique 
mass. This interpretation suggests that there are two entirely different 
neutral V-particles which must possess similar mean lifetimes, since 
comparable numbers of the two types decay in the cloud chamber. 
The evidence for the existence of neutral V-particles with mass greater 
than the proton is conclusive but the existence of neutral mesons of mass. 
about 800m, is less certain. It may be noted that the latter particles 
could be the neutral counterparts of the charged particles with masses of 
the same order of magnitude. The evidence in favour of the heavy 
charged particles is already very strong and it is probable that the 
majority of the charged V-particles observed in the cloud chamber are 
of this type. 

It is possible to explain the observed neutral V-tracks as different 
modes of decay of the same type of particle if these modes involve three 
secondary particles; for example, the following symmetrical schemes 
may be possible : 


ve Bag Paes Re tee, aves eee | 5 eee ee ie 8 (9) 
VE > NO-bar® sad. ee ra BY 
4H2 
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Leighton et al. (1951) have pointed out that the apparent Q-value of 
process (5) (calculated from data on the two charged secondaries only) is 
probably about half the apparent Q-value of process (6) and that the 
latter value should be roughly equal to the real Q of the three-body 
disintegration. This follows, since in the rest system of a V°-particle all 
three secondary particles will, on the average, have the same momenta, 
but the two lighter particles will take nearly all of the available kinetic 
energy. Therefore the apparent Q of process (6) is about double that of 
process (5). In fact, the measured Q of scheme (4) is (122+13) MeV. 
while that of scheme (3) is (46.6) MeV.; thus the observed ratio is 2-6. 
Thus the present data are not in contradiction with the prediction made 
by Leighton et al. and would be consistent with a neutral V-particle mass 
of approximately 2600m,. No evidence has, however, been found for 
the existence of either neutrons or neutral 7-mesons among the secondaries 
of neutral V-particles. 

However, the momentum spectra shown in Table VI. (a) and (b) provide 
some evidence against the explanation of the neutral V-tracks by the 
two symmetrical three-body decay schemes (5) and (6). The spectrum 
of the particles giving proton and meson secondaries, namely that of 
V{-particles, appears to include a few more low momentum examples 
than does the spectrum of the V$-particles, which decay into two mesons. 
This is just the reverse of the expected behaviour if the decays really 
followed the three-body schemes (5) and (6) because it can be shown 
that, in the great majority of cases, the nucleon has a considerably higher 
momentum in the laboratory frame of reference than that of either of 
the mesons. The momentum of the nucleon is the largest component of 
the momentum of the V-particle ; this is true for schemes (3), (5) and (6). 
Thus the spectrum of the V2-particles of scheme (6), omitting the contribu- 
tion of the neutrons, should include many more low momentum cases than 
does the spectrum of V{-particles of scheme (3), or of scheme (5) omitting 
the contribution of the neutral 7-meson. 

On the average, the two mesons in the symmetrical schemes (5) and (6) 
will have about the same momenta. The spectrum of the negative 
m-mesons can be obtained from the data on the V°-particles and if this 
spectrum is doubled an approximate spectrum for the V$-particles is 
obtained, of course, without the contribution of the neutrons. This 
calculated spectrum is given in Table VI. (c) and may be compared with 
the observed spectrum in Table VI. (b). Ifthe two symmetrical schemes (5) 
and (6) are the correct ones, then it appears that most of the V$-particles 
should be of much lower momentum than those observed. There seems to 
be no obvious explanation of the absence of low momentum V-tracks of 
the type interpreted as due to the decay of V$-particles ; their apparent 
absence can hardly be explained as being due to instrumental or selection 
bias. It is concluded that the observed spectra of the V°-particles are 


a in agreement with the consequences of the symmetrical schemes (5) 
and (6). 
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The above argument, based on the difference between the spectra of 
the V{- and V}-particles, assuming only two secondary bodies, may, of 
course, be invalidated if there is more than one type of neutral V-particle 
giving two charged mesons as decay products. 


§8. CONCLUSIONS. 


The main conclusions reached from the analysis of a small sample of 
70 V°-tracks are as follows : 


(i.) Slow protons and slow negative mesons have been identified among 
the secondary particles. These negative mesons were probably 7-mesons. 
and the corresponding positive secondaries were probably protons. 

(ii.) Several positive secondary particles, less massive than the proton, 
have been observed. In three of these cases the corresponding negative: 
particles were much lighter than the proton and were probably mesons. 

(iii.) The neutral V-tracks can be interpreted by two different schemes, 
each involving two particles, namely 

V9 > Ptta- 
and. 
Vo > at +a- (or V3 >pt+tp-). 

(iv.) No evidence has been found that suggests the presence of neutral 
secondary particles. 

(v.) Unique mass values for the two decay processes have been found ; 
they are (2203+12)m, for the V-particles and (796+27)m, or 
(705 +32)m, for the V$-particles. 

(vi.) Rather more of the first type of decay than of the second type 
have been found, but the mean lifetimes appear to be comparable. 

(vii.) An alternative explanation of the above results can be made by 
assuming two symmetrical three-body processes each involving a nucleon 
and two mesons, but as yet there is no direct evidence in favour of this 
hypothesis. Preliminary investigations have been made on the differential 
momentum spectra of the V?- and V$-particles and the results do not 
appear to be in agreement with the symmetrical three-body schemes. 
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DESCRIPTION OF THE PLATES. 


PLaTE XXXII. 
A nuclear interaction in the gas. 


‘The interaction was produced below the plate by particle (1) ; tracks (3) and (4) 
were produced, along with a heavy blob of ionization, which was probably 
due to a nuclear recoil. Track (1) cannot be measured because it is partially 
obscured by a low-energy 6-ray. The positive track (2) has a measured 
momentum of 3-5 x 108 eV./c., and an ionization of (4-5) times minimum ; 
it was probably a proton. The positive track (3) has a measured momentum 
of 6:5 x 10% eV./c. and is inclined at an angle of about 7° to.the plane 
defined by tracks (1) and (2). Momentum is apparently not conserved in 
a direction at right angles to the incident particle ; thus the event is not 
an elastic scattering. The simplest scheme appears to be : 


P++P+—> P++N0+a4, 


PLaTE XXXIV. 
An identified proton secondary from a V9-particle (Catalogue no. 47). 

‘The data for the positive track (1) are given in Table I. ; the positive particle 
was probably a proton. The momentum of the negative particle was 
(1-42 -+0-10) x 108 eV./c. and the angle of the fork was 87-5°. The mass 
of the V{-particle was (2228+10)m, (cf. Table VI.) 


PLATE XXXV. 
An identified proton from a V°-particle (Catalogue no. 56). 
The data for the positive track (1) are given in Table I. ; the positive particle 
was probably a proton. The momentum of the negative particle was 


(1-62 +-0-10) x 108% eV./c. and the angle of the fork was 44°. The mass of 
the V}-particle was (2181-++10)m, (cf. Table IV.) 
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PuaTe XXXVI. 
The decay-in-flight of the negative secondary of a V9-particle. 


This photograph is described in § 4, section (i.); the data for the negative 
secondary particle are given in Table III. Tracks (4) and (5) form a second 
V-track; the measured momenta are approximately +3108 and 
—7x 10% eV./c. and the angle of the fork is about 10°. This fork was 
probably produced by a V$-particle. 


Puate XXXVII. 
A nuclear interaction produced by the negative secondary of a V°-particle. 


This photograph is described in § 4, section (iii.). The negative secondary 
particle (2) had a momentum of about 3:5 108 eV./c. and interacted in 
the plate with the production of a slow proton. It is probable that the ~ 
particle responsible for the interaction was a negative 7-meson. 


Pirate XX XVIII. 
Two penetrating secondary particles. 


This photograph is described in § 4, section (iii.). The positive particle (1) was 
scattered through 15° in the plate and the negative particle (2) has a 
momentum of 3-3x108 eV./c. It can be shown that both secondaries 
were probably z-mesons. ; 


Pratt XX XIX, 
The decay of a V$-particle (Catalogue no. 66). 


The V°-particle was produced in a secondary nuclear interaction in the 7 mm. 
lead plate. The momentum of the positive secondary particle (1) was 
measured as 2°6x10% eV./c. and that of the negative secondary (2) as 
1:3108% eV./c. Both secondary tracks are approximately at minimum 
ionization and must have been due to mesons. The mass of the V}-particle 
is (700-+30)m,, assuming that the secondaries were 7-mesons (cf. Table V.). 


Puate XL, 
The decay of a high-energy V$-particle (Catalogue no. 53). 


The fork consists of two high-energy tracks (1) and (2) at an angle of 12°. 
The measured momentum of the positive track is 1-5 10° eV./c. and that 
of the negative track (2) is 1-6x10°®eV./c. For this event a=—0-03 and 
pr=1-66 and it is concluded that the V°-particle probably decayed into 
two -mesons (scheme (4)); the mass value of the V$-particle is 
(841-L60)m, (cf. Table V.). 
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ABSTRACT. 


An automatic Wilson cloud chamber with specially constructed 
cameras was used to obtain data on the nuclear scattering in argon of 
electrons from indium 114 and positrons from copper 62. The tracks 
were analysed by the method of projected angles. The experimental 
angular scattering distribution was compared with calculated values 
based on Dirac’s relativistic electron theory, the ratio expt/theory over 
the angular range 20°-90° being 1-02 for electrons and 0-92 for positrons. | 
Similar calculations based on the semi-classical Rutherford scattering 
formula gave for the positron case a ratio expt/theory of 0-79. This is 
regarded as evidence in favour of the Dirac theory which gives a smaller 
scattering cross-section for positrons than for electrons. 


§ 1. INTRODUCTION. 


NuMEROUS investigations have been made with the purpose of testing 
the theory of scattering of fast electrons by atomic nuclei. The experiments 
fall, in general, into two classes ; in the first the scattering material is 
in the form of a thin foil and the scattered electrons are detected by 
ionization chambers or counters, and in the second the electrons are 
scattered in the gas of a Wilson cloud chamber. 

Both experimental techniques present considerable difficulties, and 
much of the earlier work was inconclusive and inconsistent. However, 
recent investigations of the first type by Van de Graaff et al. (1946, 1947) 
and by Buechner et al. (1947) gave good agreement with the relativistic 
theory of electron scattering for the particular energies, angular range and 
scattering elements selected. 

Randels, Chao and Crane (1940, 1945) have summarized the results 
obtained by the cloud chamber method prior to 1945, and reported an 
extensive investigation covering a number of gases and a wide range of 
energy. ‘They found no evidence for the large discrepancies from the 
theoretical predictions reported by some earlier investigators. 

The theoretical cross-section for the nuclear scattering of high energy 
electrons has been derived by Mott (1929) in the form of a conditionally 
convergent infinite series. The first approximation yields the ‘‘ Mott 
formula ” which is applicable only to elements of low atomic number. 


* Now at the Glasgow Royal Cancer Hospital Research Department, 
+ Communicated by the Authors, 
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Bartlett and Watson (1940) have summed the series for the case of 
mercury (Z—80), and McKinley and Feshbach (1948) have extended the 
calculation to obtain the cross-section for scattering of fast electrons by 
all elements to an accuracy of 1 per cent. The results are given in terms 
of the ratio of the calculated scattering cross-section to the Rutherford 
cross-section (modified for the variation of mass with energy), namely :— 


Ze? \2 1—B2 
o=( =)  poseot = ae ee ee (LY 


2m,c* f 2 


The calculation shows that for elements of high atomic number a large 
departure from the Rutherford scattering cross-section is predicted. 

The scattering of fast positrons by nuclei is also of considerable interest 
since the theoretical angular distribution is modified appreciably by 
the change in sign of the Coulomb force, through the medium of the 
spin-orbital interaction. This results again in a departure from the 
Rutherford value which is quite different from that obtained for electrons. 
Following the procedure of Bartlett and Watson, Massey (1942) has 
calculated the cross-section for positron scattering by mercury nuclei 
for various energies. The calculated cross-section is lower than the 
Rutherford value over the whole angular range. The difference is still 
appreciable, even for elements of low atomic number and should be 
experimentally demonstrable. It may be noted that the method used 
by Van de Graaff et al. is not applicable to the positron case, since it 
involved the use of an artificially accelerated beam of electrons. 

No experiments on the single nuclear scattering of positrons appear to 
have been reported at the time of writing, though Fowler and 
Oppenheimer (1938) and also Lasich (1948) obtained evidence for an 
excess of electron over positron scattering in lead and gold respectively, 
under conditions in which the scattering was multiple, merging into a 
single scattering “ tail’ at large angles. 

The present paper reports the results of the analysis of 198 metres of 
electron and 163 metres of positron track obtained with a Wilson cloud 
chamber containing argon as the scattering element. 


§2. EXPERIMENTAL PROCEDURE. 


The cloud chamber method is particularly suitable for studying 
scattering phenomena in that a record is obtained of each individual 
event, and no difficulties occur such as ensuring conditions of single 
scattering and eliminating stray electrons which may arise from X-rays 
or by reflection from the walls of the vessel. Moreover, electron-nuclear 
and electron-electron collisions can easily be distinguished, and the 
energy before and after collision can be determined by measuring the 
curvature of the track in a known magnetic field. The main disadvantage 
lies in the difficulty of accumulating sufficient data to provide results of 
statistical reliability. 
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In undertaking the present investigation, it was decided (a) to design 
a cloud chamber and camera system which would be capable of 
automatically recording a large track length without sacrificing the high 
quality of track photography which is necessary for accurate work, and 
(b) to employ a method of analysis which would ensure an accurate 
comparison with theory without entailing a prohibitive amount of 
labour. 

The cloud chamber was of the moving rubber diaphragm type, 
30 cm. in diameter by 6-5 cm. deep. It was illuminated by a single 
Philips, xenon-filled discharge lamp which dissipated in each flash 
about 500 joules. The illuminated region of the chamber was restricted 
to a depth of 2-5 cm. by means of masks placed between the lamp and 
chamber. 

The cameras, which were specially designed and constructed for use 
in this type of work, were entirely automatic. The shutter was operated 
by a solenoid and plunger, and the film was wound on one frame after 
each exposure by an electric motor activated by a cam and microswitch 
incorporated in the camera. The cameras were mounted on a brass 
framework which fitted on to the upper ring of the cloud chamber ; the 
axis of one was vertical and the other had its axis inclined at 20° to the 
vertical. The lenses were a matched pair of 80 mm. focal length, 
f/4-5 Taylor-Taylor-Hobson “ Ental”’ enlarging lenses. Each camera 
could carry in a special cassette 100 feet of 60 mm. recording film. The 
film was, however, supplied in 25-foot lengths, covering about 
150 exposures, a convenient number for a normal run of the apparatus. 
The film was held flat in the position of exposure by a spring-loaded 
glass pressure plate. A circular hole in the back plate, covered during 
photography by a light-tight metal disc, allowed the insertion of a 
condensing lens system for use in the reprojection of tracks, described 
in a later section. The cloud chamber was mounted between a pair of 
large Helmholtz coils, each consisting of six flat coils of } in. wide copper 
strip insulated by “Armite” paper and energized from the D.C. mains. 
The whole assembly comprising cloud chamber, cameras, magnet coils, 
vacuum reservoir and associated electrical circuits for operating the 
chamber release valve, lamp delay, camera shutter, etc., was mounted on 
a trolley. A second trolley carried the vacuum pump, and the condensers 
and charging unit for the discharge lamp. 

The source of electrons was indium 114 prepared by neutron 
bombardment of a strip of indium metal in the Harwell pile. This 
isotope, which has an effective half-life of 48 days, proved a convenient 
and prolific source of electrons of energy up to 2 MeV. In order to achieve 
good collimation and some energy resolution of the electrons, 2nd to avoid 
the presence in the chamber of clumps of ionization due to the low energy 
y-tays emitted by the source, the latter was mounted on a-brass block 
and placed between the poles of a resolving magnet. The magnet current 
and the position of the source were adjusted until a suitable number of 
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tracks appeared in the chamber per expansion. Lead blocks placed at 
suitable points between the source and the chamber provided any further 
screening required. A shutter operated by a solenoid via a relay which 
was “slugged ” by a variable resistance to introduce a short time delay, 
prevented electrons from entering the chamber until the expansion 
was nearly complete, thus reducing the number of diffuse tracks in the 
photographs. The electrons entered the chamber through five portholes 
consisting of thin mica discs attached to brass rings which were waxed 
into holes drilled in the glass cylinder of the chamber. 

The taking of photographs was entirely automatic, the whole process 
being initiated by pressing a button on the control panel. The sequence 
of events, with approximate times, was as follows : 

1. Zero time : button pressed ; camera shutter opens. 

2. 1/50 sec. : expansion completed ; source shutter opens. 

3. 1/5 sec. : lamp flashes. 

4. 1 sec.: camera shutter closes; source shutter closes; slow 

expansions start ; film winds on one frame. 

5. 60 secs. : slow expansions completed. 
After an interval of 40 seconds to allow temperature equilibrium to be 
reached, the cycle is repeated. 

The selection of a suitable positron source presented some difficulty. 
Inspection of the table of radioactive isotopes shows that all those of 
suitable energy have such short half lives that irradiation in atomic 
piles is not practicable. It was decided therefore to use copper 62 which 
may be obtained by a y-n reaction with the stable isotope copper 63, the 
threshold energy being about 10 MeV. The maximum energy of the 
positron is 2-6 MeV., but the half life is only 10-5 minutes. The source 
was prepared by irradiating five in. side, 1/16 in. thick copper squares 
for 20 minutes in the beam of the Glasgow University synchrotron. 
After irradiation, the squares were attached to an aluminium support 
and placed one in front of each of the chamber portholes. As the activity 
weakened the squares were brought nearer the chamber, and in this way 
a useful working life of 30-40 minutes was obtained per radiation. 


§ 3. ANALYSIS OF TRACKS. 


Since stereographic pairs of photographs were taken, it would have 
‘been possible to make measurements in three dimensions. It was 
decided, however, to employ the method of projected angles used by 
Randels, Chao and Crane, and further developed by O’Ceallaigh and 
MacCarthaigh (1944) and by Barker (1948). In this method the 
projections of the scattering angles and the curvatures in the plane of the 
chamber are measured, and the scattering formula is expressed in terms 
of the projections. 

The projected form of the general scattering law has been derived 
by O’Ceallaigh and MacCarthaigh in the form of a double integral. 
This involves a parameter p=A,/2a, where 2a is the depth of the 


1140 A. F. Howatson and J. R. Atkinson on the 


illuminated part of the chamber and A, is the “ critical ’ length. Deflected 
tracks shorter than A, are rejected since such tracks would not allow 
accurate measurement of scattering angles and curvatures. Barker, 
and O’Ceallaigh and MacCarthaigh show that if the minimum track 
length condition is applied, the ratio of the number of scattering events 
in a given angular range in the projected system to the three-dimensional 
value is lowest at small angles, projection having little effect at large 
angles. As interest usually centres on large angle scattering statistics, 
this is a distinct advantage as it effectively reduces the number of small 
angle scattering events. In addition, there is the enormous increase in 
the speed and accuracy of the measurement of angles and curvatures 
when only a single plane is involved. The second, oblique photograph 
is, nevertheless, most invaluable in clarifying cases where the association 
of crossing tracks is in doubt. 

Barker (1948) and O’Ceallaigh (1950) have drawn up tables for various 
values of p from which may be calculated the number of scattering events 
predicted by theory in successive intervals of scattering angle. In the 
present work the illuminated depth was 2:5 cm. The minimum track 
length was taken to be 2-5 cm. also, so that the value of p was 1. 


§ 4. MEASUREMENT OF TRACKS. 


To make measurements on the tracks the developed films were 
returned to the cameras, which, together with the supporting framework, 
were removed from the cloud chamber and placed above a table. A full- 
size image of the vertical photograph was projected on to a sheet of 
paper on which was drawn a circle 23 cm. in diameter, the latter being 
adjusted to be concentric with the image of the chamber. Tracks which 
emerged from the windows and were of good quality were traced on 
the paper with a sharp pencil. The oblique image was also examined 
but the tracks were not traced, its function being, as previously mentioned, 
to clear up doubtful cases where tracks crossed. 

Each track was numbered and the following information entered 
in a Table at the upper right-hand corner of the sheet. 

1. The length / within the limits of the circle, determined by means 

of a map measurer. 

2. The radius of curvature, obtained by finding the best fit of a series 
of curves consisting of arcs of concentric circles engraved on a 
celluloid sheet. 

3. Nuclear deflections were noted and the angle of scattering measured 
where this exceeded a certain minimum value. 

Deflections occurring outside the circle were ignored, though parts 
of tracks outside the circle were often used in measuring angles and 
curvatures relating to events occurring inside. The drawing of the circle 
serves the same purpose as the fixing of a minimum track length ; viz., 
it ensures that no scattering events are recorded of which it is impossible 
by reason of the shortness of the scattered track to measure accurately 
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the curvature and the angle of scattering. The mean radius of curvature 
and hence the mean energy of the particles was obtained by forming 
the sum 2{l/(Hp)?} and equating it to 57/(Hp4,)?. 


§ 5. Discussion or REsuLts. 


The calculation of the scattering angle distribution in the projected 
system is based on the fact that the usual scattering formule can be 
expressed in the form :— 


P(0) dd=Nf(Z, B)F(6)2n sin 0 dd, 


where P(@)d@ is the probability that a particle moving at speed Sc through 
a medium containing N nuclei of atomic number Z per unit volume will 
be scattered once through an angle lying between @ and §+dé@ in a 
path-length of 1 cm., and where F(é) is of the form :— 


cosec* 6/2 A cosec? 6/2-+B cosec? 6/2+-C cosec 6/2 
e.g. in Mott’s formula :— 
F(6)= cosec* 6/2—? cosec? 6/2-++-7ZaB(cosec® 6/2— cosec 6/2) .  .(2) 
or, in the formula derived by McKinley and Feshbach :— 
F(0)= cosec* 6/2—? cosec? 0/2-+-7ZaB(cosec® 6/2—cosec? 6/2) . (3) 


The probability of scattering into the angular range ¢ to 6+d¢ where 
¢ is the projection of # in a plane perpendicular to the axis of the camera 
(y-axis) is (O’Ceallaigh and MacCarthaigh, 1944) :— 


2Nf (Z, B) a °u  -F(@) sin 0 
ra, ceed eng) 
where cos 6 seieN 
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Barker gives values of P,,(¢)/Nf(Z, 8) corresponding to 6=10°, 20°. . . 90°, 
for F(@)=cosec” 6/2 where n=1, 2, 3 and 4, and for values of p=0, 0-5, 
1:0 and 2-0. From these figures, by graphical integration or otherwise, 
the integrated values over suitable angular ranges may be calculated. 
O’Ceallaigh (1950) presents similar Tables containing the integrated 
values calculated for the angular ranges 20°-30°, 30°-40°, etc., for a 
number of different values of p. 

The number of scattering events observed, together with the number 
predicted from theory, in angular intervals of 10° in the range 20°—90° are 
shown in Table I. The theoretical values are based on formula (3). The 
half integers occur because angles of 20°, 30°, 40°, etc., were assigned 
half to the lower and half to the upper interval. 

The agreement with theory is excellent over the range 20°-90°. ‘The 
total number of scattering events predicted is 69-4, the observed number 
70-5, giving a ratio of experiment/theory of 1-02. This is well within the 
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range of statistical fluctuations and experimental error. The distribution 
of number of events with scattering angle also agrees very well with 
theory. In particular, there is no indication of the excessive scattering 
between 25° and 55°, of which Randels, Chao and Crane claimed to have 


TABLE I. 


Scattering of electrons and positrons in argon ; comparison of 
experimental and theoretical results. 


ELECTRONS POSITRONS 
Energy Range (MeV.) 0:2—2-0 || Energy Range (MeV.) 0-2 —2-4 


Effective Energy (MeV.) 0-83 Effective Energy (MeV.) 0-7 


Theory | Theory 


Angular Range | Expt. | Theory || Expt. (1) (2) 
20°—30° 49-5 44-9 41:5 52:8 47-1 
30°—40° 13-5 13-6 15-5 16-4 13-8 
40°—50° 3-0 5:6 3:5 6:7 5:3 
50°—-60° 2:5 2-6 3-0 3-4 2:5 
60°—70° 1-0 1-4 1-0 2-0 1:3 
70°—80° 1:0 0-8 1-0 1-3 0-8 
80°-90° 0 0-5 0 0-7 0-4 

20°—90° 70-5 69-4 65-5 83:3 71-2 


Tasue II. 


Scattering of electrons in argon ; Randels, Chao and Crane (1945). 


Energy Range 


Hin Energy Range 
0-8—3-3 MeV. 


3-3—9-3 MeV. 


Angular 
Range Pee 
. Effective Energy Effective Energy 
2-4 MeV. 4-6 MeV. 
ixpt. Theory Expt. Theory 
15°-25 30-0 30:3 52-0 32-2 
25°-35 10-5 7:2 10-5 7-6 
35°45 5:0 2-7 4-5 2:8 
45°--55° 2°5 1:3 5-0 1-4 
55 —65° 1-0 0-7 1-0 Q-7 
65 ~75 0 0-4 0 0-4 
15°—85° 1-0 0-2 0 0-2 
85°-90 0 0-1 0 0-1 


— 
oo 
ete! | 
>. | = 
—_—————————| 
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obtained strong, though not conclusive, evidence from their observations 
over several gases. The results obtained for argon by Randels, Chao and 
Crane are shown for comparison in Table IT. Table IIT. gives a summary 
of earlier results of Zuber (1938), Stepanowa (1937, 1939) and Bleuler 
(1942). 

The assessment of experimental error is difficult, since the observer’s 
judgment is necessarily involved to some extent in the selection and 
measurement of tracks. It cannot be too strongly emphasized that unless 
the quality of the tracks and of the track photography is of a high 
standard serious errors in interpretation are likely to occur. The 
standard appropriate to this type of work is such that individual back- 
ground drops should be clearly distinguishable over the whole of the 
illuminated part of the cloud chamber. Since it is only within recent 
years that rare gas discharge lamps have been developed of sufficient 
output to allow the vertical photography of individual drops and thinly 
ionizing tracks over a large cloud chamber volume, it seems probable 
that the large discrepancies observed by some earlier workers are 
attributable to experimental difficulties in the photography, and 
consequent inaccuracies in the analysis of the tracks. 


Tasxe III. 


Scattering of electrons in argon; summary of earlier measurements. 


Track 

Energy Angular Number of| LHzxpt. k 

Range (MeV.)| Range Hagen Deflections | Theory oS eae 
(metres) 

1-7-2-4 30°-180° 350 48 0-75 Zuber 1938 

0-2-1-1 20°-150° 103 308 1-0 Stepanowa | 1939 

1-5-3-0 20°-150° 130 84 2-5 Stepanowa | 1939 

0:2-3:0 20°-180° 708 158 1-5 Bleuler, 1942 

Scherrer & 
Zunti. 


In the present experiments it was found that the agreement between 
theory and experiment did not extend below 20°. Scattering angles 
were measured down to 15°, but for a calculated value of 82-1 in the 
range 15°—20° the observed number was 49 : for positrons the corresponding 
figures were 97:8 and 55. This discrepancy is attributed to small angle 
scattering which is always present to some extent, and which makes the 
detection of clear-cut scattering events difficult at angles as small as 
15° or so. The effect becomes rapidly less important as the scattering 
angle increases. 

An important objective in scattering experiments is to distinguish 
between the relativistic scattering formula which is based on the Dirac 
electron theory and the semi-classical Rutherford formula. A comparison 
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of the two for the case of scattering of 1 MeV. ee se steer 
(calculated for projected angles with p=1) is shown 1n the g. It ° gene 
shows that for electrons the difference is too small in the ae a 
angular range 20°-60° to be detected by the present aie ie : e 
corresponding curve for positrons shows, however, that the wie - 
much larger, being of the order of 20 per cent in the range 20°—60". 

is feasible therefore to attempt to distinguish experimentally between 
the two theoretical predictions, even with a scatterig element of low 
atomic number such as argon. 


electrons 


positrons 


Ratio of scattering to Rutherford scattering 


“02 rs 10° 20° 30° 40° 50° 60° 70° 30° 90° 


Scattering angle. 


1 MeV. electrons and positrons: calculated ratio of scattering to Rutherford scattering 
(projected angles, p=1) in Argon. 


The experimental angular distribution for positron scattering was 
compared with the calculated values based on (1) the Rutherford formula 
with relativistic mass (R), and (2) the formula derived on the assumption 
that the positron is a positively charged Dirac particle. This may be 
written :— 

Ri{1—P? sin? 6/2—7aBZ(sin 6/2—sin? 6/2)}, 
and differs from the corresponding formula for electrons only in the 
change in sign of the third term. These results are also shown in 


Table I. The ratio experiment/theory in the angular range 20°-90° is 
for case (1) 0-79 ; for case (2) 0-92. 
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The results obtained under similar experimental conditions for electron 
scattering in argon provide a good basis for assessing the significance 
of the above figures. The authors consider that the first value is outside 
and the second inside the limits of error appropriate to the conditions of 
measurement, which are estimated to be about +10 per cent and that the 
experiments therefore provide evidence in favour of the relativistic theory 
of scattering of positrons. 


CONCLUSIONS. 


The cross-section for the nuclear scattering in argon of high energy 
electrons over the whole of the angular range 20°-90° agrees with 
theoretical predictions, but the experiments are not sufficiently sensitive 
to distinguish between different formulae. The corresponding formule 
for positron scattering are such that it is possible to distinguish 
experimentally between predictions based on semi-classical theory which 
gives the same scattering cross-section as for electrons, and on the 
relativistic theory of Dirac which gives a smaller cross-section. The 
experimental results are in agreement with the latter theory. It follows 
that the cross-section for positron scattering is in these experiments 
shown to be significantly smaller than the cross-section for electron 
scattering. 
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CXV. Micro-slip in Metal Orystals. 
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Cavendish Laboratory, Cambridge*. 


[Received June 14, 1951.] 
[Plate XLI.] 


SUMMARY. 


Aluminium crystals, the surfaces of which have been prepared for 
examination solely by electropolishing, show extremely fine slip bands 
(displacements of 150 A. or less) after small plastic deformations of about 
0-5-2 per cent. The more familiar coarse bands involving displacements 
of the order of 5000 A. appear subsequently. However, if the metal 
surface is mechanically polished prior to electrolytic polishing, the coarse 
bands occur in the earliest stages of deformation. 

Two possibilities are considered: firstly that the micro-slip is a pheno- 
menon restricted to the vicinity of the surface and is highly dependent on 
the nature of it. Secondly it is suggested that the fine slip is typical of 
initial slip processes within the metal, but becomes almost entirely 
suppressed as the deformation proceeds. 


It has become accepted that the appearance of slip bands on the surface of 
a plastically deformed metal is evidence that the deformation is not 
homogeneous, but is concentrated on relatively few atomic planes. At 
small strains and at room temperature slip band spacings are usually of the 
order of several microns. 

During some recent work on the deformation of aluminium single 
crystals it was noticed that, on specimens which had been polished 
purely electrolytically, with no preliminary mechanical preparation, the 
slip band spacing after small strains was much less than usual. Most of 
the bands were diffuse and the finer ones could not be resolved by the light 
microscope. Specimens which had been mechanically polished, however, 
showed the usual sharp, widely spaced slip bands. This suggested that 
the condition of the metal surface was responsible for the difference, so 
a detailed study has been made of the effect of surface treatment on the 
appearance of slip bands produced by subsequent deformation. 

In a typical experiment, one face of a large, flat single crystal of 99-99 
per cent aluminium was rubbed along its length and for half its width with 
coarse 150 carborundum paper and then smoothed off with successively 
finer papers. The other half of the crystal was not abraded and the 
surface remained the same as when the crystal was formed. The whole 


* Communicated by Sir Lawrence Bragg, F.R.S. 
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crystal was then electrolytically polished in the De Sy-Haemers bath (20 
per cent by volume perchloric acid and 80 per cent absolute ethyl alcohol) 
until about 60 of metal had been removed. At this stage polishing 
scratches were no longer visible. The crystal was then deformed 0-7 per 
cent in tension and examined microscopically. In the zone which had 
been mechanically polished the slip bands were sharp and widely spaced 
(20apart). There were no fine bands. On the other hand, the zone which 
had been polished purely electrolytically showed a very large number of 
closely spaced diffuse markings parallel to the coarse bands in the other 
zone. The spacing was not measurable in the light microscope. After 
further deformation to 3-3 per cent extension, the mechanically polished 
zone showed sharp slip bands about 10, apart, while in the zone which had 
not been mechanically worked there was a high density of small, diffuse slip 
markings together with a few more sharply defined bands parallel to them. 
On further deformation the number of slip bands in the mechanically 
polished zone continually increased. In the other zone, the fine markings 
predominated, but as the deformation proceeded the more sharply 
defined bands increased in number. The appearance of the two zones 
after 5 per cent extension is shown in figs. 1 and 2 (Pl. XLI.). 

After 6-8 per cent extension the crystal was repolished electrolytically to 
remove all traces of slip bands. When the straining was resumed the new 
slip bands appearing were mainly of the coarse, sharp type in both zones. 
Thus, tensile deformation of the crystal produces eventually the same 
effect as prior mechanical polishing: the very small diffuse slip bands are 
suppressed and the more familiar widely spaced bands alone appear. 
Throughout the experiment, the appearance of the underside of the 
mechanically polished zone of the crystal, which had been electropolished 
only, was the same as that of the non-mechanically polished zone of the 
upper surface. Thus the effect of mechanical polishing did not penetrate 
_to the underside of the crystal. The depth of penetration depended on 
the severity of the mechanical polishing. This was investigated by 
removal or known thicknesses of metal by electropolishing, followed by 
re-straining and examination of the new slip bands, the total strain through- 
out remaining small. The effects of the rubbing with 150 carborundum 
paper were still very marked after 70» had been removed and X-ray 
micrographs taken at this stage showed that the fresh surface was 
considerably distorted. After 100” had been removed the effects of 
the 150 carborundum still persisted but were less obvious. The coarse 
slip lines could be produced by rubbing the surface with the finest polishing 
alumina but then the effects did not penetrate more than a few microns into 
the crystal. Annealing crystals for an hour at 600° C. completely removed 
the effects of any surface polishing treatment. 

Examination of the abraded and un-abraded zones by means of the 
electron microscope gave the results shown in figs. 3 and 4 (Pl. XLT.) 
respectively. These are from the same specimen and at the same strain 
as the light micrographs in figs 1 and 2, The replica technique used was 
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as follows: methyl-methacrylate monomer was polymerized on the 
surface of the specimen and, when hard, stripped off to form een 
mediate replica. This was then shadowed with gold palladium alloy i 
a low angle to the surface and in the direction down the slip steps. 4 
this means each slip step on the surface is revealed by its shadow, and 
from measurements of this the height of the step and the slip distance can 
be calculated. For the abraded zone the slip distance is about 5000 A. per 
band, a figure in good agreement with that calculated from the measured 
strain and the observed density of slip bands. In the case of the 
un-abraded zone, slip bands can be clearly measured with slip distances less 
than 200 A. However, in this case it appears that the bands large enough 
to be observed cannot account for the whole strain: there must therefore 
be some still smaller bands, perhaps with slip distances down to atomic 
dimensions. , 

Figs. 3 and 4 are crossed by a furrow-like structure with a mean spacing 
of about 650 A. This structure is developed in the electrolytic polishing 
bath and it has been suggested that it is the traces of planes making up a 
surface of minimum energy. 

Preliminary experiments have shown that micro-slip is also observed on 
un-abraded surfaces of copper, lead and gold. 

The problem set by the above results is: which type of slip band is 
initially operating in the interior of the metal? If, as has hitherto been 
assumed, slip takes place by large amounts on isolated planes, then the 
micro-slip must be a purely surface phenomenon. In order to be visible 
on the surface, a slip process approaching from below must either create a 
new piece of surface by breaking through, or it must stretch the existing 
surface. In either case energy is required corresponding to the increase 
of surface energy, so any condition which reduces surface energy makes it 
easier for slip to appear on the surface. Severe working of the surface 
produces a high concentration of dislocations there. Thus the surface 
energy is reduced (since less work is required to remove atoms from a 
disordered surface than from a perfect one) and coarse slip results. At an 
undisturbed surface, where the surface energy has its full normal value, a 
slip process on one slip plane and involving several hundred atoms, may 
gain energy by dissociating near the surface into a number of micro-slips 
covering a relatively wide area, Some support for this view is given by the 
observation that the fine slip bands, as seen in the electron microscope, 
usually form clusters separated by regions apparently free from slip. 

On the other hand, it is possible that the fine slip bands, which are 
always observed on the surface of specimens free from previous mechanical 
distortion, are typical of initial slip proceses within the metal. Since these 
fine bands have spacings much closer to atomic dimensions than the 
slip bands observed by previous workers, the difficulty of explaining why 
slip should be nucleated in very large amounts on isolated planes is partly 
removed. Instead, it is necessary to consider how this micro-slip is blocked, 
‘for it has been shown above that blocking occurs immediately, as a result 
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of deformation of the surface, or gradually, as a result of overall tensile 
deformation of the specimen. The most probable cause of the latter is 
that dislocations on planes as closely spaced as the micro-slip interact 
strongly, preventing further movement except on a relatively small number 
of slip planes. This mechanism can also account for the high rate of 
strain hardening in the early stages of plastic deformation when the coarse 
slip bands which are usually observed, are too far apart to interact. 
Abrasion of the surface also results in the high concentration of dis- 
locations required to stop micro-slip. Annealing the specimen at a high 
temperature causes the dislocations to migrate and results in the 
Te-appearance of micro-slip. 

An alternative source of the blocking of micro-slip is the oxide film which 
is known to be particularly coherent in aluminium. Very small slip 
‘processes could appear on the surface as diffuse bulges without breaking 
the film, but, as the deformation increased the resistance of the film would 
be overcome occasionally, thus giving the opportunity for a coarse slip 
band to develop and appear as a sharp step on the surface. With a 
mechanically worked surface, the oxide is probably less coherent and the 
slip processes would then have less difficulty in penetrating it. If this 
mechanism exists, then the appearance of micro-slip should be strongly 
dependent on crystal orientation: a small slip process approaching a 
surface normally is much more likely to break a surface film than one 
approaching obliquely. Thus, when the slip planes are almost normal to 
the surface, coarse slip should be observed at an earlier stage of deforma- 
tion than when the slip planes are obliquely inclined to the surface. 
Some support for this is given by the study of deformed polycrystals. 
However, while an oxide film may assist in blocking the micro-slip 
in aluminium, it cannot do so in gold. 

Examination of the surface alone will not decide whether initially metals 
deform entirely by micro-slip or not. A sensitive technique is needed 
which is capable of revealing very fine slip bands in the interior of deformed 
crystals. 

This work is being carried out during the tenure of an I.C.I. Fellowship 
(A.F.B.) and the Armourers and Brasiers’ Research Fellowship (R.H.), 
which we gratefully acknowledge. 
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SUMMARY. 


_ This paper introduces some ideas and general theorems pertaining to 
non-linear or mechanical systems. Thus in the case of a resistive network, 
the dissipation is divided into two parts—the “content ” and 
““eo-content ’’—which are duals of each other. The dissipation itself 
has stationary properties in linear but not (in general) in non-linear 
networks, but it can be shown that the “ content”? and “ co-content ” 
shave stationary and additive properties in the non-linear case. The idea 
of “ content ” is extended to reactive systems, and it is shown that the 
total content of any system in motion is an invariant. 

This paper is intended as a companion to the paper by E. C. Cherry. 


§1. IyTRoDUCTION. 


THE mathematical analysis of non-linear systems is difficult, and often 
tends to obscure essential behaviour. An attempt is made in this and a 
companion paper to present some new concepts and general theorems 
which might assist in forming a more physical approach to the analysis 
and synthesis of systems possessing either incidental or intentional 
non-linearities. The arguments are presented in terms of electrical 
networks, but the concepts are equally applicable to mechanical and 
electro-mechanical systems. The approach here is based upon the 
methods of generalized dynamics, which in essence is a way of obtaining 
the differential equations pertaining to a dynamical system. Much 
work has already been done on the solution of non-linear differential 
equations ; here, however, attention is concentrated, not on mathematical 
methods of solution, but rather on examining the fundamental aspects 
of the subject and trying to extract useful physical concepts from them. 

One of the first theorems relating to linear networks was due to 
Maxwell (1873), the “ Minimum Heat Theorem”. (See §2.) It can be 
shown that the theorem is not applicable, as it stands, to non-linear 
networks ; the heat generated is not in general stationary, but a related 
quantity called the ‘content’, which is proportional to the heat in 
the linear case, is shown in §4 to have a stationary property. A 
summation theorem is also given, on the ‘‘ content” of two-terminal 
networks. In §6 the idea of “ content” is extended to reactive networks, 


where it is shown te be an invariant of the motion of such systems. 
ee 
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§2. Maxweti’s Stationary Heat THEOREM. 


In any purely resistive medium the current and voltage are distributed 
according to the two continuity laws: 


| div 1=0 See ter ee 
and 


p Bds=0, Pee CARO Me tae econ ed: 


The analogous conditions in a lumped circuit are given by ‘“ Kirchhoff’s 
current law” 
At=0 at a vertex Ee bohins need Psy wate. (o) 


and “‘ Kirchhoff’s voltage law ”’ 
dv=0 round any mesh. Or. xc aera hs eS) 


Maxwell’s theorem (Maxwell 1873) for lumped circuits, put in more modern 
form, states that ¢f any linear resistive circuit is driven by voltage or current 
generators or both, then of all possible distributions of current consistent with 
Kirchhoff’s current law, that which gives a minimum to the quantity 
(W—2P,) ws the only one consistent with Kirchhoff’s voltage law, and is 
therefore the actual distribution. Here W is the total power dissipated in 
the resistors, and Py is the total power taken from the voltage generators. 
Hence for a purely current-driven network the total heat is stationary. 
Notice that Kirchhoff’s current law is observed throughout, but the 
voltage law must be violated in taking the variations. 

A network with m meshes may be considered to have m degrees of 
freedom for currents. Any m non-interdependent current variables can 
therefore be chosen as generalized coordinates to define the state of the 
network; these currents could for instance be all the mesh currents, 


or some of the branch currents. If they are i,, where r—1, 2,...m, 
then Maxwell’s theorem leads to the m equations 
0 
ice — = ° . ° . . . . 5 
Ou, (W 2Py) 0, ( ) 


which are sufficient to determine all the currents (and hence all the 
voltages) in the network. 

The dual of Maxwell’s theorem is often useful, but the two theorems 
are sometimes confused. Here Kirchhoff’s voltage law is observed 
throughout, and the current law violated; the stationary quantity is 
now (W—2P,), where P,; is the power taken from the current sources. 
Then as in equation (5), if we have n generalized voltage coordinates v, 
in a network of (n-+1) vertices, we have the n equations 

0 


55 (W—2Py)=0, 2s es ee (8) 


dj 


which again tell us everything about the network. 
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Maxwell’s theorem applies only to linear networks, and it can in fact 
be shown, by the methods of § 4, that (W —2Py) is not in general stationary 
if any of the elements are non-linear. The search for some other quantity 
which does have the stationary property forms the starting-point here. 

e 


§3. Dermirions Retatine TO Non-LiInEaR NETWORKS. 


3.1. Elements.—A very general type of 2-terminal element may be 
specified by a relation of the form 


di d% dv d?v 
£(i.Gpegue s+ gpg) =0 ots gerade (7) 


where i and v are the current and voltage at the terminals. If the 
time-derivatives are absent the element is said to be non-reactive, and is 
then specified by a relation of the form f(i, v, )=0. Problems in which 
all the elements are non-reactive can be solved at any instant in terms of 
time-invariant non-reactive elements specified by a relation of the form 
fl, v)=0. Discussion of non-reactive elements is therefore confined to 
this latter case, in which elements are specified at all times by the single 
characteristic, f(7, v)=0. Such an element is said to be “ passive” if 
the characteristic cuts the i-v axes nowhere except at the origin ; 
otherwise it is said to be “active”. It is sometimes desirable to 
distinguish between “simple” and “ non-simple”’ elements, according 
as to whether or not the current and voltage are everywhere single-valued 
functions of each other. 


3.2. Generators.—A “generator” is a particular form of active 
non-reactive element in which either the current or the voltage is absent 
from the defining function f(i, v). Thus generators are either of the 
“constant voltage’ or “constant current” type. In what follows, 
expressions such as “all the elements of the network ”’ will include the 
generators uniess stated otherwise. 


3.3. Networks.—A “ passive non-reactive network” is defined to be one 
which contains only passive non-reactive elements. It is assumed that 
if two vertices of such a network are chosen as terminals, the whole 
network between the terminals may be regarded as an element with its 
own passive non-reactive characteristic. 


3.3.1. Although Thévenin’s network theorem is well known, certain 
points which have a bearing on this discussion are worth clarifying. 
This useful theorem really consists of two parts ; as usually proved they 
are (for non-reactive networks) : 

(i.) A 2-terminal linear active resistive network can always be replaced 
by a. new. 2-terminal. network consisting only of a constant-voltage 
generator E,, in series with a linear resistor R,, so far as any linear 
element outside the terminals is concerned. 
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(ii.) Further, E, is given by the open-circuit voltage at the terminals 
of the original network, and R, is given by the resistance at the terminals 
of the original network when the voltage generators are short-circuited 
and the current generators removed. 

If the load is non-linear, the theorem may be shown still to apply, 
although the usual proof for the linear case is no longer valid. If the 
network is non-linear the theorem of course does not apply ; all that can 
now be said in place of part (i.) of the theorem is that, for external purposes, 
the 2-terminal network may be replaced by a single active non-linear 
element (the ‘‘ equivalent ’’ element). 


3.4. “ Content” and “ co-content’’—Suppose that a current 7, is 
passing through an element specified by the characteristic f(i, v)=0, 
such as that shown in fig. 1(a). The voltage across it will then be »v,, 
where f(2,, ¥,)=0. We now define a quantity G, which we shall call the 
“content ”’ of the element when it is in the particular state (i,, v,), by 
the integral 


Ge { adi: 
0 


The content has the dimensions of power. The dual quantity, obtained 
by interchanging current and voltage everywhere, we shall call the 
““co-content ”, denoted by J. Thus 


= | ~ ide, 


The integrals G and J are shown as areas on fig. 1(a). If the element is 
passive (fig. 16), G and J are respectively the areas below and above the 
characteristic curve. If the total dissipation in the passive element is W, 
it can be seen from fig. 1 (b) that 


Wise = Gd wet oem tr ws» (8) 


A linear passive element is a special kind of passive element in which 
f(t, v)=Ri—v, where R is a property (the “‘ resistance ”’?) of the element. 
In this case, G—=J—W/2 for all values of current and voltage. Fig. 1 (c) 
shows the characteristic of an element which is a constant-voltage generator. 
In this case, G=—P,,, where P, is the power supplied by the generator ; 
no meaning can be attached to J. Similarly, fig. 1(d) shows a constant- 
current characteristic. Here J=—P,, while no meaning can be attached 
to G. If the elements are non-simple, G and J must be interpreted with 
care. 


§4, THEOREMS oF “‘ ConTENT”’ AND “ Co-CONTENT ans 
4.1. In this section four network theorems are given. Theorems TH. 
and IV. are the duals of I. and IL., and will be stated without proof since 
they follow from the Principle of Duality. These theorems apply at 
every instant of time. 
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Theorem aie \. The pene of a two-terminal network, regarded as 


content 
: f co-contents 
an element between the terminals, is the sum of the tye } of the 


constituent elements. 


Theorem ny \. If in an active non-reactive network the sum a of the 


tent. F 3 , 
Whale her all the constituent elements is eapressed in terms of the 


co-contents 
current 


coordinates in the network, subject 
voltage 


defining number of generalized 


only to the restrictions of Kirchhoff’s at law, then aI is stationary for 


voltage 
the actual distribution of ae 
Fig. 1 
(Z,v)=O F%,yv)=0 
(Z,,¥,) tap oars aie sa oo (21,0) 


(b) Passive element. 


F(é,v)=0 


Se oe ea ee (2,,v)) 


(c) Constant-voltage generator. (d) Constant-current generator, 


Typical characteristics of various kinds of element, showing content, G 
? 
and co-content, J, as areas. 


4,2. To prove these theorems, Kirchhoff’s two laws are used. For 
Theorem I., consider a network such as that shown in fig. 2, in which each 
branch is a non-reactive element. These branches meet in (n+-1) vertices 
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numbered Ol pies oeysnk, =. <> .%e Lake the vertices 0 and*7 as 
terminals, and let a current to, be flowing in at vertex 0 and out at vertex n. 
Let the resulting voltage at vertex j be v,;. (This use of “ vertex voltages ” 
automatically incorporates Kirchhoff’s voltage law in the proof.) Without 
loss of generality, we may put v,—0. Also let the resulting current in 
branch jk, measured from vertex j to vertex k, be i;,. Now suppose 
that %» is increased slightly by an amount 47,. All the branch currents 
will then change by an amount 47,,, and the vertex voltages by an 
amount 4v,;. It follows from Kirchhoff’s current law that 


¥ 0, jA0, JAN, 
—Aiy, j=n. j 
Fig. 2. 


mene 


rae 
Seon tt eee 


oes 
A typical network, showing vertices j, & and branch jk. 


The change AG,, in the content of each element is given, from Taylor’s 

theorem, by 

; AG 5 .=(¥j;—0;) 40; ,+441;,(40;—A0,)+> « - 
=(v;—v;,)4i;,, to first order. 
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Now since i;,=—t,;, we have 4G;,=4G;,;, and hence, if G is the sum 
of the contents of each element, we have, to the first order, 


2AG=2 2 AG;;, 
ja hae 
=2' X (v;—0 4) A1;;, 
aK 


= SiS ops yee 2 Onigy 
j 


le k 


=2 5 v;Ai,;, (from symmetry of j and k) 
jk 


Spy E z Ai. | 
j k 


=2v,4i,, from equations (9), since v,=9, 


1. €. AG=v,4i,, to first order. 
Hence a SHUge ee RR le 


But if G’ is the content of the network, regarded as an element between 
the terminals, we have 


dG’ 
ie SEU.) sop. a Pas ee 
Hence, from (10) and (11), : 
G—G’=(G),.>— (Gee =) <n 


Now it follows from the definition of content that (G’);_9>=0. If all the 
elements are passive, then (7;,);-9=0, and so (G);_»=0. Therefore 
from (12) we have G=G’, which proves Theorem I. Notice that we have 
proved it only for passive networks ; it can be shown that the theorem 
does not hold (in general) if the network is active. . 

4.3, For Theorem II., again consider such a network as that shown in 
fig. 2, but without special terminals. Any of the elements may be active. 
It will be assumed that only one element joins any two vertices, since 
any network may be reduced to this form. If there are any constant- 
current generators they are not considered as branch elements, since 
variations cannot be made in the currents through them. As before, 
we write 7;, for the actual values of the branch currents, and v, for the 
actual values of the vertex voltages. Now consider an adjacent set of 
current values (7;,+ 4i,;,); the variations 4i,, will be equivalent to 
variations in generalized current coordinates, provided that they are 
consistent with Kirchhoff’s current law, 7. e. provided that 


3 AiO tt ee ne 
k=0 
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(The adjacent current values can only obtain if Kirchhoff’s voltage law is 
violated, and hence cannot themselves be also actual values.) Let the 
variation in the content of each element be 4G,,, and in the sum of all 
the contents be 4G. 

We have 


24G=2 2 AG,, 
ne 
=2 7 (v;—v;,)4i,,, to the first order, and hence, by the 
S} 


same reasoning as used in the proof of Theorem I., 


Gas E 3 Ai, 
j k 


j 


whence AG=—0, from: (13): 


This proves the theorem. 

4.4. It should be noted that either Theorem II. or Theorem IV. is 
sufficient to determine all the currents and voltages in a network. This 
is done by using the equations corresponding to (5) and (6) in §2. They 
are 


0G 

a= 0 eS a Vn ean gE 
and 

oJ 

Pe ke er O 
where r=1, 2,....mand j=l, 2,....n, for a network of m meshes and 


(n=1) vertices.. The expressions in (14) and (15) are written as partial 
derivatives with respect to independent generalized coordinates, such as, for 


example, mesh currents or vertex voltages. (In computing the number \ 


current : current 
coordinates, constant 
voltage voltage 


any, are not counted as branches, but they appear as making contributions 


of the branch a 
voltages 
generalized coordinates. ) 


An example of the case n=1 is given in the next section. 


of generalized } generators, if 


} when the latter are expressed in terms of the 


§ 5. A Smwete PracticaL ExamMPLe oF THEOREM IV. 


This section illustrates the rather trivial one-variable case of 
Theorem IV. Consider the “ scale-of-two ”’ circuit, whose essentials are 
shown in fig. 3 (a). The right-hand portion of this circuit can be treated 
as a passive élement, which may be shown to have the characteristic 
given by the full line in the upper part of fig. 4. If now the two halves 
of the circuit are connected together, the voltage across the right-hand 


* 
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element will take up one of the “actual values”’ consistent with 
Kirchhoff’s two laws. There are three such values, as may be seen by 
drawing the “load-line” for R, shown dotted in fig. 4. (The question 


Fig. 3. 


(a) (0) 


Essentials of ‘‘ scale-of-two ”’ circuit. Schematic of fig. 3 (a). 


Fig. 4. 


Co-CONTENT 


vO 
Co-contents J, and J, of the elements in fig. 3 (b), plotted as a function of 
the voltage variable v. The sum (J,-+-J5) is also shown. 


of stability is not discussed here since it involves times; in a physical 
circuit, the stability of each of these three positions depends ‘upon 
reactive elements inevitably present.) Now consider the network to be 


¥ 
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drawn schematically as in fig. 3 (b). This is a two-element passive network 
fed by a voltage generator ; it has only one degree of freedom for voltages, 
and so, taking the voltage at B as zero, the voltage v at A may be chosen 
as variable. If, without using Kirchhoff’s current law, the co-contents 
J, and J, are plotted as functions of v, the curves shown in the lower 
half of ey. 4 are obtained. Also shown on the same figure is J, which is 
the sum of J, and J,. It will be seen that J has stationary values at, 
and only at, the oF acho values ” of v. 


§ 6. INVARIANTS OF THE MoTION FoR GENERAL Circuits INcLUDING 
REACTANCE. 


Consider now a network consisting of general elements of the type 
defined by equation (7), but with ¢ not appearing explicitly in the defining 
function /. Suppose that a complete solution has been obtained for all 
the currents and voltages in the network. Then, for each element, an 
i-v diagram may be drawn ; this diagram is not in general a characteristic 
of the element, but represents a trajectory in a type of so-called “ phase 
space’. The change in content of the element in going from the position 
(2;, V;) to the position (t,, v,) may be defined by the integral 


tg 
G.—Gy= | var . . . . . . . (16) 
and similarly 
ee 
J2—Jy= | tdv. . . . . . . . (17) 
Vy 
We are interested here in studying only the change in content as the 
currents alter with time, so we may simply talk of the “ content ”’, G, 
without concerning ourselves with the original state to which it is referred. 
Now Theorem II. as stated applies only to non-reactive elements, but 
there is nothing in the argument to exclude reactive elements, since the 
theorem depends only upon Kirchoff’s laws. If then we write G or J for 
the total content or co-content of the network (7. e. the sum of all the 
G,,’s or J;;’s), including that of any generators, we have from (14) and 
(15) 


rete) Me eee ee esa oaer ee, CLS) 


0G 
0, 
and A 
= = fd ea, ere ae, a8, 
ao OS (19) 
where r=1, 2,....m, andj=1, 2,....n. The total time-derivatives of 
G and J are given by 
dG 0G Ga, 0G di, 0G din on 
er oto Ton die on ai 20) 
and 
2 oJ a3 dv, , a dv, oJ dv, 


=; OF bey: dit ' A, bet Sy: Sas 
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Consider now any period of time during which all the voltage and current 
generators, if any, have steady voltages and currents respectively. During 
this time, t does not appear explicitly in the defining functions of any of 
the elements, and so a = a =0. Hence from equations (18), (19), 


(20) and (21), 


dG 

dt = 0 
and i 

dt = 0. 


We can state this result as follows :— 

Theorem V. In any network of time-invariant elements (possibly 
including generators), the total content and the total co-content are invariants 
of the motion. This theorem could, for instance, be applied to any circuit 
for all time following an impulsive change in a current or voltage generator, 
or in, say, a capacitor element. Note that it is very general, and applies 
even to dissipative systems in which the total energy is not invariant. 
The content and co-content, used in the above sense, are akin to integral 
invariants ; information on these is given by Whittaker (1904). 
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SUMMARY. 


The theorems and concepts presented here are generally applicable to 
circuits of electrical or mechanical elements with non-linear characteristics. 
The quantity called the co-energy (the dual of energy) is shown to possess 
stationary properties (maximum or minimum) and _ superposition 
properties ; this is sufficient to establish the concept of an “‘ equivalent 
element ”’ for any 2-terminal system of like elements (all-inductor, all- 
capacitor, all-springs, etc.). The unfamiliar “ rectangle representation ” 
of a circuit of linear resistors is explained and extended to the non-linear 
case, including reactive elements. 

It is shown that the equations of motion of a non-linear system, 
possessing reactance, may be expressed in Lagrangian form, thus 
emphasizing the importance of co-energy and also showing that the 
Principle of Duality is applicable. . 

Finally, systems are considered possessing mutual inductance and 
moving magnetic circuits (as in rotating machines). 

This paper is intended as a companion paper to the accompanying one 


by W, Millar. 


§ 1. IvrropuctTion. 


THIS paper is concerned not with the solution of non-linear differential 
equations, but rather with certain general properties of a non-linear 
physical system. In order to develop some simple, but fundamental, 
concepts and theorems, a return is made to the methods employed by 
Maxwell, J. J. Thomson and others when first deriving the basic pro- 
perties of linear lumped networks, namely the methods of general 
dynamics and the use of simple variational principles (Jeans 1948). 
Although specific reference is made here to electric circuits, it will be 
apparent that the results apply generally to other systems, also to mixed 
systems such as ‘‘ Servo ” systems. 


§ 2. Types or Non-LINEAR CHARACTERISTICS, 


It is important to distinguish between circuit element characteristics (the 
relation between the mutually dependent quantities x, y, which defines 
the element operation), which are (a) non-linear but bilateral and 


* Communicated by the Author. 
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(b) unilateral. Elements of the latter type can transmit energy one way 
only and are essentially control elements (for example, relays, triode 
valves). 

Bilateral characteristics of elements, which are the only type that are 
considered here, again divide into two classes. The first contains simple 
characteristics which here imply a single-valued relationship between the 
two variables x, y ; the second contains non-simple characteristics which 
imply a multi-valued relationship (for example a dynatron valve, an 
inductor with hysteresis loss, or gears with backlash). The theorems 
which follow will refer directly also to the non-simple class in an obvious 
manner. 


§ 3. ENERGY AND ‘“‘ Co-ENERGY ” STORAGE. 


The power dissipated in any resistor at any instant is W=v7, the area 
of the rectangle (inst. volts current), independently of the shape of 
the v/i characteristic. In a companion paper Millar has divided this area 
into two parts, lying respectively below and above the curved characteristic 
of the element, called the content, G={v di, and the co-content, J=fi dv. 


Bigoals 


(a) Inductor. (b) Capacitor. 
ae nes ste Olas , : : ; 
Characteristics (“simple ’’) of non-linear elements, showing energies T, U and 
Co-energies T’, U’. 


In non-linear resistive circuits G and J possess certain properties. 
It may be shown that similar analogous, yet physically quite distinct 
notions arise in connection with inductor and capacitor elements ; with 
reactance elements time is involved. Thus with reference to fig. 1 (a) 
the rectangle ®;,7,,, is divided into two parts T and T’ by the characteristic 
Q is the instantaneous working point. 


- 
ie | i d®, the inductive energy (area above curve) (1) 
i} t Mey ae Swe 
r +}. ® di, called the inductive co-energy* (area below CUrVe),-. ee” (ae 


; '* . . awe . . - 
; an invented term ; It is, in a sense, dual to energy, and is a quantity of 
sufficient fundamental interest to justify the allocation of a special term 
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where O=Ng¢, the flux-linkage. Note that T, in a non-linear inductor, 
may be derived from the work expression vi dt : 

ay. t d® ae 

| T=] vidt=l it di=|id®. 2. 2 1, (3) 
Similarly with non-linear capacitor elements (fig. 1 (6)) : 


We | q dv, called the capacitive co-energy (area above curve), 


(4) 


rq 
Ua | v dq, the capacitive energy (area below curve). 
~ 0 


Again U is derivable from the work done : 

t ; a At dy 

: u=|_ Vv ae va 

Observe that in linear cases T=T’ and U=U’ and energy=co-energy. 

It will be shown later that this equality obscures the importance of the 

co-energy. 

The characteristics in fig. 1 have, for simplicity, been taken as simple 

(unique, time invariant and single-valued ; hysteresis is excluded). 


at=|" v dg. eee Se 59 


§ 4. THEOREMS CONCERNING Muuti-MisH Non-Linear REACTIVE 
CIRCUITS. 
4.1. The Energy Rectangle Theorem (fig. 2). 

A linear or a non-linear circuit of resistor elements may be represented 
by a set of contiguous rectangles fitting inside a bounding rectangle ; each 
component rectangle has sides of length v,, and 7,,, representing the 
instantaneous voltage and current on any one branch of the circuit, 
while the area (v,,2,;,) represents the instantaneous power in the branch. 
The bounding rectangle represents vp, i), the source voltage and current 
(and power). 

This is a valuable concept* which, although first put forward for linear 
resistive circuits many years ago, is not well known (Hering 1927). It is 
valid also for all-inductor and all-capacitor circuits, linear or non-linear, 
though it applies only instantaneously. Fig. 2 shows an example; by 
whatever path is chosen between the external terminals AB, the total 
P.D. is vp, while at the same time, if a dividing line XY be imagined to 
cut the circuit anywhere between the terminals AB, it always intersects 
the same total current 7,. That the set of rectangles v,,i,, always forms 
contiguous areas inside the reactangle v»i, merely illustrates the truth of 
(a) Kirchhoff’s equations and (b) the conservation of energy. « In those 
cases in which an internal branch contains a generator the component 
rectangle may have a side —v,, or —?,,, and the construction may require 
overlapping rectangles. 

* Brooks, Smith, Stone, and Tutte: “The Dissection of Rectangles into 
Squares’, Duke Math. Journ., Vol. 7, Dec. 1940. For explanation see Cherry 
(1949). 

4K2 
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If the diagram in fig. 2 (b) is imagined as being turned on its side, ig! ya 
everywhere v and 7 are interchanged, the diagram is relevant to the dua 
circuit, dual in the sense of Cauer’s construction (Cauer 1934), everywhere 
voltage and current are interchanged, as are also meshes and nodes. 

4.2, The Superposition of Energies and Co-energies. 
4,2,1,—Theorem, At any instant the power supplied by various sources to 
any linear or non-linear circuit equals the rate of change of the 
energy stored in the inductor and capacitor elements plus the power 
dissipated in the resistors. 

This follows of course from the Conservation of Energy, and is illustrated 

by the rectangle theorem referred to above. We may write 


P=F+ 5 (0-+U), a VS a 


Fig. 2. 


ee 
CURRENT 


(2) A circuit configuration, and (6) its energy rectangle representation. 


where P is the power supplied and F is the total dissipation. It will now 
be shown that the co-energies T’ and U’ possess certain properties when the 


circuits consist of any network of like-elements, i.e. all-inductor or all- 
capacitor. 


4.2,2.—Theorem, At any instant the co-energy in any “‘ sumple”’ 2-terminal 

non-linear all-inductor circuit regarded as a single element between 

the terminals is equal to the linear sum of the co-energies in the 
constituent inductors, 

For the moment we exclude hysteresis and mutual inductance and also 

assume that the cores of the inductors are immovable. A simple circuit of 
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this type is illustrated by fig. 3(a). Let vv,...... v, be the potentials of 
the various nodes of the circuit, assuming, without loss of generality, that 
U,=0. Let 7;, be the current in the branch jk and i, be the current fed to: 
the external terminals. Consider the growth of the branch currents from 
zero at i—0. Due to an increment 4i,, the co-energy in the branch jk 
will change by an increment AT‘,; the total co-energy in the whole- 
circuit will then change by AT’, where, to a first approximation, 


A Meee ce eg ey. | Dr t2 |) et 4 (2) 
jk jk 
(see fig. 1(a)), the 4 appearing by virtue of the double-suffix notation. 


But t 
®j=| (vj—v4) dt, . dies “el ere (Sy 


PITTI LLL 


(a) Inductor circuit. (6) Equivalent inductor. (c) Rectangle 
representation. 


A non-linear 2-terminal inductor circuit and the “equivalent ’’ inductor. 


which is independent of the voltage wave form, and dependent only upon 
its integral. The integration constant is zero if we assume O®=0 at t=—0. 
Then 


rt 


: t 
av=125[f'@,—a)d] din= 2 | vdt| 2 Ain], . . 0) 
jk j 70 k 


since i;,—=—?,; But from Kirchhoff’s current law 2 41,;,=0 at every 


0 


k 
node except at the external terminals, where it is At). At this point 
v;=U. Therefore 


t 
AT’ =Ai,. i v, dt. 
(0) 


This we may write as 
AT’ = Ain . Dp, ae het Wk ee eee ee | (10) 
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where ©, is a hypothetical flux associated with the integral of the applied 
voltage v). Integrating these increments of co-energy, 


Ta ("@, din keen 


~ 0 
which itself represents a co-energy expressed entirely as a furiction of the 
input quantities Dp, 2. 

This theorem suggests that it may be possible to construct the 
characteristic of the “ equivalent’? inductor point-by-point. Before so 
doing it is necessary to establish a law governing the distribution of 
currents amongst branches, in order to define the instantaneous working 
points v;,, 7;, (see sub-section 4.3, below). 

. Since both the energies T,, and co-energies T),, add linearly in this way 
to give the overall energy T and co-energy T’ of the inductor circuit, 
as measured at the external terminals, then so also do their sums add 
linearly : 
tet ies ca [Ty +752] 
jk 


‘or Al 
Doin =X 2 [Pjx Uc]. | 
j ik 


This conservation of flux-current is of considerable importance. It 
suggests that a rectangle construction, similar to that in Theorem (4.1) 
(fig. 2) is applicable to an all-inductor circuit. Thus fig. 3 (c) illustrates 
the rectangle diagram appropriate to the inductor circuit (@) at some given 
instant ; by whatever path we ‘trace between the external terminals we 
shall encounter the same total flux ©, since 


Up=Uj tet V3--+-- 
and integrating gives = @,+ 0,10,..... 

Note that each rectangle may be divided into two areas T,,, and T’,. by 
the element characteristic, corresponding to the diagram fig. 1 (a) ; one 
such rectangle division is illustrated in fig. 3(c). From Theorems 4.2.1 
and 4.2.2 these areas superimpose separately, giving the corresponding 
areas for the composite characteristic. 

Duality is now interpreted in quite a new way; turning the diagram, 
fig. 3(c), on its side interchanges flux-linkages ® and currents i everywhere, 
and converts nodes to meshes. This provides completely general rules for 
finding the “equivalent circuit ” of any loaded transformer ; alternatively 
it shows that a circuit containing only electrical constraints, as that in 
fig. 3 (a), may be assembled upon a transformer and so be given only 
magnetic constraints. This has been made the subject of a previous 
study (Cherry 1949): it is here shown to be an equally useful theory when 
the magnetic circuits have non-linear characteristics. : 


weed These theoretical points concerning all-inductor circuits may be 
a ae A have their equivalents in relation to all-capacitor circuits. 
1g. 1 (0) illustrates a typical “ simple ” non-linear characteristic, with P 


an instantaneous working-point and U’ = 
at that instant. A Me eee aa ae 
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4.2.4.—Theorem. At any instant the co-energy in any “ simple” 2-terminal 
non-linear all- -capacitor circuit, regarded as a single element 
between the terminals, is equal to the linear sum of the co-energies in 
the constituent capacitors. 


This theorem is exactly dual to Theorem 4.2.2. For proof, it is 
‘simplest to work on a mesh basis, using ij i,..... ?, a8 mesh currents and 
to consider the co-energy increment 4U/, due to a branch voltage incre- 
ment 4v,,. Again, flux ® in the above theorem is replaced by charge q 
(the dual). The proofs are then analogous. 

( Again, since both energies U,,-and co-energies U/, superimpose linearly, 
so also must their sum, 


U+U'=4 FE (U,,+Up] 


rs Vrs» 


or Ud=t 2 v 
and a rectangle theorem similar to that illustrated. by fig. 3 (c) is applicable, 


; anes Stationary Value Theorems for Non-linear Systems. 


The following theorems establish variational principles which govern the 
distribution of currents in the branches of non-linear “simple” all-inductor 
or all-capacitor circuits. These theorems are later extended (§§5 and 6) 
to circuits including mutual inductance and rotating (or sliding) cores, 
such as those of rotating machines. 


4.3.1—Theorem. Ina circuit of voltage-driven non-linear inductor elements 
the total energy T has a stationary value* with respect to voltage 
variations. 

Working on a mesh basis, let i) 7,..... i, be the Maxwell mesh 
circulating currents. These currents will take up their natural values ; 
however, imagine the voltages across elements v,, to be artificially varied 
by 4v,, to unnatural values such that around every mesh Kirchhoff’s 
voltage law still applies : 


2eAte == 0 Pe oasis pe) tas 2s 74 (13) 


(where these variations 4v,, are zero across all voltage-generator terminals) 
while the current law is violated, at the various nodes 


De One nmerard witty OM (1A) 


(the circulating currents no longer have “ natural” values). The energy 

stored in any element rs is thus varied by (7,—i,)4®,, to a first approxi- 

mation ; the total energy of the circuit is then varied by AT: = 
ZA At ON Ee WCC | GARR oe a EEN 


ne ey ee ee Se 
* }, e, maximum, minimum, or inflected. 
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But 4,,=| Av,, dt, assuming that ©,,=0 at t=0. Then 
0 


t 
AT=} DE (i,—i,) | dv, dt 
r 8s 0 


t 
i [=/ Av,, dt | , , oS ee 

r 8 40 
since v,,=—v,,. But the integrands here sum to zero, from equation (13). 
Hence AT=0, 95 5 sc 


or the energy T has a stationary value, at every instant. 


4.3.2—Theorem. Ina circuit of current-driven non-linear inductor elements 
the total co-energy T’ has a stationary value with respect to mesh- 
current variations. 


Working on a node basis, let v v,..... v, be the node potentials. 
Imagine the current i, in any branch to be artificially varied by 47,, to 
some unnatural value such that at every node Kirchhoff’s current law 
holds, while the voltage law is necessarily violated : 


DAi,=0; | Des ones one 
k 


The co-energy increment in the branch jk is 4T%;,,; summing in all the 
branches, 


AT’=}32E [ (v,—v,) dt. At...) ee 
which may be written cts 
4v= fo dt E diy |, 
since 7;,——1,;. But from equation (18) this is zero, or 


AT’=0 at every instant, a: Lottel: © irom tot aur 


meaning the co-energy 7’ has always a stationary value. 


4.3.3.—Theorem. Ina circuit of non-linear capacitor elements driven by any 
number of voltage (current) sources the currents will distribute so as 
to give the total co-energy U' (energy U) a stationary value at any 
imstant with respect to voltage (current) variations. 


The proofs of this theorem are exactly dual to the proofs given above ; 
for node we substitute mesh, voltage for current, charge for flux. 
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The properties of all-inductor or all- -capacitor circuits may be 
summarized as in Table I. : 


TABLE I. 


All-inductor All-capacitor Mixed circuit 


Voltage drive AT=0 AU’=0 4 | é (T—U") dt=0 


(voltages varied) | (voltages varied) 


Current drive AT’=0 AU=0 A | (TU) di=0 
(currents varied) | (currents varied) y 


The properties shown here for the “ mixed circuit ’’ will be demonstrated 
later in sub-section 7.3. 


Fig. 4. 


a 


(a) Transformer with non-linear core. 


(6) Transformer characteristics. 
Energy stored in a non-linear transformer. 


§ 5. IncLUSION oF MutuaL INDUCTANCE. 

If circuits contain mutual inductances they nevertheless satisfy 
Kirchhoff’s laws, and the instantaneous power acting on any branch jk 
i8 V4, 4 Where, however, v;, may be a function of more than one current. 
Consequently the Theorem 4.2.2 relating to the superposition of co- 
energies applies equally well to circuits containing transformers, with 
linear or non-linear iron cores. 

However, the quantity ‘“‘co-energy ”’ needs careful definition. Consider 
any two branches of a circuit which are coupled by mutual inductance 
possessing non-linear properties (fig. 4 (a)); the energy stored at any 
instant is 


t t 
Ty=| Vih1 dt+ | Volo dt, ° e ° ° . ° (21) 
0 0 
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where, if ®, and ®, are the total fluxes linking the primary and secondary 
coils respectively, . 
Oh = di P(t, ty) and SAE Oy(74, V2), ers) 
Assume that families of characteristics have been plotted, as in fig. 4 (6), 
giving relations between ®, and 7, for different settings of i, and also 
between ®, and i, for different settings of 7, Assume also that the ® 
values are unique functions of 7, and i, and that Ty depends only on 
the final values of i, and i, and not upon the manner in which they 
reach these values; that is, hysteresis is excluded. 
Substitution from (22) into (21) gives T,,; four components : 
i. OO, .. ph, @,, ue (ao OO, oe 
Ty= [. by Wi, - dit | (h Fi, digt joka Wi, dig+ |, a 3h di,, (23) 
If the point P, Q is the working-point (currents I, and I,) in fig. 4 (6), 
it may be reached in two particularly simple ways : (a) With 7,=0 raise 
i, to I, (points A and E on the characteristics); the stored energy is 
represented by the area OAB corresponding to the first integral in (23). 
Next raise 7, to I,, keeping 7,=I, constant (points P, Q); the energy has 
two components, represented by ABCP and EQF, corresponding to the 
2nd and 8rd integrals above; the 4th integral is zero. (6) Alternatively, with 
1,=0 raise 7, to I, (points C, G); then raise 7, to I, (points P, Q). The 
energies are respectively given by the areas OCD, DCQF and GPC, 
Equating these stored energies, 

OAB+ ABCP+ EQF=OCD-+ DCQF+ GPC, or OAPG=OCQE, (24) 
so that the shaded areas in fig. 4 must be equal. Similarly the areas 
enclosed between any corresponding pairs of characteristics, at every 
instantaneous pair of working currents, must be equal. 

Consider sets of characteristics to be plotted, as in fig. 5, at small 
current intervals di, and 4i,. The actual loci of the currents 7, and 7, 
starting from zero may be plotted on these characteristics, as illustrated by 
the dotted lines. These loci may then be approximated by repeating the 
process, described above, of increasing i, and 7, alternately, but now in 
small steps 47,, diy. If this is done the sum of the two areas T, and T, 
lying above the working loci is the total stored energy in the transformer 
at the working point P,Q. Then we may define the total co-energy as the 
sum of the areas T) and T) lying wnder the working loci. 

To see how the various theorems, described previously, may be applied 
to circuits containing mutual inductance, consider one current i, to be 
increased by di, The total increase in co-energy is AT’ =®©, . Ai, since 
AT, will be zero. Similarly, if i, is increased by Az,, then the total 
co-energy increases by AT), while 47.—0. Thus Theorems 4.2.2, 
4.3.2, relating to superposition and stationary values of co-energy, are 
still directly applicable since they depend only on Kirchhoff’s current law 


and upon expressing the total circuit co-energy increment as AT” (see 
equation (7)). 
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Similarly, with reference to fig. 5, if the flux ®, is given an increment 
A®,, the total energy changes only by i, 4®,; in a similar manner, then, 
Theorems 4.2.1, 4.3.1, relating to superposition and stationary values 
of energy, are also applicable since they depend only on Kirchhoff’s 
voltage law and upon expressing the total circuit energy as AT (see 
equation (15)). 


5.1. Equivalent Inductor Theorem. 


_ From Theorems 4.2.1, 4.2.2 and 4.3.1 or 4.3.2, above, it is clear that 
the following theorem is also valid: 


99 


5.1.1.—Theorem. For any circuit of non-linear “ simple’ inductors, with 
| two external terminals A, B, the concept of a unique “ equivalent ’”’ 
“inductor, connected between A, B, is a valid one. 

The total energy T and total co-energy T’ at any instant equal the sums 
of the instanteous energies and co-energies, respectively, of the constituent 


inductors, while T (or T’) has a stationary value with respect to voltage 


Energy and co-energy in a non-linear transformer. 


(or current) variations, if a constant voltage (or current) generator be 
connected to the terminals A, B. The theorem is valid also when mutual 
inductors are present, which may have non-linear characteristics, of the 
“< simple ”’ (non-hysteresis) class. 

This theorem suggests that point-by-point methods may be evolved, at 
least in principle, whereby the equivalent-inductor flux/current char- 
acteristics may be computed. Finally, an analogous theorem is applicable 
to any circuit of “ simple ” non-linear capacitors. 
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§6. Systems wiTH MovABLE Crrcurts—RoTatTing ELECTRICAL 
MACHINES. 


Consider any one coil of a rotating (or a linear stroke) machine, the 
reluctance of the magnetic path depending upon @, the position of some 
moving magnetic member (fig. 6). Then ®=®(i, 6). Let the coil 
current be held constant at i=i, while the rotor is moved 40; if F is the 
rotor force, F=F (i, 0), and the mechanical work done=F.40. The 
working points move from P to S, and the electro-magnetic work done on 
the circuit connected to the coil terminals is proportional to the area 
PQRS. Further, the stored inductive energy falls by an amount pro- 
portional to the difference of the areas OPQ and OSR. Thus 


F . AA=PQRS—(OPQ—OSR)=POS, . . . (25) 


being the area JT; between the adjacent characteristics. Note that in 
the linear case POS=4 (PQRS). 


Fig. 6. 


The characterstics of a coil on a rotating machine. 


Again, if we regard the co-energy of the coil as the area lying under the 
working characteristic, T’, this mechanical work F . 4@ (the area POS) 
represents a loss in co-energy which may be called AT;. However, if 6 
be held constant and 7 increased by 4i, this co-energy increases by 
AT; =®. Ai. Thus the total co-energy increment in a short interval of 
time, of the coil jk for example, is AT’, : 


The total co-energy of this coil is then found by integrating from zero: 


j bik : na os 
Th=| Degen — Bx (i, 8) « dB. ee CoG 


0 
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If the machine has a number of such coils, coupled together in some 
complex circuits, with two external terminals, the total co-energy 
increment of the entire machine is then : 


jk 
t 
But O,== | (v;—v,) dt. Using an argument similar to that used in the 
0 


proof of Theorem 4.2.2, we may write 


ct 

AV=Z |v, dt E Ai, | Lesa Ee 
7 ~0 k j k 

and again, at every node 2'4i,,—=0 except at the external terminals, where 

it equals 47), and where v;,=v). Then 


AT'=GyAiy—L EEE. AB . . . . . (29) 
ne 


t 
where ®,—| v, dt: this equals the increment in electrical co-energy D, Ai, 
0 


at the external terminals less the total mechanical work done. This is 
then a modified version of the co-energy theorem 4.2.2 for rotating 
machines ; note that in simple machines, with a single rotating core, all 
the 40,, are identical. Again, mutual inductance is usually present but, 
as deduced in the last section, makes no difference to this Theorem. 
However, note that the concept of a unique ®), 7 characteristic, as 
measured at the external terminals, in not now valid since ®, is a function 
of both i) and 0, which are independent variables. Nevertheless, if all 
6 values are equal, as in most rotating machines, (29) expresses AT’ 
entirely in terms of external quantities ©), %, 0, and so suggests that 
a family of equivalent characteristics is a valid concept. 

The stationary value Theorems, 4.3.1 and 4.3.2, may also be inter- 
preted for this case of a rotating machine, or other electro-magnetic- 
mechanical system, because the proofs of these theorems were not in any 
way restricted to time-invariant inductor elements. 


§ 7, THe GENERAL Non-Linear Circuir oF InpucToRS, CAPACITOR, 
RESISTORS AND GENERATORS. 


It appears difficult to find general theorems relating to circuits 
containing a mixture of elements. However, certain properties emerge 
from expressing the circuit equations in Lagrangian form. 


7.1. Non-Linear Circuit Equations in Lagrangian Form. 
Fig. 7 (a) represents part of some circuit of non-linear R, L, C elements 
and voltage generators. Let the coordinates 7, 7, ....%, be the Maxwell 
mesh currents. Then around, say, the rth mesh 


wr top tRUr=er ° . . . . . . . (30) 
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where ,”, cv, Rv, € are the total inductive, capacitive, resistive P.D.’s and 
generated voltages around a mesh. Consider these voltages separately ; 
firstly ,v,. Write q,=1,: 

The Ge in any inductor, rs, in the 7th mesh is a function only of the 
total current (i,—i,), where s denotes the adjacent mesh. The co-energy 
T’, is then 


i,—i,) 
| 6d 2 ee 
0 
(see fig. 7b). The total inductor co-energy in the entire circuit is then T’ : 
é_—é,) 
T=422|" "Odi, = ne 
fs O 
Fig. 7. 


= (6) 


A general non-linear circuit ; L, R,C elements and voltage generators. 


Differentiating the total circuit co-energy T’ with respect to 2, (or 2,) 
from (32) gives 


oo er ae 


since only those inductors around mesh r carry the current i,- In this 
summation s=0 relates to the external mesh ; the current i, is assumed 


zero, without loss of generality. From (33) the total mesh r inductive 
voltage is 


as asfols 
Dep dt of OD = dt () oT ite als age ‘ (34) 


Similarly for the capacitor voltage around the rth mesh ; the capacitor in 
the branch rs supports a charge (7,—,) and the voltage across this element 
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iS V,,=0V(7,—q,). The total energy (not co-energy) of all the circuit 
capacitors is then U (see fig. 7 (c)) : 


%,— Qs) 
Uap zz) Veeg weiee e (38) 
r § 0 
Then Cay x 
oq, = hae alae Of ge (Por AEE or cn ts (36) 


the total capacitive voltage round mesh r. 

Finally the resistor voltages. These may be expressed in terms of the 
resistor content and co-content (see §3). The resistor content of the entire 
circuit is G: 


(i 
eR og (87) 
Tye Sind) 
so that, differentiating, 
CGA 
Oe 23; OP aaa TR Ues . : . . . . . (38) 


$= 
Then substituting-in (30) from (34) (36) and (38), the equation of the 7th 
mesh of the circuit is 
d. (ot! Ou 0G 
which is recognized to be Lagrange’s equation of a dissipative system. If 
we define the Lagrangian function of the circuit as Z=T’(q)—U(q), the 
equation (39) may be re-written 


d (of OF. 0G 
a (aj,) ~ 39, 3 oo 
and a similar equation applies to every mesh r=1,2..... mw. If required, 


this form of the equation serves also for circuits containing movable- 
core inductors, as with rotating machines. Thus, as seen from equation 
(27), the co-energy T’ then. includes the mechanical energy and becomes 
a function of both z and 6. 


7.2. The Dual of a Non-Linear Circuit. 


The concept of a pair of dual or inverse circuits (A. Russell 1904) is, 
in the linear field, a most important one. In this present paper it 
has been noted that certain theorems may be derived from others 
directly, merely by the substitution of capacitors for inductors, of voltages 
for currents, of nodes for meshes, and vice-versa. This may now be 
generalized formally, as applying to non-linear systems. 

The Lagrangian equations (40), above, were derived from the mesh 
equations (30) of a given circuit. Imagine a topologically dual circuit 
pattern which, if drawn superposed on the original (‘‘ parent ”’) circuit has 
one of its ~ nodes lying inside one of the ~ meshes of the parent circuit 
(Cauer 1934). Then the node equations of the circuit will be 


wytotjs tet . . . . . : . : (41) 
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The capacitive co-energy U;,, the inductive energy ve and the resistive 
co-content J, in the elements connected to the node 7 are 


r(vs— vp) 7(Dj— Pr) (ese), 
Un= | me Yin» WV 5 ibaa | vin IP; Ing I, noize > QW, 
0 
if v, V,..-V;.-+++V, are the nodes potentials and if v,=@,. By differ- 
entiation 


ds/(0U.\ ae CATS Wary ment FOV. 
i art) =e 0D; i Nicos eae 


But note that, for example, 0U;,/0v;=0U'/dv;. Thus the suffixes may be 
dropped and the expressions U’, T, J adopted, which refer to the entire 
circuit. Adding these current components, as in (41), 
d (=) od ee iS 


Cee 42 
di\ dv, ‘ (42) 


OD ae 00% 
which again is a form of Lagrange’s equation (40); however, in this 
current-driven case the Lagrangian function # is 
FG=U' (ve) —T(D),. 0 te ee 0 
A comparison of equations (42) with the mesh equations (39) establishes 
the duality concept. Thus, using the symbol —, meaning “ replaced by”, 
we have that, given the set of equations (39) referring to the meshes of a | 
parent circuit, 
. dT’ /0q,>0U'/dv; 
675 dU/aq,>dT/a®, ee ony ye? 
aS, aG/aq,>dJ/dv; | 

leads directly to the node equations of the dual circuit. 

In practice this may be of limited interest since it may not be possible 
to construct physical elements having the required curvature of 
characteristics. Nevertheless, the transformations (44) define the 
theoretical relation required between the element characteristics. For 
example, if the three characteristics (b) (c) (d) in fig. 7 refer to elements of 
the parent circuit, then exchange of the ordinates and abscisse gives the 
characteristics of the corresponding elements in the dual circuit. 


7.3. Dissipationless Circuits—F orce-free Oscillations. 


If a dynamical system can be expressed in the form of Lagrange’s 
equation (40) it is a matter of simple algebra to show that, if the dissipation 
is zero (G=0) and if the generators are removed (e,=0), the force-free 
oscillations which may occur in what is now a non-linear but purely 
reactive circuit are governed by Hamilton’s Principle, 


t 
A| ¢.d=0, oa eee EO 
0 


a variational principle (Whittaker 1904), meaning that the natural wave- 
forms of the currents are such as to keep the integral of # (=T’—U or 
U’—T), in the dual, stationary at every instant. (Note that in the linear 
case T’=T and the more familiar form of Y= (T—U) appears.) 
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This variational principle has been added to the list in Table I. (sub-section 
4.3.3). This principle tells us nothing new about the non-linear circuit ; 
it is to be emphasized that it is merely the most fundamental statement of 
its behaviour, being in the form of a single equation. Ina given example, 
if the equations of the non-linear characteristics of the various elements are 
known, it may be possible to express the Lagrangian as a function thus: 


Pte ee where 11 (9,9 ;-t); 
jk 

U=322U;,, where U;,=U(7,9¢ 4), 
jk 


and substitution into ¥ and calculation of the result of arbitrary variation 
of the q,q variables leads to a set of m simultaneous equations of the 
Lagrangian form (40) (with G=0) which are the non-linear equations of 
the n-mesh circuit. 

In the linear field Hamilton’s principle (45) leads to a general formu- 
lation of the conditions for the natural modes of oscillations, as expressed 
by Foster’s Reactance Theorem. It is possible that in the non-linear 
field there may be some corresponding set of energy relations which 
defines the various characteristic types of stability or conditional stability 
that can occur ; this possibility will not be pursued here. 


REFERENCES. 


CuErRRY, HK. C., 1949, Proc. Phys. Soc. B, 62, 101. 

GuitLemin, E., 1933, Communication Networks, Vol. II. (John Wiley). 

JEANS, JAMES, 1948, Electricity and Magnetism, 5th Ed. (Cambridge : 
University Press). 

Sraruina, 8. G., 1912, Electricity and Magnztism, Ist Ed. (Bell). 

Wuittaker, E. T., 1904, Analytical Dynamics (Chap. XX.) (Cambridge : 
University Press). 


SER. 7, VOL. 42, NO. 333.—OCT. I95I1 4L 


Pewee.) 


CXVIII. A Perturbation Treatment of Closed States i Quantized. 
) Field Theorves. 


By B. TouscHEE, 
Department of Natural Philosophy, The University Glasgow*. 


[Received May 30, 1951.] 


ABSTRACT. 
The paper discusses an iteration method by means of which it Is possible: 
to obtain approximate solutions for field theoretical problems involving 
bound states. 


iB, 

Tue treatment of bound states in quantized field theories has recently 
received much attention. After the elimination of divergences by the 
work of Tomonaga, Schwinger, Feynman, Dyson (1949) and others this 
problem which is related to the convergence of the whole series of 
approximations progressing in powers of the coupling constant, has 
become the central problem of field theory. Dyson (1951) attributes the 
difficulty of the problem to two facts, (i.) the free particle nature of the 
initial and final states in the customary S-matrix formalism and (ii.) the 
impossibility of representing the wavefunction of a bound state in terms 
of a convergent power series in the coupling constant. A failure of the 
formalism which apart from the renormalization procedure is essentially 
Born’s approximation method may also be expected in those cases in 
which bound states manifest themselves in their virtual form as in the 
theory of nuclear resonances. 

Recently Ferretti (1951) and Dyson (1951) have attempted to overcome: 
the difficulties associated with the occurrence of bound states; Ferretti 
by considering adiabatic changes of the coupling constant, and Dyson 
by considering changes of the coupling constant in a fictitious world 
coincident with the real world at time ¢. At t= oo Ferretti thus enforces. 
the justification for using free particle expansions and Dyson succeeds .in 
doing the same at least in the fictitious world. In this note an alternative 
method is suggested which, though its convergence has yet to be demon- 
strated and its relativistic properties have to be investigated, allows one 
to treat stationary states and is capable of leading to quite definite 
predictions. Considerations to a similar end have been recently investi- 
gated by Power (1951), who could showy that after minor modifications 
the two methods of approach are essentially equivalent. The method 
used here consists in adding a suitable operator representing an interaction 


energy to the Hamiltonian of the free particles (nucleons and mesons) 
ra ee ee 
* Communicated by the Author. 
{ Private communication, 
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and subtracting it from the customary interaction energy, which represents 
the possibility of creating or annihilating Bose particles by a corresponding 
transition of the Fermi field. The interaction operator thus introduced is 
subject to the additional requirement that it should leave the number of 
Bose particles unaltered. It may be such that bound states occur between 
nucleons or mesons, or even between nucleons and mesons. Orthodox 
perturbation theory is then applied to the new perturbation energy, 
v7. e. to the creation-annihilation term minus the supplementary interaction 
term. Finally one tries to choose the supplementary interaction in such 
a way that perturbation theory does not give any contribution to that 
part of the total interaction energy which leaves the number of Bose 
particles unaltered. The latter requirement can only be satisfied by a 
method of successive approximations similar to the iteration method leading 
to the determination of the Hartree field in atomic theory. The adjustment 
of the auxiliary interaction can be carried out in steps proceeding with the 
even orders of perturbation theory. The point of the method is that 
wavefunctions of the right asymptotic behaviour—7. e. that of bound 
states—are forced on the formalism before the perturbation treatment 
begins. 
§ 2. 

In this section the formalism will be outlined for the first stage of the 
process, which leads to an adaptation of the wavefunction to the second 
order of the customary perturbation formalism. Let 


eH Pe ee tee se A 1) 


-be the Hamiltonian of the system. HH, is the Hamiltonian of the free 
particles : H,=H,,+H, where H,, is the Hamiltonian of the free nucleon 
field and H, the Hamiltonian of the mesons. From (1) one obtains the 
Schrédinger-equation 

(in natural units, A=c=1). We now introduce the additional interaction 
term V—which may absorb some of the renormalization terms—and write 


iy=(H,+V)y+(H'—V)xy. 2 . «© 2 ee (8) 
Taking H,+V as the zero order approximation whose Schrédinger- 
equation 


Wee (a eve me meme cece 6. xo (A) 
can be solved by orthodox methods, we put 
aOR ee Re rec s\ of ip teh gD) 


Inserting this into equation (2) one obtains the wave-equation of a modified. 
interaction representation : 
WeeKSuwithe Keed hove aie. 5 (8) 
Equation (6) differs from the corresponding equation (V=0) occurring in 
the usual formulation of the interaction representation by the properties of 
4L2 
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the states to which the indices of the S-matrix refer. To see this more 
clearly it is useful to fix the formal meaning of the quantities occurring in 
equations (1) to (6). The Schrodinger-vector k= x(q, #) is @ function of the 
physical variables q of the system and the time ¢. The solutions of 
equation (4) are of the form ¢y(q, t)= x(q) exp (—tExt), where the Ey are 
the energy values of the free particles in the interaction field V. There may 
bea discrete spectrum superimposed on the continuous spectrum of the centre 
of mass motion. (5) reads explicitly Saale IC exp (—?Eyt) Sy(t). Ifthe 


states represented by the solutions of (2) are labelled by an index «— 
distinguishing the true states of the system—y(q, t) will be of the form 
W (7) exp (—iE,£) and S turns into a matrix S\(¢) which thus represents 
the coefficients in an expansion of the true solutions in terms of the solutions 
of the model system (4). (In the customary interaction representation the 
true states and the free particle states are usually labelled by the same 
index so that the mixed nature of the S matrix is masked.) Our equation (6) 
therefore differs from the corresponding equation of the interaction 
representation in that in (6) the ‘‘ rows ” of S refer to the (possibly closed) 
states N of the model system and therefore do not differ qualitatively, 
z.e. in their asymptotic behaviour from the true states as in the 
interaction representation. 
The typical matrix-elements of K are of the form 


Kory exp [t(Ey—Ey)t]}, . . . . . 2 2 {¥) 
If it is assumed that V can be chosen ocg? one may write 
K=K%+K® with K®=¢-1H’¢, K®%=—d¢-1Vd, . . (8) 


where the upper indices indicate the lowest power of g, 7. e. the way in 
which an operator vanishes if g>+0. Thus in general F™=g"F(g) with 
F(0)A0, because of the implicit dependence of 4 on g. The successive 
approximations of S become : 


S@—T, 
nt 
Sl) ee | K(¢,) dt, 


rt t “th 
Si] KE) de, +(—H# | a KG) [aK . 0) 


The first equation establishes a 1-1 correspondence N<+« between the true 
states of the system and the states of the model system. 

We now make use of our freedom to choose V. Clearly the ideal choice 
of V should ensure that the perturbation method when applied to 
equation (6) yields most rapidly converging results. However this aim 
can hardly be formulated in a mathematical fashion. For definiteness 
we therefore try to determine V in such a way that the pure interaction 
terms arising in the second approximation in equation (9) vanish. For 
this we remember that H’ will only contain terms corresponding to the 
emission and absorption of a meson. The same will hold for K® since V 
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leaves the number of mesons unaltered. It will therefore be possible to: 
split the second term on the right-hand side of the last equation (9) into 
two terms 


t “ty rt 
| aK, | dtxK(t,)= | dt,I+I), . . . (10) 


where I represents all those matrix elements for which the number of 
mesons does not change and II. those matrix elements for which the 
number of mesons changes by 2. If therefore 


Kee ee ae her fie 68 cata. 


all the interaction terms will vanish in the second approximation and the 
interaction potential is “adapted to the second order of perturbation 
theory ”’. 

Equation (11) is an implicit relation from which V has to be determined. 
It appears that this can be done in rapidly converging steps, for the 
interaction energy I draws most of its strength from the highly excited 
states of the system—owing to the short range of nuclear forees—so that 
the deformation of the particle wavefunctions caused by the presence of 
the interaction energy I should be small. The self-consistent ‘‘ potential ”’ 
V can thus be determined in the following manner. Let 7, .... denote 
free particle states. As a first step in the process leading to the self- 
consistent V satisfying (11) one chooses V/,,,=iI,,,.. Inserting this V 
into equation (4) one obtains a new set of states n’, .... with energies 
E,,, .... lis now determined in terms of the new eigenstates n’. A new 
V” is determined from (11), viz.: Vjj..=tI qn and the process is 
repeated. If this procedure converges it will lead to a limiting set of 
states N=lim (n’, ....; ”, ...., .-..) and to a limiting potential 
V=lim (V’, V”, ....), which when inserted into equation (8) will satisfy 
the condition (11). A rough estimate of the convergence of the procedure 
shows that, for example, in the case of the nucleon—nucleon interaction, 
successive steps of the approximation will add a correction of the order 
E,/» where E, is the binding energy of the “ deuteron” and y the rest 
energy of the meson. This rapid convergence stems from the fact that the 
energy denominators of the customary second approximation formalism 
are of the order » and is in this way related to the short range 1/. of 
nuclear forces. This seems to indicate that the present method may not 
be applicable to quantum electrodynamics where the range of forces is oo. 

In second approximation V will—in addition to self energy-terms— 
contain only two terms corresponding to the interaction between two 
nucleons and between a nucleon and a meson respectively. These 
interactions can be represented by the two 2-vertex Feynman graphs 


XN 
N 
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The internucleon term represents the two-body contribution to nuclear 
forces. The second graph will cause a distortion of the mesonic waves | 
in the neighbourhood of a nucleon. It represents the perturbation energy 
in terms of which the scattering of mesons by nucleons may be described. 
In the case of scalar charged mesons it will correspond to an attractive 
force between protons and 7+-mesons and neutrons and a~-mesons 
respectively. The forces between protons and 7~-mesons and neutrons 
and 7+-mesons are of opposite character. It should be noted that these 
forces cannot be described in terms of a simple potential function V(r), 
in which r is the distance between the meson and the nucleon. It will 
rather be of the character V(r, r’) attributing a certain energy to the 
annihilation of a meson at r’ and its re-emergence atr. In k-space the 
first ‘“‘ Hartree-approximation ”’ for infinite mass nucleons will be 
2 

V' (deh) = (+h), ss (18) 
where L3 is the volume of the normalization cube. The states containing 
one positive meson may be expressed in terms of the one meson Hilbert- 
vector : 

pineal ge A Ve 8. SO, pee ck 


Here y, is the Hilbert-vector of the vacuum, a* the operator 
representing the creation of a positive meson of momentum k and u(k), the 
one meson Schrédinger function normalized to Y|u(k)P=1. In our 

k 


approximation energies and wavefunctions of the meson in the presence 
of an co-mass nucleon are the eigenvalues E and solutions u(k) of the time 
independent Schrédinger-equation : 

2 Bxcelcog— EB) V Gee UE) 0.9 ee 
For large energies V(k, k’) will represent a small perturbation and Born’s 
approximation will give approximate expressions for the distortion of the 
wavefunctions near the nucleon. The distortions may become considerable 
for low meson energies. The possibility of bound states can be 
investigated by using Hyleraas’ variational method. It turns out that 
at least for g?>7* the system becomes unstable, giving bound states of 
infinite negative energy. This can be shown by using the trial function 
u(k)=Ne-“/k (where N is a normalization factor, and « is the parameter 
which is to be varied). The energy E comes out smaller than —const/«. For 
values of g less than a certain critical value—of which 7 is a rough 
estimate—no bound states will occur. This can be made plausible by 
the following consideration. The energy of the perturbation will be of the 
order g®« and will be 40 only in a neighbourhood of dimensions 1/y of 
the nucleon. This energy has to be balanced against the energy required 
for the localization of a meson in this neighbourhood. This energy is n. It 
thus appears that the possibility of obtaining energy-values less than 
be only depends on the magnitude of g. The critical value of g represents 
the limit between weak and strong coupling theories. Restricting 
ourselves to weak coupling theories, we will expect no bound states 
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between mesons and nucleons and equation (15) will only be used to 
account for the distortion of (low energy) meson waves. 

This distortion may contribute to some of the well-known asymmetries 
observed in meson-production processes. 

When the self-consistent potential V has been determined up to the 
second approximation one arrives at a definite set of states N, . . . which 
one may call the set of states adapted to the 2nd order of perturbation 
theory. Each state will represent a definite pattern of motion—including 
that typical of scattering experiments—of bare particles interacting 
via the “ phenomenological” potential V. Some of the states N may 
have the asymptotic behaviour of closed states. The first-order matrix- 
elements in (9) will describe the absorption or emission of onemeson coupled 
to a nuclear transition MN. A special case of such a process has been 
‘treated by Gunn, Power and Touschek (1951). Their results show that 
the deformations caused by the assumed phenomenological internucleon 
potential may greatly affect the cross section for meson production. 
The high energy region of this process has received exhaustive treatment 
by C. Morette, whose results, however, cannot be extrapolated into the 
threshold region where—owing to the possibility of the formation of 
closed states between the nucleons—Born’s approximation is bound to 
break down. From the point of view of the present investigation the 
G.P.T.-treatment of the problem is still incomplete owing to the neglect of 
a possible distortion of the meson waves. 

The remaining terms of the second approximation to equation (6) will 
describe processes involving the emission and absorption of two mesons. 

-The third approximation will give further terms to the matrix-elements 
for the emission and absorption of one meson, in addition to three meson 
terms. It is clear that the terms of the first approximation together 
with those 1-meson terms arising from the third approximation should 
add up to C. Morette’s results in the range of validity of Born’s 
approximation. The terms arising in the third order will therefore 
describe the transition from the strong interaction low energy region 
treated in G.P.T. to the region of validity of Born’s approximation. 
The energy at which this transition is expected to occur will be of the 
order of the average well depth above the threshold. Recent experiments 
lie just in that region and therefore cannot be directly compared with any 
of the present calculations on the subject of meson production in proton- 
proton collisions. 

§ 3. 


So far we have only adapted the phenomenological potential V to the 
second approximation of perturbation theory. A possible method for 
the continuation of the process would be to leave V unaltered and 
continue the solution of equation (8) by orthodox methods and with a 
fixed potential and in terms of the stationary states determined in the 
second approximation. This treatment will be successful only if the 
higher even approximations do not lead to new forces sufficiently strong 
to produce bound states which have not already been accounted for in 
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the second approximation. The fourth-order formalism will introduce: 
interaction terms of the form 


i] 
: ‘aie 
a \f-=- (16) 
ne 1 


in addition to interaction (and self-energy) terms already treated in 
the second approximation. The first graph represents three nucleon 
forces, the second the interaction between a meson and two nucleons, 
the third the interaction between two mesons and the fourth the 
interaction between two mesons and a nucleon. The order of magnitude 
of the forces arising from the graphs (16) could be estimated by computing 
the matrix-elements corresponding to the graphs (16) in terms of the. 
second-order states N. If any of these graphs gives a contribution 
sufficient to lead to the formation of bound states a new supplementary 
potential V has to be introduced, capable of producing such states and 
fitted in such a manner that it compensates for all the even interaction 
terms including those of the fourth order. This means that the whole 
process described in §2 has to be repeated. As we have seen, this. 
retracing of the steps need not be compulsory, and it will indeed be 
unnecessary if the potential V adapted to the second approximation 
exhausts all the possible closed states of the system. The possibility 
of this situation is related to the short range of nuclear forces resulting 
from the non-vanishing mass of the meson. It is more likely in meson- 
theory than in quantum electrodynamics where the infinite range of 
nuclear forces entails an accumulation of bound states (level-density= 
infinite) at the boundary between the discrete and continuous spectrum. 

A possible objection to the suggestions outlined in this paper may be 
raised against its lack of relativistic invariance. Though, of course, the 
total results obtained in the case of a convergent procedure should be 
relativistically invariant this is not true for the individual contributions 
SM, S®, ... to the S-matrix. The reason for this is the obviously 
non-relativistic character of bound states between bare particles. It 
appears, however, that this lack of invariance is due to the physical 
interpretation of scattering experiments in terms of such bound states, 
which with the introduction of V is imposed upon the formalism. An 
investigation of the relativistic properties of the method will be attempted. 
elsewhere. 
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CXIX. CORRESPONDENCE. 


On Loading Nuclear Emulsions with Wires. 


By M. Danysz and G. Yexurtet, 
H. H. Wills Physical Laboratory, University of Bristol*. 


[Received July 17, 1951.] 


A MEHTOD of loading nuclear emulsion with wires, similar to that recently 
published by Meulmans, Occhialini and Vincent (1951), has been developed 
independently in this laboratory. Nylon threads and copper wires of 
thickness 25, 40 and 60 microns have been used. In general, for substances 
which do not react with the emulsion, there is no difficulty in developing 
the plates. For substances which do react, such as copper, we have found 
that satisfactory protection of the emulsion can be given by coating the 
wires with a thin layer of enamel. 

The most important technical problem met in processing such plates 
is the elimination of the distortion introduced by the presence of the 
wires. This arises from volume changes of the emulsion, and movements 


Fig. 1. 
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of the wires relative to it, during processing. The degree of distortion 
increases rapidly with the thickness of the wires employed, and it arises 
from twe principal causes :— 

(1) “Longitudinal” distortion is introduced by a lengthwise displacement 
of the wires with respect to the emulsion. The wires, when first placed 
in the emulsion, are under a slight tension, and it is important that this 
should be released in the early drying stage, so that it is not “frozen in”’. 
Should any tension remain when the plates are dry, it will be released 
during processing, and thus cause a displacement both of the emulsion 
and of the tracks in the immediate vicinity of the wires. 

(2) “Volume” distortion is caused by the shrinkage of the emulsion 
round the wires after processing, and is the more serious. 

To eliminate both longitudinal and volume distortion, the following 
methods of loading and threading the emulsion were adopted. The wire 
or thread is first evenly wound on a lathe, round a prepared frame (aa), 
fig. 1, to which it is stuck with shellac. The frame itself is made of 
glass 1 mm. thick and 1 cm. wide, and half the threads are cut away, 
only those attached to one face of the frame being retained. 
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The plates upon which the emulsion is poured are of specially treated 
glass supplied by Messrs. Ilford Ltd. On such a plate is placed, first, 
another glass frame (bb), with inside area equal to that of the desired 
emulsion surface, the upper and lower surfaces of (bb) having been first 
smeared with vaseline. The larger frame (aa) is now placed over (bb) 
so that the wires rest on its upper surface without bending or extra 
tension. A third frame, (cc), exactly similar to (bb), and of which the 
surfaces have also been smeared with vaseline, is placed over the top 
of (bb). 

The emulsion is now poured on the glass plate—within the area 
defined by the frames (bb) and (cc)—until its upper surface is just level 
with the top surface of (cc). When the emulsion gels, usually after 1} 
hours, the wires are cut away from the frame with a sharp razor blade, 
and the frames removed. The exact stage at which this is done is critical, 
for the emulsion must be sufficiently solidified to be able to support the 
wires, but not yet so rigid that it “freezesin”’ the tensions. If the emulsion 
is now left to dry it shrinks to a thickness of approximately 300 microns. 

The most serious distortion, volume distortion, arises because the 
gelatine shrinks as the result of the removal of the large quantities of 
silver halide during processing, whereas the wires are unaffected. The 
only way to avoid it is to prevent the shrinkage of the emulsion, and we 
have succeeded in doing so by the following method :— 

In effect, we replace the silver salt by a suitable transparent substance, 
and for this purpose we introduce resin into the emulsion, a method which 
has been developed in collaboration with Mr. J. E. Hooper. After 
washing, the water in the emulsion is replaced by alcohol by immersion in 
a succession of baths of gradually increasing concentration of alcohol. 
The plates are then soaked in a solution of resin in alcohol for a time of 
the order of two days, and then allowed to dry. 

The final thickness of the resulting emulsion appears to depend on the 
strength of the resin solution used and the period of soaking. Concentra- 
tions of 30 gm. of resin in 100 c.c. of alcohol, and soaking for two days 
give a final emulsion thickness close to that before processing ; but even 
greater thickness can be achieved, so that the emulsion is expanded 
instead of contracted. 

Nylon loaded plates exposed to cosmic radiation and treated by these 
methods have been examined, and the results show that both types of 
distortion had been almost completely eliminated. We believe that, 
apart from its particular applications to the present problem, the impreg- 
nation of emulsions with resins and other substances has other more 


general advantages and applications. These will be described in a later 
communication. 
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Neutron Emission from Nuclei Excited by High Energy Protons. 


By D. M. Sxyrng, 
Atomic Energy Research Establishment, Harwell, 
and 
W.S. C. WiLtiaMs, 


University College, London*. 
[Received August 1, 1951.] 


NucueaR reactions induced by high energy protons are thought to take 
_ place in two stages (Serber 1947). In the first stage the incoming particle 
makes collisions with individual nucleons in the target nucleus, and some 
of the products of these collisions pass directly out of the nuclear field. 
These “ prompt ” particles have high kinetic energy and they are emitted 
predominantly in the direction of motion of the incoming particle. After 
their departure the residual nucleus is left highly excited, and in the 
second stage of the reaction the excitation energy is dissipated by the 
emission of neutrons, charged particles and eventually y-rays. These 
“evaporation ”’ particles have relatively low energies, and they are 
emitted isotropically because of the comparatively long lifetime of the 
intermediate nuclei. Evaporation protons have often been studied 
before (e.g. Harding, Lattimore and Perkins 1949), but there has been 
little investigation of the neutrons following high energy reactions. 
This letter reports preliminary results of an investigation of evaporation 
neutrons from carbon and tungsten targets bombarded in the Harwell 
cyclotron. Values have been obtained for the differential cross-section 
for neutron evaporation as a function of the kinetic energy of the emitted 
neutron. If some assumptions are made as to the value of the proton 
inelastic cross-section of the target nuclei, the average number of neutrons 
emitted per excited nucleus can be determined also. 

The two targets were exposed in turn to the internal proton beam, 
whose maximum energy was calculated to be 171 MeV. Slowing down 
in the targets reduces the effective primary proton energy to about 
157 MeV. in both cases. Neutrons emitted backwards (7. e. at 180° to 
the direction of the ingoing proton beam) were collimated by 2 in. diameter 
holes in concrete blocks and entered photographic plates 18-9 metres 
from the target. The distributions in energy of these neutrons were 
determined by means of the tracks of recoil protons produced in the 
photographic emulsions (Champion and Powell 1944). The numbers 
of primary protons striking the targets in the two runs were determined 
measuring the amount of “Be produced by the “C(p, 3p3n)‘Be reaction, 
whose cross-section is known with an accuracy of 10 per cent (Dickson 
and Randle 1951). Several thin plates of carbon were incorporated in 
the tungsten target for this purpose. Two subsidiary experiments were 
carried out to check that the neutrons did, in fact, come directly from 
the target. Negligible numbers of recoil proton tracks were observed 
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in the plates (a) when the target was out of line with the peace 
holes, (b) when the target was realigned and the collimating holes w 
<ed with copper rods. ; 

eee re cross-section, o, for neutron evaporation from Vs 
tungsten target is shown in fig. 1. The value given for ee wit 
energy less than 0:5 MeV. may be as much as 30 per cent low ( eae 
and Reines 1950). A change of this magnitude will not, however, 
significantly affect the integrated cross-section. The total cross-section, 
for neutron energies up to 11 MeV. is found to be (6-010°9) barns 
(1 barn=10-24 cm.2), and the mean energy of an evaporation ae 
is 2-6 MeV. There are 675 tracks in the spectrum below 11 MeV. ; only 
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two tracks corresponding to neutron energies greater than this were- 
observed. The shape of this spectrum agrees fairly well with that to be 
expected from evaporation theory (Weisskopf 1937), if account is taken 
of the fact that more than one neutron can be evaporated from each 
excited nucleus. ; 

In order to calculate the average number of neutrons evaporated 
from an excited nucleus, it is necessary to know the cross-section for 
inelastic collisions of protons with the target nuclei. Such values are 
not available, but it is known that the total cross-section of carbon for 
156 MeV. neutrons is (0-330 +0-003) barns (Taylor, Pickavance, Cassels 
and Randle 1951 a), and that of tungsten for neutrons of the same energy 
is 3-1 barns (Taylor, Pickavance, Cassels and Randle 1951 b). If it is 
assumed that one-half the neutron total cross-section is inelastic, and that 
at these energies the inelastic neutron and proton cross-sections are equal, . 
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the number of neutrons evaporated per excited nucleus in the tungsten 
target is equal to 4-0. If the binding energy of the most loosely bound 
neutron in a nucleus in the neighbourhood of tungsten is taken to be 
6 MeV., the average energy lost in neutron evaporation by a nucleus 
in the tungsten target is 34-4 MeV. An estimate of the mean excitation 
energy of the evaporating nucleus may be obtained if it is assumed that 
one proton with an average kinetic energy of 10 MeV. is evaporated 
per nucleus, and that 7 MeV. are lost in other processes, including 
y-emission (Le Couteur 1950). If the binding energy of the evaporated 
proton is taken as 7 MeV., the value (34--17+7)=—58 MeV. is obtained 


Fig. 2. 
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for initial excitation energy—i. e. 0-37 of the energy of the primary 
protons. Goldberger (1948) calculates that the average excitation energy 
of a heavy nucleus when bombarded with 90 MeV. neutrons is 42-5 MeV. 
It is hoped to investigate, by a method similar to that described here, 
the energy spectrum of the evaporated protons, and the cross-section 
for this process, in order to obtain more precise information on the 
average energy of the excited nucleus. 

Results obtained with the carbon target are shown in fig. 2. 
Curve A shows the differential cross-section for neutron emission, assuming 
that the evaporating nucleus was at rest in laboratory space. In fact, 
on account of its small mass, it may have been recoiling with sufficient 
velocity to affect the observed spectrum. It is difficult to estimate an 
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average value for the recoil velocity, but this will certainly be much less. 
than that corresponding to the full transfer of momentum from the 
incident proton to the nucleus as a whole. Curve C, in fig. 2, is derived 
on the basis of full momentum transfer, and is therefore greatly over- 
corrected for recoil. Curve B, calculated on the assumption of 30 per cent 
momentum transfer, probably represents a more reasonable corrected 
spectrum. The presence of a relatively large number of high energy 
neutrons makes it difficult to give unambiguous values for the total 
“evaporation ” cross-section and for the mean neutron energy. The 
figures given below were obtained by cutting off the spectra at 14 MeV. 


Curve A | Curve B | Curve C 


‘ 


Total cross-section for “‘ evaporation ” of : 
neutrons (barns) 0-027 0-038 0-068 


Mean kinetic energy of neutrons (MeV.) 4-5 4-6 5-4 


Average number of neutrons per excited 
nucleus 0-16 0-23 0-41 


The general conclusions to be drawn from these data are fortunately 
not affected by the considerable uncertainties in the corrections which 
have been applied to the observed spectrum. The average number of 
neutrons “ evaporated ”’ is less than one per excited nucleus, and there are 
relatively more high energy neutrons than would be expected from 
a nucleus in thermodynamic equilibrium. 
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The Photodisintegration of °B into Neutron, Proton and two «-Particles. 


By M. J. Brinkworts and KE. W. Tirrerton*, 
Atomic Energy Research Establishment, Harwellt. 


[Received July 27, 1951.] 


In the investigation of the reactions B(y, t)2 4He (Titterton and Calcraft 
1951) and 1°B(y, d)2 4He (Brinkworth and Titterton 1951) with the 
*Li(p, y) radiation of 14:8 and 17:6 MeV. energy, three-pronged stars 
have been observed in Ilford E, boron loaded emulsions, which fail 
to give a momentum balance for either of these reactions, or for the 
reaction C(y, 3x). The impossibility of a momentum balance is 
immediately obvious in some cases, when all three tracks fall within an 
angle less than 180°. It appears that at least one neutron is emitted in 
each of these events. The disintegrations are consistent with the reaction 


0B hy>2 *Het1H-+—8-07MeV. . . . . (1) 


which was reported by Goward, Titterton and Wilkins (1950). 

Sometimes the proton track may be distinguished from the «-particle 
tracks by its lower grain density, but in many cases this is impossible, due 
to the intensifying action at development of the boron loading. In such 
cases, it is only when one track is too long to be an «-particle, from 
energy considerations, that the proton can be identified. Consequently, 
the proton energy distribution is known only for energies greater than 
about 2 MeV. 

Of the 182 observed events which have been ascribed to reaction (1), 
101 are such that the proton can be identified. In these cases the neutron 
energy can be deduced by application of the conservation of momentum, 
but the probable errors in the observed energies of the proton and 
a-particles are too great for this to be a satisfactory method. In 
consequence, only those events with K,—E,>15 MeV. can definitely be: 
assigned to the 17-6 MeV. y-ray ; those with E,—E,<15 MeV. may be 
due to either the 14:8 or 17-6 MeV. y-rays. It is therefore difficult to 
obtain the variation of cross-section with energy, but the experiment 
indicates that it does not fall off rapidly from 14-8 to 17-6 MeV. as does 
the B(y, d)2 4He cross-section (Brinkworth and Titterton 1951). This 
rapid fall of the (y, d) cross-section may be due to the competition of 
reaction (1) with increasing energy. The mean cross-section of reaction (1) 
for the 7Li(p, y) radiation is estimated to be (3-2261-3) x 10° 78:cm.3: 

The selection of events consequent on the inability to identify low 
energy protons in the present experiments renders it difficult to decide 
on the course of the reaction. 
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Neglecting a three-body reaction, the possible modes of disintegration 
which need to be considered are :— : 

(i.) Initial (y, 2) reaction, the residual excited °B decaying through 
either 8Be* or ®Li*. If the reaction proceeded through the ground state 
of °B, the event would be too small to be observed. 

(ii.) Initial (y, p) reaction, the residual *Be* decaying through either 
8Be* or ®He*. The level of 9Be at 2:42 MeV. (Hornyak, Lauritsen, 
Morrison and Fowler 1950) would result in a proton energy of 5:3 or 
7-8 MeV. for the 14:8 and 17-6 MeV. y-rays respectively. Two events 
have been obtained with protons of the latter energy. The only three 
events with proton energy of approximately 5-3 MeV. have each too 
great a total energy to be caused by the 14:8 MeV. y-ray. 

(iii.) Initial (y,«) reaction, leaving ®Li*. Although the reaction 
10B(y, d)2 *He does not appear to proceed through *Li (Brinkworth and 
Titterton 1951) the reaction (1) may well do so, since it is known that 
the reaction Li®+hv+4*He+1H-+ 1m does occur (Titterton and Brinkley 
1951), 

(iv.) Emission of a deuteron in the singlet state, dissociating into 
neutron and proton, since neutron scattering in ortho- and para-hydrogen 
shows that the singlet state is virtual. The residual *Be* would yield 
two «-particles. 

An attempt to decide between these alternatives can be made only 
when it is possible to identify the protons in all events so that complete 
distributions of particle energies or Q values may be examined. 

Members of the nuclear emulsions group at the Atomic Energy Research 
Establishment are thanked for microscope work. 
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